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ABSTRACT 

This  experimental  study  was  conducted  at  the  Naval  Post- 
graduate School  to  investigate  the  heat  transfer  character- 
istics of  flat  vertical  plates  in  a  rectangular  gas  - f luidized 
bed.   The  primary  objective  was  to  determine  what  effect 
variations  in  the  bed  width-to-height  ratio  had  on  heat  trans- 
fer to  the  vertical  flat  plates  forming  the  container  walls. 
The  experiment  was  conducted  using  a  specially  heated  and  in- 
strumented fluidized  bed  equipped  with  a  movable  side  wall 
which  permitted  modification  of  the  bed  geometry.   As  the 
width  of  the  bed  was  adjusted,  the  settled  bed  height  was  main- 
tained at  a  constant  level  by  the  addition  or  removal  of  bed 
material . 

A  secondary  objective  of  the  study  was  to  determine  the 
effect  of  variations  in  the  fluidization  gas  flow  rate  on  heat 
transfer  to  the  bed  walls.   Flow  rates  ranging  from  fixed  bed 
fluidization  to  pneumatic  conveying  were  studied. 

Pressure  drop  measurements  as  well  as  visual  observations 
were  used  to  determine  minimum  fluidization  flow  rates  for 
each  configuration.   Heat  transfer  coefficients  were  calculated 
for  each  flow  rate  and  bed  geometry  using  temperature  data  ob- 
tained from  a  computer  controlled  thermocouple  network.   In 
addition,  a  study  was  made  of  variations  in  wall  temperature 
with  changes  in  vertical  position. 
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A       Ammeter  -  Figure  19 

2 
A,       Heater  Wall  Bottom  Section  Surface  Area  (in  ) 

? 
A       Copper  Plate  Surface  Area  (in") 

A,      Distributor  Surface  Area  (in") 

2 
A,-      Front/Rear  Wall  Surface  Area  (in  ) 
fr 

2 
A.       Heater  Insulation  Surface  Area  (in") 
l 

2 
A  ,      Heater  Mounting  Block  Surface  Area  (in  ) 

2 
A       Heater  Wall  Side  Section  Surface  Area  (in  ) 
s 

2 
A       Heater  Wall  Top  Section  Surface  (in  ) 

C        Specific  Heat  (Btu/Lbm-°F) 

D       Particle  Diameter  (in) 

P 

E       Heater  Voltage  (Volts) 

2 
G       Superficial  Gas  Mass  Velocity  (Lbm/Hr-Ft  ) 

g       Proportionality  Constant  in  Newton's  Second  Law 
c      of  Motion  (Lbm-Ft/Lbf-Sec2) 

G  r  Superficial  Mass  Velocity  at  Minimum  Fluidization 

mr      (Lbm/Hr-Ft2) 

2 
G       Optimum  Superficial  Mass  Velocity  (Lbm/Hr-Ft  ) 

H       Fluidized  Bed  Height  (in) 

h       Heat  Transfer  Coefficient  (Btu/Hr-Ft 2 -°F) 

h  ,      Heat  Transfer  Coefficient  Copper  Plate  to  Bed 
CD      (Btu/Hr-FtZ  -°F) 

•hmoV    Maximum  Heat  Transfer  Coefficient  (Btu/Hr-Ft   -°F) 
max 
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hmba  Heat  Transfer  Coefficient  Heater  Mounting  Block 
to  Atmosphere  (Btu/Hr-Ft2  -°F) 

Hg  Expanded  Bed  Height  (in) 

I  Heater  Current  (Amps) 

kc  Thermal  Conductivity  of  Copper  Plate  (Btu/Hr- Ft -°F) 

kfr  Thermal  Conductivity  of  Front/Rear  Walls  (Btu/Hr- Ft  -  °F) 

k-  Thermal  Conductivity  of  Heater  Insulation 
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L^  Thickness  of  Heater  Insulation  (in) 

L  ^  Thickness  of  Heater  Mounting  Block  (in) 

m  Air  Mass  Flow  Rate  (Lbm/Hr) 

P  Atmospheric  Pressure  (in.  Hg  absolute) 

a 

q  Total  Thermal  Energy  (Heat  Flux)   Losses  to 
Atmosphere  (Btu/Hr) 

q  Thermal  Energy  Carried  Out  of  Bed  by  Outgoing 


Lao 


Airstream  (Btu/Hr) 


q.       Total  Heat  Flux  Entering  Bed  from  All  Sources 
D       (Btu/Hr)  ' 

q.       Thermal  Energy  Lost  to  Atmosphere  from  Wall  Section 
Below  Heater  (Btu/Hr) 

q,.      Thermal  Energy  Entering  Bed  from  Wall  Section  Below 
Heater  (Btu/Hr) 


Total  Heat  Flux  Entering  Bed  from  All  Right  Hand 
Sources  (Btu/Hr) 


^br 

q,  -,      Total  Heat  Flux  Entering  Bed  from  All  Left  Hand 


Sources  (Btu/Hr) 


q  ,      Heat  Flux  into  Bed  from  Copper  Plate  (Btu/Hr) 

q-      Thermal  Energy  Lost  to  Atmosphere  through  Front/ 
±ra    Rear  Walls  (Btu/Hr) 

q.       Electrical  Energy  into  Heater  (Btu/Hr) 

q       Total  Heat  Flux  Leaving  Bed  (Btu/Hr) 

q  ,      Total  Heat  Flux  into  Bed  from  Sources  Other 
than  the  Copper  Plate  (Btu/Hr) 

q  ,      Thermal  Energy  Lost  through  Heater  Mounting 
mD      Block  (Btu/Hr) 

q       Thermal  Energy  Lost  to  Atmosphere  from  Wall  Sections 
sa     Adjacent  to  Heater  (Btu/Hr) 

q  ,      Thermal  Energy  Entering  Bed  from  Wall  Sections 
S       Adjacent  to  Heater  (Btu/Hr) 

q       Thermal  Energy  Lost  to  Atmosphere  from  Wall  Section 
Above  Heater  (Btu/Hr) 

q..      Thermal  Energy  Entering  Bed  from  Wall  Section  Above 
tD     Heater  (Btu/Hr) 

R(N)     Average  Heater  Row  Temperature  (°F)  (N  indicates 
row  number) 

R       Particle  Reynolds  Number 
e  J 

R  _r    Particle  Reynolds  Number  at  Minimum  Fluidization 
em±  J 

T(N)     Thermocouple  Temperature  Reading  (°F)  (N  indicates 
thermocouple  number) 

T  Ambient  Temperature  (°F) 

T  .  Temperature  of  Air  Entering  Bed  (°F) 

T  Temperature  of  Air  Leaving  Bed  (°F) 

T,  Average  Bed  Temperature  (°F) 

T,  Average  Heater  Temperature  (  F) 

T  ,      Average  Temperature  of  Heater  Mounting  Block 
Outer  Surface  (  F) 

T       Average  Temperature  of  Copper  Plate  Outer  Surface  (°F) 
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U  Superficial  Gas  Velocity  (Ft/Sec) 

V  Voltmeter  -  Figure  19 

V  Air  Volumetric  Flow  Rate  (Ft3/Hr) 
W  Bed  Width  (in.  ) 

X  Highest  Number  Thermocouple  Probe 

GREEK  LETTER  SYMBOLS 

AP  Pressure  Drop 

AP  r  Pressure  Drop  at  Minimum  Fluidization 

m±  r 

e  Particle  Voidage 

e  r  Particle  Voidage  at  Minimum  Fluidization 

m±  ° 

y  Dynamic  Viscosity  (Lbf-Sec/Ft  ) 

p  Air  Density  (Lbm/Ft  ) 

3. 

p  Particle  Density  (Lbm/Ft  ) 
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I.   INTRODUCTION 

A  fluidization  bed  in  its  simplest  form  consists  of  a  con- 
tainer filled  with  small  particles  which  are  subjected  to  an 
upward  flow  of  fluid  emerging  from  the  bottom  of  the  container. 
When  the  upward  viscous  force  (drag)  exerted  on  the  particles 
by  the  moving  fluid  balances  or  exceeds  the  downward  force 
of  gravity,  the  particles  become  suspended  and  begin  to  move 
about  the  container  in  a  fluid-like  manner,  hence  the  term 
"fluidization."   This  churning  mixture  of  particles  and  gas 
is  often  referred  to  as  the  "emulsion  phase."   When  the 
velocity  of  the  fluid  is  increased  to  the  point  that  the 
particles  are  carried  out  of  the  container,  a  transition  is 
made  to  a  different  form  of  fluidization  known  variously  as 
pneumatic  conveying,  entrained  flow  fluidization,  or  carryover. 
When  the  velocity  of  the  fluid  is  too  low  to  counteract  the 
gravitational  force  and  cause  fluidization,  the  bed  is  referred 
to  as  being  "packed"  or  "fixed"  and  the  process  is  known  as 
flow  through  a  porous  medium  or  fixed  bed  fluidization. 

Distinctly  different  modes  of  behavior  are  observed  depend- 
ing on  whether  the  fluid  passing  through  the  container  is  a 
liquid  or  a  gas.   Liquid  fluidized  beds  are  often  referred  to 
as  homogeneous  or  particulately  fluidized  beds.   This  type  of 
bed  tends  to  display  a  very  smooth  transition  through  the 
various  stages  of  fluidization  from  fixed  to  flui'dized  to 
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carryover.   Concurrently  a  very  smooth  transition  in  heat 
transfer  characteristics  and  hydrodynamic  behavior  is  observed. 
On  the  other  hand,  gas  fluidized  beds  tend  to  display  a  more 
erratic  performance.   Most  researchers  attribute  the  difference 
in  behavior  to  the  differing  particle  to  fluid  density  ratios. 
Gas  fluidized  beds  are  known  as  "aggregate"  or  "bubbling 
fluidized  beds."   In  this  type  of  bed,  once  fluidization  has 
begun,  the  gas  usually  forms  into  bubbles  which  displace  the 
bed  particles  as  the  bubbles  rise.   The  appearance  of  a  bubbling 
bed  is  remarkably  similar  to  that  of  a  boiling  liquid,  as  shown 
in  Figure  1.   During  conditions  of  very  active  fluidization 
the  bubbles  sometimes  form  together  or  coalesce  to  create  a 
large  bubble  whose  diameter  approaches  the  width  of  the  bed, 
as  shown  in  Figure  2.   This  condition  is  known  as  "slugging" 
or  "heaving"  and  can  significantly  effect  bed  performance. 
Another  phenomenon  of  gas  fluidized  beds  is  shown  in  Figure  3. 
It  is  a  condition  known  as  "channelling."   When  this  occurs, 
the  gas  travels  through  a  fairly  confined  area  forming  a  gas 
filled  tube  within  the  bed  rather  than  passing  evenly  through- 
out the  cross  section.   This  behavior  is  believed  to  be 
caused  by  poor  distributor  design  and,  like  slugging,  can 
strongly  effect  the  heat  transfer  characteristics  of  the  bed. 

Fluidized  beds  offer  a  wide  variety  of  beneficial  char- 
acteristics such  as  high  rates  of  heat  transfer,  high  thermal 
inertia,  large  solid  to  fluid  contact  area,  and  isothermal 
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bed  temperatures.   In  addition,  once  fluidized,  the  pressure 
drop  across  a  bed  remains  relatively  constant  despite  varia- 
tions in  fluid  flow  rate.   The  pressure  drop  is  a  function 
of  bed  height  and  particle  density  and  is  independent  of 
particle  size.  [Ref.  1]   This  permits  increasing  the  surface 
area  by  using  smaller  particles  without  paying  a  penalty  in 
increased  pressure  loss. 

Fluidized  beds  also  suffer  from  certain  limitations.   The 
size  of  particles  which  can  be  fluidized  is  limited  to  a 
specific  range  extending  from  approximately  0.001  inch  in 
diameter  to  about  0.25  inch  diameter,  depending  on  the  density 
of  the  particles  and  the  fluid.  [Ref.  2]   In  addition,  some 
materials  cannot  be  fluidized  regardless  of  size  because  of 
their  tendency  to  clump  together  and  form  particles  too  large 
to  fluidize,  or  because  of  their  tendency  to  break  down  into 
particles  which  are  so  small  that  they  become  entrained  in  the 
fluidizing  medium  and  are  carried  from  the  bed.   Another  dis- 
advantage is  that  fluidization  requires  the  expenditure  of 
energy  to  force  the  fluidizing  medium  through  the  bed  at  the 
required  velocity. 

In  general  the  beneficial  aspects  of  fluidized  beds  far 
exceed  the  limitations  and  as  a  result,  fluidized  bed  prin- 
ciples have  been  utilized  in  many  diverse  applications 
over  the  years.   Early  mining  engineers  used  fluidization 
to  separate  solids  of  different  densities  by  suspending 
them  in  rising  streams  of  water.  [Ref.  1]   Prior  to  World 
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War  II,  chemical  engineers  discovered  that  fluidized  beds 
made  outstanding  chemical  reactors  because  of  the  extremely- 
large  contact  area  between  the  particles  and  the  fluid. 
For  example,  a  container  measuring  one  foot  by  one  foot  by- 
one  foot,  filled  with  0.01  inch  diameter  spherical  particles 
yields  approximately  4,000  squares  feet  of  surface  area. 
This  characteristic  of  fluidized  beds  has  been  utilized  to 
great  advantage  in  a  wide  number  of  chemical  processes,  most 
notable  the  catalytic  cracking  of  petroleum.  [Ref.  3] 
Because  the  particles  within  a  fluidized  bed  are  constantly- 
churning  about  and  mixing,  the  temperature  at  one  point  in 
the  bed  is  usually  in  very  close  agreement  with  the  tempera- 
ture at  any  other  point.   Industry  has  capitalized  on  this 
isothermal  property  by  using  fluidized  beds  to  improve  the 
heat  treatment  of  metals  in  such  processes  as  quenching, 
tempering,  and  carburizing.   As  a  combined  result  of  the 
large  surface  area  and  the  thorough  mixing  of  the  particles, 
fluidized  beds  possess  impressive  heat  transfer  capabilities. 
Consequently,  the  power  generation  industry  has  begun  utiliz- 
ing this  technology  to  create  a  new  generation  of  heat 
transfer  devices  such  as  steam  generators  which  efficiently 
burn  poor  quality  fuels  while  simultaneously  yielding  very 
low  levels  of  harmful  emissions.   These  devices  are  capable 
of  burning  high  sulfur  coal,  pelletized  garbage,  and  even 
raw  crushed  shale,  just  as  it  comes  from  the  ground.   Another 
application  is  seen  in  waste  heat  recovery  units  which  are 
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fluidized  by  the  exhaust  gases  from  conventional  boilers  or 
gas  turbines.   The  thermal  energy  which  would  otherwise 
have  been  wasted  is  transferred  to  a  secondary  fluid  cir- 
culating through  tubes  emersed  within  the  bed.   This  re- 
covered energy  is  then  utilized  to  increase  overall  plant 
efficiency.  [Ref.  4] 

The  use  of  fluidized  bed  technology  has  expanded  greatly 
in  recent  years  as  hundreds  of  new  applications  have  been 
discovered,  ranging  from  advanced  ship  propulsion  plants  to 
sophisticated  solar  energy  conversion  and  storage  devices. 
As  the  utilization  of  these  devices  increases  so  does  the 
quest  for  a  better  understanding  of  the  many  complex  phenom- 
enon associated  with  f luidization . 

Fluidized  bed  research  has  been  conducted  throughout  the 
world  and  more  has  been  written  about  fluidization  than  can 
be  read  in  a  lifetime,  however  much  remains  to  be  learned 
in  order  to  unlock  the  full  potential  of  this  fascinating 
process.   This  thesis  represents  a  modest  attempt  to  add  to 
the  knowledge  of  fluidization  in  the  area  of  heat  transfer 
to  vertical  flat  plates  in  gas  fluidized  beds. 
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II.   OBJECTIVES  AND  METHOD  OF  APPROACH 

A.  GENERAL  STATEMENT  OF  THE  PROBLEM 

The  bulk  of  fluidized  bed  heat  transfer  research  conducted 
to  date  has  been  related  to  the  study  of  heat  transfer  to  or 
from  objects  immersed  within  the  bed,  particularly  tubes  and 
tube  bundles.   Comparatively  little  research  has  been  done 
on  the  study  of  heat  transfer  to  the  walls  of  the  containing 
vessel  itself.   In  addition,  a  large  majority  of  the  test 
apparatus  utilized  thus  far  to  collect  experimental  data 
have  been  cylindrical  in  configuration.   Very  few  studies 
have  been  conducted  on  heat  transfer  within  rectangular 
fluidized  beds;  a  geometry  which  is  being  used  with  increasing 
frequency  in  a  wide  variety  of  applications. 

The  principle  objective  of  this  thesis  was  to  gain  a 
greater  understanding  of  the  heat  transfer  to  the  flat  ver- 
tical plates  which  form  the  walls  of  a  rectangular  gas 
fluidized  bed  and  to  determine  what  effects,  if  any,  varia- 
tions in  bed  geometry  had  on  heat  transfer  performance. 

A  secondary  objective  was  to  study  the  variations  in  heat 
transfer  which  occurred  as  the  fluidization  gas  flow  rate 
was  adjusted  over  a  wide  range. 

B.  METHOD  OF  APPROACH 

This  investigation  was  primarily  experimental  in  nature 
with  limited  utilization  of  analytical  methods  for  comparison 
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purposes.   To  acquire  the  necessary  data  for  this  study,  a 
computerized  data  acquisition  system  was  linked  to  a  network 
of  up  to  75  copper-constantan  thermocouples  installed  in  a 
specially  instrumented  and  heated  fluidized  bed.   This  system 
was  used  to  obtain  a  thermal  "snapshot"  of  the  test  apparatus 
which  showed  the  temperature  at  numerous  key  locations  at 
a  specific  point  in  time.   This  data  was  then  combined  with 
measured  values  of  input  power  and  incoming  air  flow  rates 
to  produce  a  calculated  heat  transfer  coefficient.   Experi- 
mental data  was  obtained  for  a  number  of  different  air  flow 
rates  and  bed  conditions  beginning  with  fixed  bed  fluidiza- 
tion  and  progressing  through  to  the  onset  of  entrained  flow 
f luidization .   Once  data  had  been  collected  for  the  entire 
range  of  flow  rates  ,  the  bed  geometry  was  altered  by  moving 
an  adjustable  side  wall  within  the  test  apparatus  so  that 
the  width  of  the  rectangular  bed  was  increased.   Additional 
particles  were  then  added  to  the  bed  so  that  the  bed  depth 
remained  at  the  same  level  as  in  the  previous  configuration. 

The  data  collection  process  was  then  repeated  in  its 
entirety  for  the  new  configuration.   A  total  of  four  differ- 
ent bed  geometries  were  studied  in  the  course  of  this  work. 
Over  5,700  temperature  readings  were  recorded  for  eighty  dif- 
ferent test  runs.   All  experimental  data,  including  narrative 
comments  regarding  observed  bed  behavior  were  recorded  on 
magnetic  disks.   The  recorded  data  was  automatically  pro- 
cessed by  computer  and  then  printed  out  in  such  a  manner 
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that  bed  performance  for  a  particular  test  run  was  readily 
apparent.   The  calculated  results  for  the  individual  runs 
were  then  analyzed  to  determine  overall  trends. 

In  addition  to  the  heat  transfer  investigation,  an  ex- 
perimental study  of  bed  pressure  drop  characteristics  was 
conducted  for  each  geometric  configuration  in  order  to 
establish  the  point  of  fluidization  onset  and  to  confirm 
proper  operation  of  the  fluidizing  apparatus. 
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Ill .   EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

A.   DESCRIPTION  OF  FLUIDIZATION  APPARATUS 

The  f luidization  apparatus  consisted  of  a  rectangular 
plexiglas  box  filled  with  0.012  inch  diameter  silica  glass 
beads  as  shown  in  Figure  4.   The  box  was  equipped  with  a 
movable  side  wall  which  permitted  the  inside  dimensions  of 
the  bed  to  be  varied  from  6.25  inches  wide  by  6  inches  deep 
by  18  inches  high  to  24  inches  wide  by  6  inches  deep  by  18 
inches  high.   The  movable  wall  was  equipped  with  felt  gaskets 
which  prevented  the  escape  of  beads  through  the  gap  between 
the  movable  side  wall  and  the  fixed  front  and  rear  walls. 
0.5  inch  thick  clear  plexiglas   (grade  GM)  was  used  through- 
out the  construction  and  all  components  were  assembled  using 
plastic  laminating  cement  with  wood  screw  reinforcements. 
The  floor  of  the  apparatus  also  served  as  the  distributor 
plate  whose  function  was  to  evenly  distribute  the  flow  of 
incoming  air  over  the  entire  cross  section  of  the  bed.   The 
distributor  plate  was  constructed  from  a  plexiglas  sheet 
with  0.125  inch  diameter  holes  drilled  0.25  inches  on  center. 
Secured  over  this  plate  was  a  140  mesh  stainless  steel  wire 
screen  to  prevent  the  beads  from  falling  down  into  the  holes 
and  to  provide  an  even  finer  distribution  of  air. 

Located  beneath  the  distributor  plate  was  an  air  chamber 
or  plenum  which  served  to  ensure  that  the  air  pressure  on 
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the  underside  of  the  distributor  plate  was  the  same  at  any 
point.   This  was  accomplished  by  first  directing  the  incoming 
air  through  a  2.5  inch  inside  diameter  plexiglas  distribution 
tube  centered  lengthwise  within  the  24  inch  by  6  inch  by  6 
inch  plenum.   This  tube  had  0.5  inch  diameter  holes  drilled 
on  1  inch  centers  on  its  underside  and  served  to  evenly  direct 
the  incoming  air  downward  over  the  entire  length  of  the  plenum. 
In  order  to  prevent  the  escape  of  air  when  the  bed  was  con- 
figured to  a  width  less  than  the  24  inch  width  of  the  distri- 
butor plate,  the  movable  side  wall  was  equipped  with  a 
horizontal  baffle  plate  which  laid  flat  over  the  top  of  the 
unused  portion  of  the  distributor  plate.   Once  the  desired 
position  of  the  movable  wall  was  set,  the  baffle  plate  was 
securely  clamped  to  the  distributor  plate  and  the  edges  were 
sealed  with  silicon  caulking.   After  passing  through  the 
distributor  plate  and  finally  the  bed  itself,  the  fluidizing 
air  escaped  from  the  apparatus  via  one  or  more  of  three 
plexiglas  exhaust  tubes  positioned  on  the  top  of  the  apparatus. 
These  tubes  were  quipped  with  wire  mesh  covers  to  prevent  the 
loss  of  glass  beads  due  to  carryover. 

In  order  to  measure  the  air  pressure  at  various  depths 
the  apparatus  was  equipped  with  static  pressure  taps  posi- 
tioned along  the  rear  wall  with  2  inch  vertical  spacing 
between  centers.   Additionally  taps  were  installed  just 
below  the  distributor  plate  and  immediately  above  it.   These 
taps  were  connected  by  flexible  plastic  tubing  to  a  manifold 
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which  permitted  selective  connection  of  the  taps  to  a 
Meriam  Model  33KA35  30  inch  water  manometer  as  shown  in 
Figure  5. 

The  fluidizing  apparatus  was  supplied  with  air  from  a 
Spencer  axial  flow  turbo  compressor  rated  at  500  SCFM. 
The  fluidizing  air  was  not  filtered  or  dehydrated  prior 
to  entry.   The  volumetric  flow  rate  of  air  entering  the 
apparatus  was  measured  by  a  calibrated  Fisher  and  Porter 
flowmeter  rated  at  34  CFM. 

B.   DESCRIPTION  OF  HEAT  TRANSFER  APPARATUS 

The  test  apparatus  was  equipped  with  two  instrumented 
flat  plate  heater  assemblies  each  measuring  5  inches  wide 
by  10  inches  high.   These  assemblies  were  mounted  flush  with 
the  inside  surface  of  both  the  movable  side  wall  and  the 
fixed  side  wall  opposite  it.   Each  heater  assembly  was 
powered  by  a  Watlow  strip  heater  rated  at  250  watts.   The 
heater  consisted  of  a  grid  of  nickel  alloy  resistance  wire 
embedded  in  a  thin  sheet  of  rubberized  supporting  material 
and  backed  by  a  0.5  inch  thick  sheet  of  silicone  sponge 
rubber  insulation.   This  insulation  was  in  turn  bonded  to 
a  0.5  inch  thick  clear  plexiglas  mounting  block.   The  face 
of  the  assembly  consisted  of  a  0.25  inch  thick  copper  plate 
mounted  flush  with  the  heating  surface  of  the  strip  heater. 
In  order  to  avoid  the  use  of  screws  or  other  conventional 
attachment  methods  which  may  have  altered  or  impeded  the 
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heat  flow  from  the  strip  heater  to  the  copper  plate,  the 
two  were  held  together  using  a  strongback  arrangement  as 
shown  in  Figures  6  and  7.   This  mounting  method  ensured 
that  the  strip  heater  was  held  firmly  against  the  copper 
plate,  thus  helping  to  minimize  thermal  contact  resistance 
between  the  two.   To  further  reduce  this  resistance,  a  ther- 
mally conductive  paste  was  applied  to  the  heating  surface 
of  the  strip  heater  prior  to  positioning  the  copper  plate 
over  it.   Power  was  supplied  to  the  left  hand  heater  by  a 
Lambda  regulated  DC  power  supply  model  LK345A,  while  the 
right  hand  heater  was  supplied  by  a  Hewlett  Packard  DC 
power  supply  Model  6296A.   Both  power  supplies  were  rated 
at  60  Volts  and  2  Amps  and  were  each  equipped  with  an  in- 
dividual voltmeter  and  ammeter. 

All  thermocouples  used  in  this  apparatus  were  10  gauge, 
ANSI  Type  T,  copper-constantan ,  exposed  junction  type.   Each 
heater  assembly  was  equipped  with  twenty  thermocouples 
arranged  in  five  rows  of  four  thermocouples  each  as  shown 
in  Figures  8  and  9.   These  thermocouples  were  mounted  to  the 
copper  plate  in  such  a  manner  so  as  to  be  least  disruptive 
of  the  heat  distribution  patterns  within  the  plate.   The 
thermocouple  junctions  were  soldered  into  small  shallow  holes 
drilled  into  the  back  of  the  plates  so  that  the  top  of  the 
junction  was  flush  with  the  plate  surface.   The  thermocouple 
leads  were  then  routed  straight  back  through  holes  drilled 
into  the  strip  heater  and  its  mounting  block.   Figure  10 
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Figure  6.   Heater  Assembly  Mounting  Diagram;  Side  View 
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Heater  Assembly  Mounting  Diagram;  Rear  View 
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Figure  8.   Right  Hand  Heater  Thermocouple  Arrangement 
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Figure  9.   Left  Hand  Heater  Thermocouple  Arrangement 
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Figure  10.   Photograph  of  Heater  Thermocouple  Mounting  -Procedure 


37 


shows  the  thermocouple  mounting  procedure.   X-ray  photographs 
were  used  to  determine  the  exact  location  of  the  resistance 
wires  embedded  within  the  heaters  so  as  to  avoid  striking 
them  while  drilling  the  holes.   This  thermocouple  arrangement 
made  it  possible  to  read  the  plate  inside  surface  temperature 
immediately  adjacent  to  the  surface  of  the  strip  heater  while 
causing  minimum  disruption  of  the  heat  flux  patterns  within 
the  plate  itself.   A  determination  of  the  plate  outside 
surface  temperature  as  seen  by  the  bed  was  then  accomplished 
by  means  of  simple  calculations  as  shown  in  section  HIE. 

The  temperatures  within  the  bed  itself  were  monitored  by 
15  retractable,  24  inch,  stainless  steel  sheathed,  exposed 
junction,  Omega  Engineering  Type  T  thermocouple  probes.   The 
position  of  these  probes  was  held  fixed  at  a  bed  depth  of 
7  inches  above  the  distributor  plate  for  all  data  runs.   The 
number  of  bed  probes  actually  utilized  depended  on  the  geo- 
metric configuration  being  studied,  as  shown  in  Figure  11. 

Incoming  air  temperature  was  monitored  by  3  retractable, 
4  inch,  stainless  steel,  sheathed  thermocouple  probes  posi- 
tioned within  the  air  plenum  as  shown  in  Figure  12.   Out- 
going air  temperature  was  monitored  by  3  similar  probes 
mounted  one  each  in  the  center  of  the  3  air  exhaust  tubes , 
as  shown  in  Figure  13.    Because  of  the  positioning  of  these 
tubes  with  respect  to  the  movable  side  wall  only  one  thermo- 
couple was  used  to  read  the  air  outlet  temperature  for  all 
data  runs. 
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NOTE:   Arrangement  is  as  viewed  from  top  front 


Figure  12.   Air  Inlet  Thermocouple  Probe  Arrangement 
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NOTE:   Arrangement  is  as  viewed  from  top  front 


Figure  13.   Air  Outlet  Thermocouple  Probe  Arrangement 
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Right  Hand  Heater  Mounting  Block  Thermocouple 
Arrangment 
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Figure  16.   Right  Hand  Heater  Sidewall  Thermocouple  Arrangement; 
Rear  View 
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Figure  17.   Left  Hand  Heater  Sidewall  Thermocouple  Arrangement; 
Rear  View 
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Three  insulated  thermocouples  were  mounted  to  the  outside 
surface  of  each  of  the  plexiglas  heater  mounting  blocks,  as 
shown  in  Figures  14  and  15,  in  order  to  determine  the  mounting 
block  outer  surface  temperatures.   These  temperatures  were 
used  to  estimate  the  amount  of  heat  lost  through  the  back  of 
the  heater  assemblies  as  explained  in  section  HIE.   Thermo- 
couples were  also  mounted  on  the  inside  and  outside  surfaces 
of  the  plexiglas  sidewalls,  as  shown  in  Figures  16  and  17, 
in  order  to  estimate  the  losses  to  the  atmosphere  and  to  the 
bed  itself  from  the  thin  plexiglas  wall  strips  adjacent  to 
the  edges  of  the  heater  assemblies.   A  final  set  of  thermo- 
couples was  mounted  to  the  inside  and  outside  surfaces  of 
the  container  rear  wall,  as  shown  in  Figure  18,  in  order  to 
obtain  an  estimate  of  losses  through  the  rear  wall.   Because 
of  symmetry  it  was  assumed  that  the  losses  through  the  front 
wall  were  identical  to  the  rear  wall  losses. 

All  thermocouples  used  in  the  test  apparatus  were  sampled 
by  an  HP  3497A  data  acquisition  system  controlled  by  an  HP  85 
desk  top  computer.   This  system  was  capable  of  sampling  and 
storing  the  output  of  seventy  five  thermocouples  within  seven 
seconds.   Because  of  the  steady  state  conditions  in  existance 
during  each  of  data  runs,  the  seven  second  sampling  time  was 
considered  brief  enough  to  treat  the  thermocouple  readings 
as  through  they  had  all  been  taken  simultaneously. 

Once  all  of  the  thermocouples  had  been  sampled,  the  HP  85 
retrieved  the  stored  thermocouple  voltages  and  converted  them 
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to  temperature  values  using  equations  previously  developed. 
These  temperature  values,  along  with  other  stored  parameters,. 
were  then  recorded  on  magnetic  disks  using  a  Hewlett  Packard 
Model  HP  82901M  dual  flexible  disk  drive.   This  recorded 
data  was  later  processed  and  permanently  displayed  using  an 
HP  82905B  printer  and  an  HP  7225A  plotter. 

Computer  programs  used  to  control  data  acquisition  and 
data  processing  are  shown  in  Appendixes  J  and  K,  respectively. 
The  arrangement  of  the  heat  transfer  data  collection  and  pro- 
cessing system  is  shown  in  Figure  19  while  Figures  20  and  21 
show  photographic  views  of  the  entire  experimental  system. 

C.   EXPERIMENTAL  APPARATUS  CALIBRATION 

The  inlet  air  flowmeter  was  calibrated  using  an  ASME 
Herschel -  type  venturi  with  an  inlet  diameter  of  4.26  inches 
and  a  throat  diameter  of  2.13  inches.   Experimental  data 
(shown  in  Appendix  B)  was  collected  for  nine  different  data 
points  and  from  this  data,  volumetric  air  flow  rates  were 
calculated  using  equations  shown  in  Appendix  H.   A  plot  of 
air  flowrate  vs  flowmeter  reading  is  shown  in  Figure  22.   A 
least  squares  curve  fit  was  applied  to  the  data  points  to 
obtain  a  linear  equation  which  converts  flowmeter  readings 
to  air  flow  rates  in  CFM.   This  equation  was  then  incorporated 
into  the  data  reduction  computer  program  shown  in  Appendix  K. 
The  rotameter  calibration  data  shows  that  a  100*  flowmeter 
reading  equates  to  34. .0  CFM. 
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Figure  18.   Rear  Wall  Thermocouple  Arrangement 
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Heat    Transfer   Data   Collection   System   Diagram 
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Thermocouples  number  1  through  67  were  calibrated  using 
a  Rosemount  Commutating  Bridge  Model  920A  and  a  Rosemount 
Calibration  Bath  Model  913A  using  liquid  nitrogen  and  water. 
The  HP  5497A  and  the  HP  85  were  included  in  the  system  during 
calibration,  as  shown  in  Figure  23.   Data  was  collected  for 
twenty  different  bath  temperatures  and  a  least  squares  curve 
fit  was  applied  to  the  data  to  obtain  a  fourth  order  poly- 
nomial function  for  each  of  the  thermocouples.   The  constants 
from  these  functions  were  then  incorporated  into  the  data 
acquisition  computer  program  shown  in  Appendix  J  to  allow 
automatic  conversion  of  thermocouple  output  voltage  to  tem- 
perature, in  degrees  Fahrenheit.   A  sample  of  the  calibration 
data  in  shown  in  Appendix  D. 

Thermocouples  number  68  through  79  were  not  calibrated 
but  were  instead  assigned  constants  from  the  standard  power 
expansion  for  Type  T  thermocouples.  [Ref.  5]   They  were  then 
compared  in  a  water  bath  of  known  temperature  to  confirm 
proper  operation. 

The  water  manometer  used  to  determine  static  pressure 
drop  within  the  bed  was  regularly  rechecked  for  zero  readings 
and  was  readjusted  as  necessary. 

D.   EXPERIMENTAL  PROCEDURE 

Preparation  for  the  taking  of  heat  transfer  data  was 
accomplished  by  adjusting  the  movable  side  wall  to  the  proper 
position  for  the  geometry  to  be  studied.   The  baffle  plate 


51 


was  then  clamped  into  position  and  the  joints  were  sealed 
with  silicone  rubber  caulking.   This  caulking  was  allowed 
to  harden  for  a  minimum  of  twenty  four  hours  to  ensure  a 
leak  tight  seal.   Carefully  weighed  quantities  of  glass 
beads  were  then  added  to  the  bed  through  the  right  hand 
air  exhaust  tube  until  the  depth  of  the  bed  reached  12.5 
inches  above  the  distributor  plate.   The  air  supply  turbo 
compressor  was  then  lined  up  and  energized.   Following  this, 
the  air  inlet  control  valve  was  opened  and  adjusted  until 
the  desired  air  flow  rate  registered  on  the  flowmeter.   The 
heater  power  supplies  were  then  energized  and  adjusted  to 
the  desired  level. 

The  power  supply  controls  were  set  so  that  the  voltage 
going  to  each  of  the  heaters  was  the  same.   Current  into 
the  heaters  was  determined  by  the  resistance  of  the  nickel- 
heating  element  which  was  in  turn  determined  by  the  temper- 
ature of  the  heaters.   Once  the  voltage  was  equilized,  the 
bed  was  allowed  to  come  to  steady  state  conditions.   Prior 
to'  the  taking  of  data,  the  system  was  checked  for  thermal 
equilibrium  by  monitoring  the  temperature  at  specific  points 
using  a  special  computer  program  which  sampled  thermocouples 
at  five-second  intervals.   Once  it  was  confirmed  that  the 
heater  assemblies  and  the  bed  itself  were  at  a  steady  tem- 
perature, the  data  acquisition  program  was  initiated.   At 
the  beginning  of  the  program  the  computer  requested  that 
the  operator  provide  certain  information  regarding  the  data 
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run  to  be  recorded.   The  first  item  of  information  requested 
was  the  run  number.   ThLs  number  was  structured  so  as  to  be 
as  informative  as  possible.   The  first  two  digits  represented 
the  bed  width  in  inches.   These  two  digits  were  followed  by 
a  letter  indicating  which  of  the  heaters  was  energized;  "L" 
designated  left,  "R"  stood  for  right,  and  "B"  designated  both 
This  letter  was  followed  by  two  digits  indicating  the  air 
flowmeter  reading.   Finally,  a  single  letter  was  used  for 
special  purpose  designations  such  as  a  repeated  run  without 
a  change  in  configuration  or  flow  rate.   The  next  item  of 
information  requested  by  the  program  was  the  date  and  time. 
This  was  followed  by  requests  for  bed  width  and  height,  air 
flow  rate,  ambient  temperature  (read  on  a  mercury  in  glass 
thermometer  positioned  near  the  test  apparatus),  voltage 
and  current  readings  from  the  two  power  supplies  to  each  of 
the  heaters,  the  number  of  bed  thermocouple  probes  in  use, 
and  the  amount  of  bed  expansion,  if  any.   Finally,  general 
comments  concerning  observed  bed  behavior  were  entered. 
As  the  push  of  a  button,  all  of  the  thermocouples  were 
then  sampled  and  the  corresponding  temperatures  were  then 
recorded  on  disc  along  with  the  other  data  entered  from  the 
keyboard.   Once  the  recording  process  was  completed,  the 
system  was  ready  to  begin  another  data  run.   If  the  next 
run  to  be  studied  was  for  the  same  bed  configuration,  all 
that  had  to-be  done  was  adjust  the  air  flow  rate  to  the  new 
reading  and  wait  for  the  temperatures  to  stabilize.   In 
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some  cases,  the  power  supplied  to  the  heaters  was  also  ad- 
justed in  order  to  maintain  as  large  a  temperature  differ- 
ential as  possible  between  the  incoming  air  and  the  heaters. 
In  order  to  avoid  excessive  softening  of  the  plexiglas 
adjacent  to  the  heaters,  the  maximum  heater  temperature  was 
limited  to  approximately  one  hundred  eighty  degrees 
Fahrenheit . 

After  data  was  collected  over  the  entire  range  of  flow 
rates  for  a  given  geometry,  the  caulking  was  removed  and  the 
movable  wall  was  repositioned  for  the  next  configuration. 
The  entire  data  collection  process  was  then  repeated.   In 
all,  four  different  geometric  configurations  were  studied. 
Bed  widths  of  6.25  inches,  8  inches,  10  inches,  and  12  inches 
were  used.   In  addition,  data  runs  were  taken  in  the  6.25 
inch  configuration  with  the  bed  empty.   Data  was  also  col- 
lected with  the  underside  of  the  distributor  to  the  left  of 
the  movable  wall  masked  off  as  a  check  to  see  if  the  baffle 
plate  was  performing  effectively.   A  total  of  eighty  data  runs 
were  completed  during  this  study.   Experimental  heat  transfer 
data  as  well  as  calculated  results  are  shown  in  Appendix  C. 

Pressure  drop  data  was  collected  for  each  bed  configura- 
tion by  aligning  the  valve  manifold  (shown  in  Figure  5)  so 
that  the  static  pressure  tap  immediately  above  the  distri- 
butor plate  was  connected  to  the  30  inch  water  manometer. 
The  fluidizing  air  flow  rate  was  then  adjusted  to  the  desired 
level  and  the  pressure  reading  at  the  base  of  the  bed  was  read 
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on  the  manometer  and  recorded.   The  manifold  was  then  re- 
aligned so  that  the  static  pressure  tap  near  the  top  of 
the  bed  housing  was  connected  to  the  manometer.   This  read- 
ing was  then  read  and  recorded.   By  substracting  this  read- 
ing from  that  obtained  at  the  base  of  the  bed,  the  bed 
pressure  drop  was  determined.   Pressure  drop  data  was  ob- 
tained in  this  fashion  over  the  entire  flow  range  for  each 
of  the  bed  configurations.   During  the  pressure  drop  data 
runs,  the  power  to  each  of  the  heaters  was  secured. 

E.   EXPERIMENTAL  DATA  PROCESSING 

All  raw  data  was  processed  using  an  HP  85  desk  top 
computer.   The  final  results  consisted  of  six  pages  of 
printed  data  for  each  run.   The  first  page  was  a  listing  of 
experimental  conditions  and  comments  recorded  during  the 
data  acquisition  phase  of  the  study. 

The  second  page  was  a  listing  of  all  thermocouple  out- 
puts and  their  corresponding  temperatures. 

The  third  page  showed  a  temperature  profile  of  the  right' 
hand  heater  surface  as  seen  from  within  the  bed.   In  addition, 
an  average  value  of  the  thermocouple  readings  in  each  hori- 
zontal row  was  calculated  and  displayed.   Thermocouple  number 
four  was  found  to  be  malfunctioning,  therefore  it's  output 
was  disregarded.   In  a  separate  display,  the  average  row 
temperatures  were  plotted  against  bed  depth  to  show  how  the 
heater  temperatures  varied  from  top  to  bottom. 
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The  fourth  page  of  the  results  printout  showed  the  left 
hand  heater  temperature  profile  in  a  format  identical  to 
that  of  page  three. 

The  fifth  page  showed  the  temperature  profile  of  the  bed 
itself,  as  measured  by  the  thermocouple  probes  immersed 
within  it.   The  profile  is  arranged  as  seen  from  above  the 
bed. 

The  sixth  and  final  page  contained  the  calculated  results. 
The  first  value  displayed  on  this  page  was  the  average  heater 
temperature  (Til .   This  value  was  a  simple  arithmetic  average 
of  the  five  horizontal  row  temperatures  previously  calculated. 


Th  =  {R(l)  +  R(2)  +  R(3)  +  R(4)  +  R(5)}/5  (3.1) 


The  average  bed  temperature  (T,)  was  similarly  calculated 
using  the  temperatures  displayed  on  page  five  of  the  printout 


Tb  =  (T(50)  +  T(51)  +  T(52)  +  T(53)  +  T(54)  +  T(56) 

+  +  T(X0}/(X-5O)  (3.2) 


X  =  Highest  number  bed  probe 

The  electrical  energy  into  the  heater  (qie)  was  calculated 
from  the  voltage  and  current  measurements  read  off.   The  heater 
power  supplies  as  shown  below: 
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qie  ■  I  x  E  x  0.0569  (  Wagji;Min)  x  60(^1)  (3.3) 


An  estimate  of  the  heat  less  from  the  heater  to  the  atmos- 
phere, through  the  heater  mounting  block  (q  ,  )  was  obtained 
as  follows:   First  an  average  temperature  (T  ,  )  was  calculated 
for  the  three  thermocouples  attached  to  the  outer  surface  of 
the  heater  mounting  block. 


Tmb  =  Cn41)  +  T^42)  +  T(43)}/3  (3.4) 


Next  it  was  assumed  that  the  thermal  contact  resistance 
between  the  copper  plate  and  the  strip  heater  was  negligible, 
therefore  the  temperature  of  the  strip  heater  surface  was 
equal  to  TV .   The  contact  resistance  between  the  strip  heater 
and  the  insulation  was  also  assumed  to  be  negligible. 

The  heater  insulation  and  mounting  block  was  then  con- 
sidered to  be  a  composite  wall  as  shown  in  Figure  24.   Fourier's 
Law  of  Conduction  was  then  utilized  to  calculate  the 'heat  flux 
through  the  mounting  block. 

Th  "  Tmb  (3.5) 


lmb      L.        L  , 

(_i_)  +  (-JSL_) 


k-A.      k  .  A  , 
11      mb  mb 
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L ■  =  L  ,  =  0.5  inches 
1    mb 

? 

A.  =  A  ,  =  5  in.  x  10  in.  =  50  in 
i     mb 

ki  =  °'°25  HrBFt-0F 

kmb=  °-112  Hr-Ft-UF 


Next  a  convection  heat  transfer  coefficient  fh  .  )  from 

v  mba 

the  mounting  block  to  the  atmosphere  was  calculated  using 
Newton's  Law  of  Cooling. 


hmba  = C3.63 

A  .  (T  .  -  T  ) 
mb v  mb     a' 


Then  an  estimate  was  made  of  the  heat  lost  to  the  atmos- 
phere from  the  hot  plexiglas  side  wall  sections  adjacent  to 

the  side  edges  of  the  heater  (q   ) .   This  estimate  was  based 

s  a 

on  the  side  wall  outer  surface  temperature  and  used  Newtons 
Law  of  Cooling  and  the  convection  heat  transfer  coefficient 
previously  calculated  for  the  heater  mounting  block.   Since 
the  side  wall  sections  and  the  mounting  block  were  in  such 
close  proximity  it  was  assumed  that  the  coefficients  would 
be  in  close  agreement  with  one  another. 

q    =  h  .   x  A   x  (T(72)  -  T  }  (3.7) 

Hsa    mba    s  a 

1  Ft 
A   =  0.5  in.  x  10  in.  x  Tini — . — T  x  2 
s  144  m^ 

(one  strip  on  each  side  of  heater) 
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A  similar  loss  estimate  was  then  calculated  for  the  heat 
flux  to  the  atmosphere  from  the  side  wall  sections  .adjacent  to 
the  top  edge  of  the  heater  (q   ) . 

*ta  =  hmba  x  At  x  {T^68^  "  Ta}  (3.8) 

At  =  0.5  in.  x  6  in.  x  i-E£ — - 

144  in 

Yet  another  estimate  was  then  calculated  for  the  loss  to 
the  atmosphere  from  the  side  wall  section  adjacent  to  the 
bottom  edge  of  the  heater  (a,  ). 


<»ba  "  hmba  x  Ab  x  {TC71'  '  V  (3.9) 


A,  •  =  0.5  in.  x  6  in. 


1  FtZ 


144  in- 


All  atmospheric  losses  were  then  totaled. 

^a  =  ^mb  +  ^sa  +  ^ta  +  %a  C3-10) 

The  difference  between  the  electrical  power  input,  q  -, 
and  losses  from  the  heater  to  the  atmosphere,  q  ,  is  the  heat 
flux  transferred  to  the  bed.   In  order  to  identify  that  portion 
transferred  directly  from  the  heated  copper  plate  to  the  bed, 
an  estimate  of  the  heat  flux  from  the  heated  plastic  strip 
surrounding  the  copper  plate  to  the  bed  must  be  made.   This 
is  done  by  assuming  that  same  heat  transfer  coefficient  applies 
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to  both  the  copper  plate  and  to  the  plastic  strip.   The  method 
selected  was  an  iterative;  wherein  it  was. first  assumed  that 
the  heat  flux  passed  through  the  copper  plate,  from  which  the 
heat  transfer  coefficient  was  calculated.   Using  this  coeffi- 
cient, the  heat  flux  from  the  heated  plastic  strip  was  cal- 
culated which  in  turn  yielded  a  new  heat  flux  from  the  copper 
plate.   The  process  was  repeated  until  the  heat  flux  values 
did  not  change  significantly.   Thus  following  this  procedure 

the  heat  flux  into  the  bed  from  the  heater  q  .  was  assumed 

ncb 

to  be  equal  to  the  incoming  electrical  power  q  ,  minus  the 
losses  from  the  heater  to  the  atmosphere  q  . 


[cb 


=  q 


ei 


(3.11) 


Once  the  heat  flux  through  the  copper  plate  q  ,  was  known 
the  plate  outside  surface  temperature  T   seen  by  the  bed 
could  be  calculated. 


°*cb  x  L< 


Tp   Th  "  (jT 


(3.12) 


L   =  0.25  in.  x  (  *  Ft) 
c  k12  in' 


k      =    232    -     (0.032    x    (Th    -    70)}     Btu/Hr-Ft-°F 


1    Ft 


A„    =    5    in.    x    10    in.    x    (ttti — ■ — 2) 
c  v144    mL 
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Once  the  outside  plate  surface  temperature  was  known,  the 
convection  heat  transfer  coefficient  from  the  plate  to  the 
bed  (h  ,  )  could  be  calculated. 


cb 


Lcb 


A„(T   -  T\J 
c  v  p     bJ 


(3.13) 


Using  the  plate  to  bed  heat  transfer  coefficient  (h  ,  ) 
just  determined,  an  estimate  was  made  of  the  heat  flux  which 
reached  the  bed  from  the  hot  plexiglas  strips  adjacent  to  the 
heater  sides  (q  ,  )  . 


q  K  =  h  K  x  A   x  {T(7  3)  -  T,  } 
Msb     cb     s  b 


(3.14) 


Similar  estimates  were  then  prepared  for  the  flux  entering 
the  bed  from  the  hot  sidewall  section  adjacent  to  the  top  edge 
of  the  heater  (qtu)  as  well  as  the  section  adjacent  to  the 
bottom  edge  of  the  heater  (q,  ,  ) . 


«tb  =  hcb  x  At  x  {T(69)  "  V 


(3.15) 


%b    -   hcb  x  Ab  x  {T(70)  "  V 


(3.16) 


The  heat  flux  into  the  bed  from  all  sources  other  than  the 
heater  plate  itself  were  then  added  to  obtain  a  total  value 

^ob}- 
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%b   =  qsb  +  qtb  +  %b  (3'17) 

The  original  copper  plate  heat  flux  value  was  then  cor- 
rected to  allow  for  that  amount  (q  ,  )  which  did  not  pass 
through  the  plates. 


q-.~q.-q      -a*  (3.18) 

Hcb    Mei    Ma    Mob  ^     * 

It  was  then  necessary  to  recalculate  all  of  the  parameters 
dependent  on  this  value.   This  process  was  repeated  in  an  iter- 
ative fashion  until  the  calculated  flux  through  the  plate  did 
not  change  signf iciant ly .   Once  the  iterative  process  had  been 
completed,  the  final  values  were  printed  out  and  the  entire 
calculation  process  was  repeated  for  the  left  hand  heater. 
The  left  hand  heater  calculations  were  similar  to  those  used 
for  the  right  hand  with  the  exception  of  variable  names  and 
thermocouple  numbers.   Because  of  a  limitation  on  the  total 
number -of  thermocouples  which  could  be  controlled  by  the 
HP  3497A,  there  were  no  thermocouples  available  to  measure 
the  inside  and  outside  surface  temperature  of  the  plexiglas 
sidewall  section  adjacent  to  the  bottom  edge  of  the  left 
hand  heater.   This  prevented  the  calculation  of  a  loss 
estimate  based  on  a  direct  measurement.   To  compensate  for 
this,  it  was  assumed  that  the  losses  from  the  sidewall 
sections  adjacent  to  the  top  and  bottom  edges  of  the  heater 
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were  identical  to  one  another.   Once  the  left  hand  heater 
computations  were  completed  and  the  results  displayed,  cal- 
culations were  begun  on  an  energy  balance.   The  objective 
was  to  determine  the  total  amount  of  energy  entering  the  bed 
and  then  compare  that  with  the  amount  calculated  to  be  leavins 
in  order  to  validate  the  data  and  the  calculations  used.   The 
first  parameter  to  be  determined  was  the  density  of  the  air 

within  the  bed  p  . 

a 


1.32  5  x  P 

pa  = 


(T  +459.69)  (3.19) 

a 

P   =  29.92  in.  Hg 

a  ° 

T   =  T  (55) 
a     ^ 


Since  the  air  entering  a  fluidized  bed  is  known  to  assume 
the  bed  temperature  within  a  very  short  distance  of  entry, 
the  air  outlet  temperature  was  used  as  the  bed  air  temperature 

Using  this  value  of  density,  the  amount  of  energy  carried 
out  of  t-he  bed  by  the  fluidizing  air  stream  (qao)  was 
calculated. 


^ao    =   °   Cp    (Tao    "    Tai} 


(3.20) 


q         =   V  x    p   ,     x    60    Min/Hr    x    0.241    x    (T         -    T    .)  (3.21) 

Hao  ab  a<J  <il 

T    .  =    {T(47)    +    T(48)    +    T(49)}/3 

3. 1 

T  =    T(55) 
ao 
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The  heat  loss  through  the  front  and  rear  walls  (q,   )  was 

to  nfra' 

then  estimated.   Only  the  rear  wall  was  instrumented,  however,. 
because  of  symmetry,  it  was  assumed  that  the  loss  through 
the  front  wall  would  be  identical. 


kfr  x  Afr  x  {T(75)  -  T(7d)} 
qfra  =  -S d^ (5.22) 

Lfr 

kr        =  0.112  Btu/Hr-Ft-°F 
f  r 

1  Ft 
Ar    =  12  in.  x  W  x  (t^ — ■ —  )  x  2  (to  account  for 
f r  v12  in. 

both  front  and  rear  walls) 

W     =  Bed  Width 

Lr        =  0. 5  in. 

f  r 


The  total  heat  flux. leaving  the  bed  q   was  then  calculated. 
It  was  assumed  that  the  only  flux  paths  out  of  the  bed  were 
via  the  front  and  rear  walls  and  via  the  outgoing  air  stream. 
The  side  walls  were  considered  to  be  at  a  sufficiently  elevated 
temperature  as  to  block  the  flow  of  heat  outward  from  the  bed. 
The  floor  of  the  bed  was  considered  to  be  free  of  losses -since 
any  heat  passing  through  it  would  be  returned  to  the  bed  by 
the  fluidizing  air  stream  as  it  passed  through  the  distributor. 

q   =  q    +  qr  (3.23) 

Mo    Moa    Mfra 

The  total  energy  into  the  bed  from  all  right  hand  sources 
was  then  calculated  (q,  ) .   This  included  both  the  flux 
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through  the  plate  and  the  flux  through  the  sidewalls 
adjacent  to  the  heater. 

%r    =  %i    "  ^a  ^'24) 

A  similar  calculation  was  performed  for  the  left  hand 
sources  (q,  J  and  then  the  total  energy  into  the  bed  from 
all  sources  was  calculated  (q^) • 

%     =  %T     +  %1  (5'25) 

The  superficial  gas  velocity  (U)  (in  Ft/Sec)  was  calculated 
using  the  following  equation: 

V  x  (4-^) 
U  =   5U  Sec  (3.26) 

Ad 

1  Ft 
A,  =  Distributor  Area  =  6  in.  x  W  x  (^44  in2) 

The  superficial  gas  mass  velocity  (G)  was  then  calculated. 


G  =   Lbm9   =  2L.  =  u  (-^-)  x  o      0*2)    x  (-5600  Sec) 
Hr-FtZ     Ad       Sec      a   FtJ        Hr 


G  =  U  x  p   x  3600  (3.27) 

cL 
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The  final  calculation  of  the  program  was  the  determination 
of  the  particle  Reynolds  Number  (Re) . 


p  V  D 
Re  =  £ 


G  x  D 


3600  x  u  x  g 


(3.28) 


u  =  3.96  x  10"7  Lbf-Sec/Ft2 

g  =  32.174  Lbm-Ft/Lbf-Sec2 

D  =  0.012  in.  x  (A    *Z) 
p  12  in. 


Following  the  completion  of  this  calculation  the  energy  balance 
section  of  page  six  was  printed  out  and  the  program  was 
terminated . 
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IV.   PRESENTATION  AND  DISCUSSION  OF  RESULTS 

A.   FLOW  PATTERNS 

The  particle  flow  patterns  observed  in  the  course  of  this 
study  were  in  general  agreement  with  those  described  in  the 
literature  for  a  bubbling  bed.   The  particles  moved  upward 
near  the  centerline  of  the  bed  as  did  most  of  the  bubbles. 
Once  the  particles  reached  the  surface  of  the  bed  they  were 
thrown  outward  from  the  center  and  then  proceeded  to  travel 
downward  along  the  bed  walls  as  shown  in  Figure  25.   The 
area  of  greatest  activity  was  roughly  cylindrical  in  shape 
with  the  corner  areas  exhibiting  considerably  less  particle 
motion.   Observed  particle  stagnation  regions  are  shown  in 
Figure  26.   No  indications  of  channeling  were  observed  visually 
at  any  time.   At  high  fluidization  flow  rates  in  the  narrow 
bed  configuration,  slugging  was  occasionally  observed.   This 
was  characterized  by  the  appearance  of  large  bubbles  and  the 
sudden  lifting  of  the  bed  surface  to  a  level  two  to  three 
inches  above  its  previous  position.   Once  the  bubble  reached 
the  top  of  the  bed,  the  surface  would  suddenly  collapse  and 
resume  its  former  level  until  the  next  bubble  appeared. 

The  particles  within  the  bed  were  observed  to  behave  in 
a  manner  much  like  a  viscous  fluid.  The  layer  of  particles 
directly  adjacent  to  the  walls  showed  little,  if  any,  move- 
ment.  Because  of  their  transparent  composition  it  was  possible 
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Figure  25.   Observed  Particle  Circulation  Flow  Paths 


:•  V.' 
:■'.  '■ 


Figure  26 


Observed  Particle  Stagnation  Regions 


69 


to  observe  that  the  next  layer  of  particles  was  moving  at  a 
slightly  higher  velocity  and  that  the  third  layer  inboard 
was  moving  even  faster.   This  movement  was  not  steady  but 
tended  to  occur  in  surges.   It  is  possible  that  a  static 
charge  contributed  to  the  "No  slip"  condition  observed  at 
the  boundary.   After  the  fluidizing  air  flow  was  secured  upon 
completion  of  a  test  run  it  was  noted  that  numerous  particles 
were  clinging  to  the  container  walls  in  the  upper  portion  of 
the  bed.   In  an  attempt  to  eliminate  this  condition,  an  anti- 
static film  was  sprayed  on  the  upper  portions  of  the  bed  walls 
For  a  short  period  of  time,  the  walls  remained  free  of  cling- 
ing particles,  however,  the  effect  did  not  last  for  more  than 
two  or  three  data  runs.   Following  this,  no  further  attempts 
were  made  to  alter  this  phenomenon.   Research  by  Miller  and 
Logwinuk  [Ref.  2]  suggested  that  the  presence  of  such  a  static 
charge  could  contribute  to  erratic  heat  transfer  results,  how- 
ever for  the  purpose  of  this  study,  any  affect  was  considered 
to  be  negligible .  . 

Closer  examination  of  the  bubble  flow  path  showed  that 
many  bubbles  tended  to  originate  from  the  left  side  of  the 
bed  near  the  bottom  of  the  movable  side  wall.   As  a  result, 
there  was  a  shift  in  overall  particle  activity  to  the  left 
of  center  and  a  corresponding  alteration  of  the  quiescent 
zone  near  the  bottom  of  the  bed  such  that  the  majority  of 
stagnant  particles  were  located  in  the  right  hand  corner. 
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Figure  27.   Photograph  of  Off-Center  Bubble 
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This  phenomenon  was  particularly  noticeable  in  the  wide  bed 
configurations.   An  example  of  this  off-center  activity  is 
shown  in  Figure  27. 

Examination  of  the  apparatus  revealed  that  a  stream  of 
air  was  being  injected  into  the  bed  from  the  area  between 
the  baffle  plate  and  the  distributor.   Air  from  the  inlet 
plenum  was  passing  through  the  unused  portion  of  the  dis- 
tributor to  the  left  of  the  movable  wall.   Once  the  air  was 
through  the  distributor,  further  upward  progress  was  blocked 
by  the  baffle  plate  and  the  sealed  joints.   There  was  however, 
a  very  thin  horizontal  passage  between  the  baffle  and  the  top 
of  the  distributor  as  shown  in  Figure  28.   The  air  flowed 
along  this  path  until  it  emptied  into  the  bed  itself.   Because 
of  this  "peripheral  injection"  there  was  an  uneven  distribution 
of  fluidizing  air  and  a  resulting  shift  in  the  center  of  activ- 
ity toward  the  left  hand  side.   In  order  to  determine  the  effect 
of  this  phenomenon  an  entire  set  of  data  was  taken  with  the 
unused  portion  of  the  distributor  masked-off  from  the  under- 
side so  that  no  air  could  enter.   During  these  ma'sked  data 
runs  it  was  observed  that  the  bubble  path  was  very  well 
centered  as  was  the  quiescent  zone  near  the  bottom  of  the  bed. 
All  other  data  runs  were  taken  with  the  distributor  unmasked 
in  order  to  study  the  difference  in  heat  transfer  performance 
between  the  two  sides.   A  detailed  study  of  particle  cell 
flow  patterns  for  this  apparatus  was  presented  by  Morgan 
[Ref.  2]  in  1981.  • 
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B.   PRESSURE  DROP  DATA 

Measurements  were  taken  of  the  pressure  drop  across  the 
bed  for  all  flow  rates  and  bed  configurations.   This  data 
can  be  found  in  Appendix  A.   Plots  were  constructed  of  LOG 
AP  vs.  LOG  G  for  each  bed  configuration  in  order  to  study 
the  variation  in  pressure  drop  with  changes  in  the  super- 
ficial mass  velocity.   These  plots  are  shown  in  Figures  29 
through  32.   An  ideal  pressure  drop  curve  is  shown  in  Figure 
33.   Each  of  the  plots  display  the  same  basic  characteristics 
of  the  ideal  curve  in  that  the  pressure  drop  increases 

steadily  until  the  point  of  minimum  fluidization  (G  r)  is 

1  r  mf 

reached.   From  that  point  onward,  the  pressure  tends  to 
remain  relatively  constant  despite  further  increases  in  air 
flow  rate.   This  is  because  of  the  fact  that  once  the  bed  is 
fluidized,  any  further  increase  in  the  flow  rate  imparts  an 
increased  lift  to  the  particles  because  of  the  viscous  drag 
force  exerted  on  them  by  the  moving  air  stream.   As  the  lift 
increases,  expansion  of  the  bed  occurs. and  voidage  (c)  in- 
creases.  As  the  voidage  increases,  the  interstitial  spacing 
between  the  particles  also  increases.   This  results  in  a 
reduction  of  the  interstitial  gas  velocity,  which  in  turn, 
lowers  the  drag  exerted  on  the  particles  until  the  gravita- 
tional force  and  the  lift  are  once  again  in  balance.   As  a 
result  of  the  increased  voidage,  the  pressure  drop  across 
the  bed  remains  constant  or  may  even  decrease. 
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Examination  of  the  pressure  drop  curve  for  the  6.25  inch 
configuration  (Fig.  29)  shows  a  deviation  from  ideal  behavior 
in  that  the  pressure  drops  sharply  after  reaching  the  point 
of  minimum  fluidization.   A  pressure  drop  of  this  sort  is 
normally  associated  with  channeling.   As  indicated  previously, 
no  visual  indications  of  channeling  were  detected,  however 
it  is  possible  that  "intermittent  channeling"  was  taking 
place.   This  type  of  channeling  would  be  difficult  to  differ- 
entiate from  rapid  bubbling  because  the  channel  does  not 
extend  through  the  entire  depth  of  the  bed  as  it  does  in 
"through  channeling."   An  alternate  explanation  for  the  drop 
in  pressure  is  the  effect  of  particle  interlocking  which  some- 
times manifests  itself  during  the  initial  fluidization  of  a 
bed.   It  is  noted  that  subsequent  data  runs  did  not  demonstrate 
such  behavior. 

Based  on  the  curve  shown  in  Figure  29  the  point  of  minimum 

fluidization  for  the  6.25  inch  configuration  occurred  at  G=159.6 

2 
Lbm/Hr-Ft   (26%  flow).   It  was- was  noted  the  visual  observation 

of  fluidization  onset  occurred  at  the  same  air  flow  rate  as 

was  predicted  by  this  curve. 

The  pressure  drop  curve  for  the  8  inch  configuration 

(Fig.  30)  is  in  very  close  agreement  with  the  ideal  curve 

with  no  indications  of  channeling  or  slugging.   Based  on  this 


curve  the  point  of  minimum  fluidization  occurred  at  G=132.6 


7 
Lbm/Hr-Ft   (281  flow).   Fluidization  was  visually  observed 


2 
to  commence  at  137  Lbm/Hr-Ft   (29%  flow). 
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Figure    29. 


Pressure  Drop  vs  Superficial  Mass  Velocity; 
6.25  in.  Bed 
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Figure  30.   Pressure  Drop  vs  Superficial  Mass  Velocity; 
8.0  in.  Bed 
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Figure  32 


Pressure  Drop  vs  Superficial  Mass  Velocity; 
12.0  in.  Bed 
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Figure  33.   Pressure  Drop  vs  Superficial  Mass  Velocity; 
Ideal 
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The  curve  for  the  10  inch  configuration  (Fig.  31)  was 
smoother  than  the  previous  configurations  and  matched  the 
ideal  curve  closely  except  for  a  gradual  decrease  in  pres- 
sure following  fluidization  onset.   The  reason  for  this 
decrease  was  not  readily  apparent.   The  curve  shows  that 
fluidization  commenced  at  G=148  Lbm/Hr-Ft2  (40%  flow). 
The  observed  starting  point  of  fluidization  was  also  at 
40%  flow. 

The  12  inch  configuration  curve  (Fig.  32)  was  also  very 
smooth  and  similar  to  the  ideal.   Like  the  10  inch  curve, 
a  slight  decrease  in  differential  pressure  was  noted.   The 
starting  point  of  fluidization  (G  r)    indicated  by  the  curve 
was  G=146.5  Lbm/Hr-Ft"  (48%  flow).   Fluidization  was  visually 
observed  to  commence  at  the  same  point. 

A  combined  graph  showing  the  pressure  drop  curves  for  all 
four  configurations  is  shown  in  Figure  34.   With  the  exception 
of  the  8  inch  configuration,  they  show  a  slight  gradual  de- 
crease in  the  minimum  fluidization  superficial  mass  velocity 

(G  r)  as  the  bed  width  is  increased.   A  similar  trend  was 
mf 

observed  by  Morgan  [Ref.  2]  who  attributed  this  to  a  reduction 
in  the  effect  of  the  wall  frictional  drag  as  a  fraction  of 
the  overall  drag  experienced  by  the  bed. 

The  average  bed  pressure  drop  at  minimum  fluidization 
(P  r)  for  all  configurations  was  approximately  16.0  in.  H90 
or  0.577  psi.  •  Research  by  Trivedi  et  al  [Ref.  6]  indicated 
that  the  pressure  drop  across  a  fluidized  bed  should  equal 
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the  weight  of  the  bed  material  per  unit  cross  sectional  area. 
The  settled  bed  density  for  the  glass  beads  used  in  this  ex- 
periment was  96  Lbm/Ft   which  resulted  in  a  weight-to-area 

2 
ratio  of  0.681  Lbm/In  .   Thus,  the  pressure  drop  observed 

was  slightly  less  than  the  theoretical  value.   Recalculation 

of  the  weight-to-area  ratio  based  on  the  expanded  bed  density 

yielded  essentially  the  same  results. 

C.   TEMPERATURE  DISTRIBUTION 

The  temperature  distribution  across  each  of  the  heater 
surfaces  was  examined  for  all  data  runs  and  two  distinct 
trends  were  found  to  be  evident.   First;  the  four  thermo- 
couple readings  in  any  horizontal  row  were  found  to  be  in 
very  close  agreement  with  one  another  in  all  cases  as  shown 
in  Appendix  E.   This  even  distribution  of  temperature  hori- 
zontally across  the  plate  suggested  that  the  local  rate  of 
heat  transfer  from  the  plate  was  the  same  for  any  given  point 
along  the  plate  surface  at  a  given  bed  depth.   This  conclusion 
was  based  on  the  assumption  that  the  heat  flux  into  the  plate 
was  also  horizontally  uniform.   The  second  trend  observed  was 
that  there  was  a  variation  in  heater  temperature  from  the  top 
of  the  plate  to  the  bottom  and  that  this  "temperature  profile" 
displayed  a  distinct  and  predictable  change  in  shape  as  the 
air  fluidization  flow  rate  was  increased.   Figures  35  through 
44  illustrate  these  trends.   At  low  air  flow  rates  the  point 
of  greatest  temperature  was  always  at  the  top  of  the  plate. 
Moving  downward,  the  temperature  then  gradually  dropped  in 
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a  non-linear  but  smooth  fashion  to  a  level  approximately  three 
to  five  degrees  less  than  that  observed  at  the  top  of  the 
plate.   Examination  of  the  plots  of  temperature  vs  bed  depth 
for  the  low  flow  rate  runs  showed  all  of  the  curves  to  have 
the  same  shape  with  each  curve  displaying  a  positive  slope. 
It  was  noted  that  as  the  flow  rates  increased,  the  slope  of 
these  curves  became  more  positive  until  they  assumed  a  ver- 
tical orientation  and  then  reversed  their  slope  to  become 
negative.   It  was  further  noted  that  this  slope  reversal 
occurred  at  or  near  the  point  of  fluidization  onset.   The 
change  in  slope  of  these  curves  meant  that  the  peak  tempera- 
ture within  the  plate  was  moving  downward  from  the  top  of 
the  plate  toward  the  bottom  as  the  fluidization  air  flow  rate 
was  increased.   A  possible  explanation  of  such  behavior  lies 
in  the  fact  that  prior  to  fluidization  onset,  when  the  par- 
ticles were  not  in  motion,  the  heater  was  primarily  cooled 
by  an  upward  flow  of  air  from  the  distributor.   Because  of 
this  upward  flow  path,  the  temperature  of  the  heater  tended 
to  increase  with  an  upward  movement  along  the  plate.   Once 
the  bed  was  fluidized,  however,  a  significant  amount  of 
cooling  was  then  provided  by  the  moving  particles  sweeping 
along  the  plate  surface.   Since  the  particles  were  moving 
downward  along  the  wall  a  thermal  gradient  was  introduced 
which  was  opposite  to  that  caused  by  the  upward  flow  of  air. 
The  peak  temperature  within  the  plate  would  consequently  be 
observed  at  the  point  where  these  two  gradients  met.   If  the 
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cooling  due  to  the  particles  was  greater  than  that  provided 
by  the  upward  air  flow,  the  peak  temperature  would  be  found 
toward  the  bottom  of  the  plate.   Conversely,  if  the  air  flow 
cooling  predominated,  the  peak  temperature  would  be  seen  in 
the  upper  section  of  the  plate.   An  increase  in  the  fluidiza- 
tion  air  flow  rate  would  provide  more  upward  flowing  cooling 
air  while  simultaneously  increasing  the  downward  moving 
particle  activity. 

Experimental  results  showed  that  an  increase  in  air  flow 
rate  caused  the  peak  temperature  to  move  downward.   This  sug- 
gested that,  as  the  flow  rate  was  increased,  the  increase  in 
particle  cooling  predominated  over  the  increase  in  air  cooling 
The  reason  for  this  predominance  may  have  been  that  the  par- 
ticles had  a  much  higher  specific  heat  capacity  than  the  air 
and  a  lower  mean  velocity.   Consequently,  they  were  able  to 
remove  a  larger  quantity  of  heat.   Another  possible  cause  for 
the  particle  cooling  predominance  as  flow  rates  were  increased 
was  that  the  motion  of  the  cooling  air  may  have  been  altered 
once  fluidization  commenced  so  that  more  of  the  upward  moving 
gas  was  concentrated  in  the  center  of  the  bed.   Although  there 
was  no  objective  evidence  that  this  was  occurring,  it  may  have 
been  that  any  increase  in  flow  rate  was  merely  diverted  to 
the  center  of  the  bed  and  consequently  did  not  increase  the 
amount  of  cooling  air  available  near  the  wall. 

Another  trend  which  was  observed  while  studying  the  heater 
temperature  distributions  was  that  as  the  bed  width  increased, 


the  shift  in  the  orientation  of  the  temperature  curve  became 
less  pronounced.    In  the  6  inch  configuration  the  peak 
temperature  point  moved  from  the  top  of  the  bed  downward  to 
the  very  bottom  as  fluidization  activity  increased.   The 
maximum  temperature  differential  observed  between  the  top 
and  the  bottom  of  the  plate  was  approximately  17  degrees. 
In  the  8  inch  configuration,  the  peak  temperature  point 
also  transited  down  the  entire  length  of  the  plate,  however 
it  did  so  at  a  slower  rate  than  in  the  6  inch  configuration. 
The  maximum  top- to -bottom  differential  was  essentially  the 
same  as  that  seen  for  the  6  inch  bed.   With  the  bed  con- 
figured to  the  10  inch  width,  the  peak  temperature  did  not 
move  down  the  entire  length  of  the  plate,  but  instead  stopped 
at  a  level  approximately  3  inches  above  the  bottom  edge.   The 
shift  in  temperature  orientation  was  even  more  gradual  than 
in  the  8  inch  configuration.   The  maximum  observed  tempera- 
ture differential  was  approximately  11  degrees.   While  in 
the  12  inch  configuration,  very  little  change  was  observed 
in  the  temperature  distribution  as  the  flow  rate  was  in- 
creased.  The  peak  temperature  shifted  from  the  top  of  the 
plate  to  just  above  center,  however  the  maximum  temperature 
differential  of  6  degrees  was  significantly  less  than  that 
observed  in  the  other  configurations. 

In  general,  the  right  hand  heater  temperature  distribu- 
tion proved  to  be  more  sensitive  to  changes  in  flow  rate 
than  did  the  left  hand  heater. 
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A  characteristic  unique  to  the  left  hand  heater  was  that 
immediately  prior  to  commencement  of  the  downward  transit  of 
the  peak  temperature  point,  the  temperature  differential 
between  the  top  and  the  bottom  of  the  plate  suddenly  increased 
to  approximately  three  times  its  previous  value.   As  the  flow- 
rate  was  further  increased,  this  differential  abruptly  dropped 
back  down  to  its  former  level.   This  trend  is  shown  in  Figures 
40  through  43.   It  is  believed  that  this  fluctuation  could 
be  related  to  a  possible  reorientation  of  air  flow  patterns 
as  fluidization  commences.   This  phenomenon  was  not  evident 
in  the  right  hand  heater  temperatures. 

Examination  of  the  data  in  Appendix  C  shows  that  for  each 
flow  rate  and  bed  configuration,  the  average  temperature  of 
the  left  hand  heater  was  lower  than  for  the  right  hand  heater. 

In  many  cases,  the  difference  in  cooling  of  the  plates 
was  quite  significant.   During  the  data  runs  with  the  distri- 
butor masked-off,  it  was  observed  that  the  left  and  right 
hand  heaters  were  being  cooled  to  approximately  the  same 
degree,  since  the  average  heater  temperatures  were  in  close 
agreement  with  one  another. 

The  temperature  profile  of  the  bed  itself,  as  measured 
by  the  immersed  probes,  showed  a  non-uniform  distribution 
prior  to  the  commencement  of  fluidization  for  all  configura- 
tions.  The  temperatures  were  highest  near  the  heater  sur- 
faces where  particle  heating  was  taking  place  by  conduction. 
Because  of  the  relatively  low  rate  of  heat  transfer  the 
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NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  35.   Right  Hand  Heater  Temperature  Profiles;  6.25  in 
Configuration 
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NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  36. 


Right  Hand  Heater  Temperature  Profiles;  6.25  in 
Masked  Configuration 
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NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  37.   Right  Rand  Heater  Temperature  Profiles; 
Configuration 


.0  in 


90 


Temp, 


NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  38.  .  Right  Hand  Heater  Temperature  Profiles;  10.0  in 
Configuration 
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NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  39.   Right  Hand  Heater  Temperature  Profiles;  12.0  in 
Configuration 
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NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  40. 


Left  Hand  Heater  Temperature  Profiles;  6.25  in 
Configuration 
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Temp. 


NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  41. 


Left  Hand  Heater  Temperature  Profiles;  6.25  in 
Masked  Configuration 
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NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  42 


Left  Hand  Heater  Temperature  Profiles; 
Configuration 
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NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  43.   Left  Hand  Heater  Temperature  Profiles;  10.0  in 
Configuration 
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NOTE:   Numerical  values  represent  flowmeter  readings  as  a 
percentage  of  maximum  rotometer  capacity. 


Figure  44.   Left  Hand  Heater  Temperature  Profiles;  12.0  in 
Configuration 
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temperature  tended  to  decrease  toward  the  center.   Immediately 
upon  commencement  of  f luidization ,  however,  the  temperature 
distribution  within  the  bed  became  extremely  uniform  and 
remained  so  as  fluidization  flow  rates  increased.   This  was 
directly  attributable  to  the  thorough  mixing  induced  by  the 
fluidized  particle  motion.   In  most  cases,  the  temperature 
of  the  outgoing  airstream  was  within  two  or  three  of  the 
temperature  of  the  bed  itself. 

D.   HEAT  TRANSFER  COEFFICIENTS 

In  order  to  establish  a  performance  baseline,  data  was 
collected  for  a  bed  from  which  all  of  the  particles  had  been 

removed  and  an  average  heat  transfer  coefficient  of  3.8  Btu/ 

2  o 
Hr-Ft  -  F  was  calculated.   All  remaining  runs  were  accom- 
plished with  the  apparatus  filled  to  a  settled  bed  depth  of 
approximately  12  inches.   The  resulting  heat  transfer  data 
for  all  runs  is  tabulated  in  Appendix  C  and  plots  of  heat 
transfer  coefficient  vs.  superficial  mass  velocity  are 
shown  in  Figures  45  through  53.   Each  of  these  plots  demon- 
strates the  same  basic  trend  in  that  initially,  during  low 
gas  flow  rate  conditions,  the  heat  transfer  coefficient  was 
correspondingly  low  (approximately  4.5  Btu/Hr-Ft  -  F)  and 
remained  low  until  a  flow  rate  approaching  the  minimum 
fluidization  superficial  mass  velocity  (G  f)  was  attained. 
From  this  point  onward,  the  heat  transfer  coefficient  in- 
creased in  a  roughly  linear  fashion  as  the  gas-  flow  rate  was 
increased . 
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Figure  45.   Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  6.25  in.  Configuration 
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Velocity;  6.25  in.  Masked  Configuration 
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Figure  47.   Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  Right  Hand  Heater,  6.25  in. 
Configuration 

101 


or 

LU 

H- 

<: 

LU 

X 

o 

z 

cc 

<: 

o 

cr 

X 

►— 

a 

ID 

i— 

»— 

GD 

3 

u. 

i— - 

CD 

LU 

£T 

i— i 

—1 

t— 

or 

CO 

►— 

o 

z 

M 

CO 

o 

a 

"— < 

CO 

I-* 

a 

> 

H- 

a 

< 

LU 

a 

-C 

£T 

iC 

LU 

3 

CO 

x 

LD 

<: 

CO 

»— « 

gg 

<: 

Lt_ 

2 

3l 

z 

3 

o 

B 

u 

D 

+ 

3 

* 

CJ 

z 

►— 1 

1 

tr> 

cm 

CD 


*+ 


I     i     1     i     I 


I     l     I 


I     i     I 


'■■•'■'•'•' 


-t  012 
-  S02 


-  002  % 
S6I  g 


<\J 


H06I 
S81 
08  T 
SZI 
OZI 
-  991 
-09T 
-SSI 
OS  I 
S*J 
0*1 
SET 
0£T 
S2I 
02T 
STT 

on 

SOI 
00T 


CD 


0(0(MOD^tO(D(\JGD'^fO(OC\JCD^rO 

(d  33a-2-JJ-UH/ni8)  M 

Figure  48.   Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  Left  Hand  Heater,  6.25  in. 
Configuration 
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Figure  49.   Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  8.0  in.  Configuration 
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Figure  50.   Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  10.0  in.  Configuration 
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Figure  51.   Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  12.0  in.  Configuration 
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Figure  52.   Heat  Transfer  Coefficient  vs  Superficial  Mass 

Velocity;  Right  Hand  Heater,  All  Configurations 
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Figure  53.  -Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  Left  Hand  Heater,  All  Configurations 
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A  second  significant  trend  which  was  observed  was  the 
fact  that  the  heat  transfer  coefficient  for  the  left  hand 
heater  was  significantly  higher  than  for  the  right  hand 
heater,  particularly  as  the  bed  width  was  increased.   In 
the  wide  bed  configurations  the  heat  transfer  coefficient 
for  the  right  hand  heater  showed  very  little  improvement 
in  response  to  an  increase  in  flow  rate.   With  a  high 
height-to-width  ratio  (narrow  bed)  ,  the  heat  transfer 
coefficient  for  the  left  hand  heater  was  roughly  twice 
the  value  for  the  right  hand  heater  at  high  gas  flow  rates. 
With  a  low  height-to-width  ratio  (wide  bed)  ,  the  left  hand 
heat  transfer  coefficient  exceeded  that  for  the  right  by 
a  factor  of  15  to  1 . 

The  exception  to  this  trend  occurred  during  the  data 
runs  accomplished  with  the  distributor  masked-off.   Under 
this  condition,  the  heat  transfer  coefficients  for  the  left 
hand  and  right  hand  heaters  were  in  much  closer  agreement 
with  one  another  than  in  all  other  cases.   The  effect  of 
masking  the  distributor  was  to  increase  the  heat  transfer 
coefficient  for  the  right  hand  heater  by  approximately  10% 
while  simultaneously  reducing  the  heat  transfer  coefficient 
for  the  left  hand  heater  by  approximately  25%.   It  would  be 
expected  that,  with  the  distributor  masked-off,  the  heat 
transfer  coefficients  would  be  the  same  for  both  the  left 
hand  and  right  hand  heaters.   Although  there  is  a  much 
closer  match  than  in  the  unmasked  configuration,  it  is 
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observed  that  the  left  hand  heater  still  achieved  higher  rates 
of  heat  transfer.   There  are  several  possible  explanations 
for  this  difference.   The  first  of  these  is  that  the  unused 
portion  of  the  distributor  may  not  have  been  completely  sealed 
by  the  masking  and  as  a  result,  some  peripheral  injection 
could  have  still  have  been  taking  place.   Secondly,  it  was 
noted  that  there  was  a  small  amount  of  gas  leakage  past  the 
movable  side  wall  gaskets  which  may  have  altered  the  gas  and 
particle  flow  patterns  along  the  left  wall.   This,  in  turn, 
may  have  altered  the  heat  transfer  coefficients.   Thirdly, 
it  is  possible  that  the  movable  side  wall  may  not  have  been 
perfectly  vertical.   Although  care  was  taken  to  accurately 
position  and  secure  the  side  wall,  deviations  from  a  true 
vertical  orientation  of  as  much  as  1°  are  considered  possible. 
Research  conducted  in  the  Soviet  Union  by  Filippovskii  and 
Baskakov  [Ref.  7]  has  indicated  that  plate  angle  has  a  con- 
siderable influence  on  heat  transfer  coefficients.   Because 
of  the  small  angle  involved  in  this  case,  the  contribution 
of  this  effect  is  not  considered  to  be  the  principle  cause 
of  the  difference  in  performance  between  the  two  plates. 
A  final  effect  which  may  have  made  a  minor  contribution  to 
this  difference  is  the  fact  that  the  losses  from  the  heater 
backing  to  the  atmosphere  were  lower  for  the  left  hand  heater 
than  for  the  right.   As  a  result,  more  of  the  electrical 
energy  going  into  the  heater  was  directed  through  the  copper 
plate  and  on  to  the  bed,  thus  resulting  in  a  minor  improvement 
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in  the  heat  transfer  coefficient.   Because  of  the  confisura- 
tion  of  the  apparatus,  the  left  hand  heater  mounting  block 
was  enclosed  within  the  unused  portions  of  the  apparatus 
walls.   Because  of  this,  there  was  a  stagnation  of  the  air 
in  the  vicinity  of  the  mounting  block  as  well  as  an  elevation 
of  the  local  ambient  temperature.   These  effects  combined  to 
reduce  the  amount  of  heat  lost  through  the  mounting  block. 
Because  of  the  specialized  conditions  which  were  found  to 
exist  near  the  left  hand  wall,  the  performance  of  the  right 
hand  heater  shall  be  considered  as  being  representative  of 
a  typical  flat  plate. 

The  magnitude  of  the  heat  transfer  coefficients  calculated 
as  a  result  of  this  study  were  within  the  range  of  values 
found  in  the  literature,  however  it  should  be  noted  that 
there  is  a  wide  variation  in  the  coefficients  found  by 
other  researchers  due  to  the  effect  of  differing  bed  sizes 
and  geometries,  particle  size  and  composition,  fluidizing 
gas  properties,  superficial  mass  velocities,  and  bed  operat- 
ing temperatures. 

There  is  also  a  wide  variety  of  empirical  correlations 
available  for  the  prediction  of  heat  transfer  coefficients, 
however  each  of  these  serves,  for  the  most  part,  only  to 
describe  the  particular  test  apparatus  and  conditions  from 
which  they  were  derived.   For  this  reason,  no  attempt  was 
made  to  fit  the  data  to  the  correlations  currently  available 
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nor  was  an  attempt  made  to  develop  a  new  correlation  to 

describe  the  performance  of  the  apparatus  used  in  this 

experiment . 

A  review  of  the  literature  shows  that  numerous  researchers 

have  observed  heat  transfer  curves  similar  to  that  shown  in 

Figure  54.   This  typical  curve  displays  a  peak  heat  transfer 

coefficient  (h    )  at  some  optimum  superficial  mass  velocity 
max  r        r 

(G   .).   Flow  rates  above  this  optimum  value  resulted  in  a 
reduction  in  the  heat  transfer  coefficient.   This  has  been 
attributed  to  the  fact  that  as  the  voidage  is  increased, 
there  is  a  point  at  which  the  beneficial  contribution  of 
increased  particle  activity  is  offset  by  the  negative  effect 
of  reducing  the  solids  concentration  within  the  emulsion. 
In  this  study  such  an  optimum  flow  rate  was  not  achieved. 
It  is  believed  that,  due  to  the  design  of  the  fluidization 
apparatus ,  the  range  of  gas  flow  rates  was  restricted  by 
the  relatively  shallow  "freeboard"  region  above  the  level 
of  the  settled  bed.   When  the  bed  was  undergoing  active 
fluidization,  particles  within  the  "splash  zone"  reached 
the  top  of  the  container  and  blocked  the  air  exit  ports . 
This  was  believed  to  be  taking  place  at  gas  flow  rates 
below  the  optimum  level.   It  is  possible  that  a  reduction 
in  the  settled  bed  depth  and  a  consequent  increase  in  the 
height  of  the  freeboard  region  would  permit  higher  super- 
ficial mass  velocities  and  result  in  higher  levels  of  heat 
transfer . 
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Figure  54.  ,  Heat  Transfer  Coefficient  vs  Superficial  Mass 
Velocity;  Typical 
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The  difference  in  the  level  of  heat  transfer  between 
the  left  and  right  hand  walls  is  clearly  related  to  the 
high  level  of  cooling  observed  at  the  left  hand  wall. 
This,  in  turn,  is  related  to  the  large  amount  of  particle 
activity  observed  to  the  left  of  center.   The  exact  mech- 
anism of  heat  transfer  which  takes  place  at  the  wall  cannot 
be  stated  with  certainty.   Although  there  are  many  theories 
as  to  the  mechanism  of  heat  transfer  within  a  fluidized 
bed,  the  bulk  of  the  models  postulated  to  date  fall  into 
one  of  two  broad  categories.   The  first  of  these  models 
proposes  that  the  heat  transfer  is  primarily  through  a  gas 
film  adjacent  to  the  heat  transfer  surface.   The  heat  is 
transferred  by  steady  state  conduction  through  the  film  to 
the  downward  moving  particles.   This  film  is  kept  thin  by 
the  scrubbing  action  of  the  solid  particles  which  come  in 
contact  with  the  surface.   This  reduction  in  the  boundary 
layer  thickness  is  thought  to  be  the  key  to  achieving  high 
rates  of  heat  transfer.   The  second  model  theorizes  that 
heat  transfer  is  primarily  by  unsteady  conduction  to  "packets 
which  are  made  up  of  several  solid  particles  grouped  together 
for  a  short  period  of  time.   These  packets  are  thought  to 
be  in  contact  with  the  heat  transfer  surface  for  only  brief 
periods  of  time  before  they  move  away  from  the  surface  and 
are  replaced  or  "replenished"  by  fresh  packets.   As  a  result, 
there  is  always  a  high  temperature  differential  between  the 
surface  and  the  packets  and  this,  in  turn,  yields  high  levels 
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of  heat  transfer.   A  significant  factor  shared  by  both 
the  packet  model  and  the  film  conduction  model  is  that 
they  both  rely  on  the  free  motion  of  particles  to  effect 
high  rates  of  heat  transfer.   It  is  believed  that  in  this 
apparatus  a  high  level  of  particle  activity  near  the  wall 
promotes  the  frequent  exchange  of  particles  from  within  the 
bed  core  where  it  is  comparatively  cool.   As  the  frequency 
of  particle  or  packet  replenishment  increases,  the  tempera- 
ture differential  between  the  wall  and  the  particles  also 
increases.   This,  in  turn,  increases  the  overall  heat  transfer 
Conversely,  when  there  is  comparatively  little  particle  activ- 
ity near  the  wall,  the  replenishment  rate  is  low  and  the  same 
particles  remain  in  contact  with  the  hot  wall  for  longer 
periods  of  time  as  they  move  downward.   Because  of  the  long 
residence  time,  the  temperature  differential  is  diminished 
and  the  heat  transfer  coefficient  is  reduced.   In  this  study, 
the  level  of  particle  activity  near  the  right  hand  wall  was 
comparatively  low  and  so  was  the  heat  transfer  coefficient. 
A' study  of  the  temperature  profiles  for  the  right  hand  plate 
shows  a  steady  increase  in  plate  temperature  with  a  downward 
change  in  vertical  position.   This  profile  is  consistent 
with  that  expected  for  long  particle  residence  times  since 
the  downward  moving  particles  would  experience  a  steady  in- 
crease in  temperature  as  they  approached  the  bottom  of  the 
plate.   It  is  noted  that  when  the  configuration  of  the  bed 
is  tall  and  narrow,  the  bubble  flow  path  is  very  close  to 

114 


the  bed  walls.   This  causes  a  high  degree  of  particle  activity 
near  the  wall  and  a  corresponding  increase  in  the  particle 
replenishment  rate.   As  the  bed  height-to-width  ratio  is 
decreased  by  increasing  the  bed  width,  the  bubble  flow  path 
moves  further  away  from  the  wall.   This  causes  a  reduction 
in  the  amount  of  activity  near  the  wall  and  a  corresponding 
decrease  in  the  replenishment  rate.   This  scenario  is  con- 
sistent with  the  observed  reduction  in  heat  transfer  coef- 
ficient as  the  bed  width  was  increased.   Conversely,  the 
amount  of  activity  seen  by  the  left  hand  plate  was  consis- 
tently high  as  a  result  of  the  peripheral  injection  effect 
described  previously.   The  degree  of  activity  and  consequently 
the  heat  transfer  coefficient,  remained  unaffected  by  an  in- 
crease in  bed  width.   In  addition,  the  heater  temperature 
profile  showed  no  sharp  transition  from  top  to  bottom,  thus 
indicating  that  the  particle  residence  time  was  short. 
Based  on  the  findings  of  this  study,  it  may  be  concluded 
that  peripheral  injection  is  an  effective  means  of  promoting 
particle  activity  or  turbulence  along  the  wall  of  the  con- 
tainer.  The  idea  of  turbulence  promoters  in  fluidized  beds 
is  not  a  new  one.   Early  research  by  the  Babcock  and  Wilcox 
Company  [Ref.  8]  included  the  use  of  spirally  shaped  strips 
to  create  a  rotational  motion  which  forced  the  solids  toward 
the  walls  by  centrifugal  force.   The  result  was  a  significant 
increase  in  heat  transfer.   Recently  Sokolov  et  al  [Ref.  9] 
in  the  Soviet  Union  studied  the  use  of  angled  louvers 

115 


positioned  within  the  bed  in  order  to  redirect  the  motion  of 
the  particles  outward  toward  the  walls,  thus  destroying  the 
stagnant  layer  which  was  insulating  the  surface.   The  result 
was  a  marked  improvement  in  heat  transfer  performance. 
Although  more  study  is  obviously  warranted,  the  use  of  pe- 
ripheral injection  would  appear  to  offer  benefits  similar 
to  other  enhancement  techniques  while  eliminating  the  need 
for  in-bed  structures  which  occupy  valuable  space  and  which 
require  periodic  replacement  as  a  consequence  of  particle 
errosion.   In  addition,  the  amount  of  turbulence  created  at 
the  wall  by  the  injection  technique  could  easily  be  adjusted 
to  optimum  levels  by  controlling  the  flow  rate  of  the  in- 
jection gas.   Such  a  system  could  be  used  to  offset  the  low 
rates  of  heat  transfer  at  the  walls  of  wide  fluidized  beds. 

As  discussed  previously,  an  energy  balance  was  conducted 
in  order  to  confirm  the  validity  of  the  numerical  results. 
The  results  of  this  energy  balance  were  disappointingly 
poor.   In  most  cases,  the  measured  heat  flux  leaving  the 
bed  with  the  outgoing  airstream  was  found  to  be  only  half 
of  that  calculated  to  be  entering  the  bed  via  the  electrical 
heaters.   It  is  believed  that  the  disagreement  was  primarily 
due  to  incorrect  measurement  of  the  fluidizing  gas  flow 
rate.   Despite  the  poor  energy  balance,  the  conclusions 
reached  as  a  result  of  this  study  are  considered  to  be  valid 
since  they  were  based  on  observed  trends  rather  than  on  the 
magnitude  of  specific  numerical  results. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

Based  on  the  experimental  results  of  this  study  the 
following  conclusions  may  be  drawn: 

1.  The  plate-to-bed  heat  transfer  coefficient  for  a  flat 
vertical  plate  serving  as  a  containing  wall  decreases  with 
a  decrease  in  the  bed  height-to-width  ratio. 

2.  The  plate-to-bed  heat  transfer  coefficient  increases 
with  an  increase  in  superficial  mass  velocity  above  the  mini- 
mum required  to  initiate  fluidization. 

3.  The  plate-to-bed  heat  transfer  coefficient  remains  small 
and  relatively  constant  at  superficial  mass  velocities  below 
the  minimum  required  to  initiate  fluidization. 

4.  The  plate-to-bed  heat  transfer  coefficient  increases 
with  an  increase  in  the  particle  activity  or  turbulence  in 
the  vicinity  of  the  wall. 

5.  The  particle  turbulence  near  the  vertical  plate,  and 
consequently  the  heat  transfer  coefficient,  may  be  enhanced 
by  the  injection  of  air  from  beneath  the  plate  in  a  direction 
perpendicular  to  the  primary  fluidizing  air  flow  path. 

6.  The  temperature  profile  of  a  vertical  flat  plate  may  be 
used  to  provide  insights  as  to  the  nature  of  particle  activity 
in  the  vicinity  of  the  plate. 


117 


7.  The  particle  flow  pattern  within  a  rectangular  fluidized 
bed  is  in  very  close  agreement  with  that  described  in  the 
literature  for  cylindrical  fluidized  beds  with  upward  motion 
of  the  particles  taking  place  along  the  centerline  of  the 
container  and  downward  return  flow  of  the  particles  taking 
place  along  the  periphery. 

8.  Slugging  or  heaving  of  the  bed  was  more  prevalent  in  the 
narrow  bed  configuration  as  opposed  to  the  wide  bed  configura 
tion  where  this  phenomenon  was  not  observed. 

9.  Fluidized  bed  particle  motion  along  a  non-moving  boundary 
matches  the  behavior  of  a  viscous  fluid  in  that  a  "no-slip" 
condition  is  observed. 

10.  The  bed  pressure  drop  characteristics  for  a  rectangular 
fluidized  bed  match  those  found  in  the  literature  for  a  cy- 
lindrical fluidized  bed  in  that  the  pressure  drop  across  the 
bed  increases  linearly  with  an  increase  in  air  flow  rate 
until  the  point  of  fluidization  onset.   Once  fluidization 
has  commenced,  the  pressure  drop  remains  constant. 

11.  The  temperature  of  the  heater  surface  at  a  given  depth 
remains  constant  for  any  point  in  a  horizontal  direction. 

12.  The  bed  temperature  will  assume  a  uniform  temperature 
distribution  once  fluidized. 

B.   RECOMMENDATIONS  FOR  FURTHER  STUDIES 

In  order  to  better  pursue  the  ojbectives  of  this  study  at 
sometime  in  the  future,  the  following  recommendations  and/or 
suggestions  are  made: 
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1.  Replace  the  three  circular  air  exhaust  tubes  on  the  top 
of  the  f luidizat ion  apparatus  with  a  rectangular  exhaust 
port  which  runs  the  entire  twenty-four  inch  length  of  the 
bed.   This  would  provide  a  non- restrict ive  exit  path  for 
the  fluidizing  air  which  would  minimize  the  amount  of  swirl- 
ing and  mixing  in  the  region  above  the  bed. 

2.  Repair  severed  wire  on  thermocouple  #15. 

3.  Replace  malfunctioning  thermocouples  #4  and  *40. 

4.  Install  a  minimum  of  four  new  thermocouples  to  monitor 
air  exit  temperature.   At  present  only  one  thermocouple  can 
be  used  for  this  purpose  (T(55))  for  bed  configurations  of 
12  inches  or  less. 

5.  Provide  a  more  positive  means  of  securing  the  baffle 
plate  to  the  distributor  plate  so  as  to  prevent  leak-by  into 
the  bed. 

6.  Install  a  bypass  valve  around  the  existing  air  inlet  con- 
trol valve  to  permit  fine  adjustment  of  fluidization  flow  rates 

7.  Install  two  additional  thermocouples  in  the  left  hand 
movable  wall  1/4  inch  below  the  bottom  of  the  heater  in 
order  to  obtain  a  more  accurate  estimate  of  wall  losses. 

8.  Install  two  additional  thermocouples  in  the  front  wall 
of  the  bed  so  as  to  obtain  a  more  accurate  estimate  of  wall 
losses . 

9.  Substitute  the  existing  power  supplies  with  units,  which 
are  capable  of  providing  additi-onal  current  and  which  are 
equipped  with  controls  to  vary  power  as  well  as  voltage. 
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This  would  permit  equalization  of  input  power  to  the  two 
heaters  and  allow  the  desired  temperature  differential 
between  the  plate  and  the  bed  to  be  maintained. 

10.  Install  a  movable  wall  in  the  air  plenum  beneath  the 
distributor  plate  so  that  the  flow  path  followed  by  the  in- 
coming air  is  identical  for  the  left  and  right  sides. 

11.  Prior  to  the  commencement  of  future  data  runs,  apply 
caulking  material  to  the  movable  side  wall  so  as  to  prevent 
leakage  of  fluidizing  air  past  the  gaskets. 

12.  Modify  the  fluidization  apparatus  so  that  a  controlled 
amount  of  injection  air  may  be  applied  to  each  of  the  heater 
plates.   Such  a  modification  would  permit  additional  study 
of  the  effects  of  the  peripheral  injection  phenomenon. 

C.   CLOSING  REMARKS 

Gas  fluidization  is  a  fascinating  process  which  promises 
to  provide  a  multitude  of  beneficial  applications.   Full 
enjoyment  of  the  potential  which  exists  depends  on  a  more 
complete  understanding  of  the  many  complex  phenomenon 
associated  with  fluidization.   It  is  hoped  that  this  thesis 
has,  in  some  small  way,  contributed  to  that  understanding 
or  helped  to  stimulate  additional  thought  and  research  in 
this  exciting  field  of  study. 
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APPENDIX  A. 
EXPERIMENTAL  PRESSURE  DROP  DATA 
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APPENDIX  B. 
EXPERIMENTAL  ROTOMETER  CALIBRATION  DATA 


**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

20  PERCENT  READING 


**EXPER I MENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  -  72.0  DEG  F 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =  10.70  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  =  .010  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  =  .334 


**CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  = 
VENTURI  INLET  GAGE  PRESSURE  = 
VENTURI  INLET  ABSOLUTE  PRESSURE  = 
VENTURI  DIFFERENTIAL  PRESSURE  = 
VENTURI  INLET  AIR  DENSITY  ■ 
VENTURI  PRESSURE  DROP  RATIO  (X)  = 
EXPANSION  FACTOR  <Y>  = 
ASSUMED  REYNOLDS  NUMBER  = 
DISCHARGE  COEFFICIENT  (C)  = 
VELOCITY  OF  APPROACH  FACTOR  (E)  = 
FLOW  COEFFICIENT  <K>  = 
AIR  MASS  FLOW  RATE  = 
AIR  VOLUMETRIC  FLOW  RATE  = 
VENTURI  INLET  REYNOLDS  NUMBER  = 


14.6259  PSI 

.38654  PSI 

14.23932  PSI 

.00030  PSI 

.07230  LBM/FT3 

. 00002 

1.00 

3000 

.889 

1 . 03280 

.91816 

.01021  LBM/SEC 

8.469  CFM 

3002 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

25  PERCENT  READING 


♦♦EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =   12.35  IN  H20 

VENTURI  DIFFERENTIAL  PRESSURE  =  .011  IN  RED  OIL 

SPECIFIC  GRAVITY  OF  RED  OIL  =  .834 


**CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  =  14.6259  PS I 

VENTURI  INLET  GAGE  PRESSURE  =  .44614  PS  I 

VENTURI  INLET  ABSOLUTE  PRESSURE  =      14.17972  PS I 

VENTURI  DIFFERENTIAL  PRESSURE  =  .00033  PS  I 

VENTURI  INLET  AIR  DENSITY  =  .07200  LBM/FT3 

VENTURI  PRESSURE  DROP  RATIO  (X)  =            .00002 

EXPANSION  FACTOR  <Y>  =  1.00 

ASSUMED  REYNOLDS  NUMBER  =  3150 

DISCHARGE  COEFFICIENT  (C)  =  .391 

VELOCITY  OF  APPROACH  FACTOR  (E)  =           1.03230 

FLOW  COEFFICIENT  (K)  =  .91970 

AIR  MASS  FLOW  RATE  =  .01070  LBM/SEC 

AIR  VOLUMETRIC  FLOW  RATE  =  3.916  CFM 

VENTURI  INLET  REYNOLDS  NUMBER  =  3147 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

30  PERCENT  READING 


**EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  =  72.0  DEG  P 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =   14.75  IN  H20 

VENTURI  DIFFERENTIAL  PRESSURE  =  .015  IN  RED  OIL 

SPECIFIC  GRAVITY  OF  RED  OIL  =  .834 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  =  14.6259  PS I 

VENTURI  INLET  GAGE  PRESSURE  =  .53284  PSI 

VENTURI  INLET  ABSOLUTE  PRESSURE  =      14.09302  PSI 

VENTURI  DIFFERENTIAL  PRESSURE  =  .00045  PSI 

VENTURI  INLET  AIR  DENSITY  =  .07156  LBM/FT3 

VENTURI  PRESSURE  DROP  RATIO  (X)  =             .00003 

EXPANSION  FACTOR  (Y)  =  1.00 

ASSUMED  REYNOLDS  NUMBER  =  3700 

DISCHARGE  COEFFICIENT  (C)  =  .895 

VELOCITY  OF  APPROACH  FACTOR  (E)  =            1.03280 

FLOW  COEFFICIENT  (K)  =  .92435 

AIR  MASS  FLOW  RATE  =  .01252  LBM/SEC 

AIR  VOLUMETRIC  FLOW  RATE  =  10.497  CFM 

VENTURI  INLET  REYNOLDS  NUMBER  =  3682 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

35  PERCENT  READING 


♦♦EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =  17.50  IN  H20 
VENTURI  DIFFERENTIAL  PRESSURE  =  .020  IN  RED  OIL 
SPECIFIC  GRAVITY  OF  RED  OIL  =  .834 


*  CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  =  14.6259  PS I 

VENTURI  INLET  GAGE  PRESSURE  =  .63219  PS  I 

VENTURI  INLET  ABSOLUTE  PRESSURE  =      13.99367  PSI 

VENTURI  DIFFERENTIAL  PRESSURE  -  .00060  PSI 

VENTURI  INLET  AIR  DENSITY  =  .07105  LBM/FT3 

VENTURI  PRESSURE  DROP  RATIO  (X)  =            .00004 

EXPANSION  FACTOR  <Y)  =  1.00 

ASSUMED  REYNOLDS  NUMBER  -  4250 

DISCHARGE  COEFFICIENT  <C>  =  .399 

VELOCITY  OF  APPROACH  FACTOR  (E)  =           1.03280 

FLOW  COEFFICIENT  (K)  =  .92797 

AIR  MASS  FLOW  RATE  =  .01446  LBM/SEC 

AIR  VOLUMETRIC  FLOW  RATE  =  12.211  CFM 

VENTURI  INLET  REYNOLDS  NUMBER.  =  4253 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

40  PERCENT  READING 


**EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =   20.30  IN  H20 

VENTURI  DIFFERENTIAL  PRESSURE  =  .025  IN  RED  OIL 

SPECIFIC  GRAVITY  OF  RED  OIL  =  .334 


**CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  =  14.6259  PSI 

VENTURI  INLET  GAGE  PRESSURE  =  .73334  PSI 

VENTURI  INLET  ABSOLUTE  PRESSURE  =      13.39252  PSI 

VENTURI  DIFFERENTIAL  PRESSURE  =  .00075  PSI 

VENTURI  INLET  AIR  DENSITY  =  .07054  LBM/FT3 

VENTURI  PRESSURE  DROP  RATIO  (X)  =             .00005 

EXPANSION  FACTOR  <Y)  =  1.00 

ASSUMED  REYNOLDS  NUMBER  -  4750 

DISCHARGE  COEFFICIENT  (C)  =  .901 

VELOCITY  OF  APPROACH  FACTOR  (E)  =            1.03230 

FLOW  COEFFICIENT  <K>  =  .93036 

AIR  MASS  FLOW  RATE  =  .01616  LBM/SEC 

AIR  VOLUMETRIC  FLOW  RATE  =  13.744  CFM 

VENTURI  INLET  REYNOLDS  NUMBER  =  4753 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

45  PERCENT  READING 


**EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  -  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =   24.30  IN  H2Q 

VENTURI  DIFFERENTIAL  PRESSURE  =  .031  IN  RED  OIL 

SPECIFIC  GRAVITY  OF  RED  OIL  =  .334 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  =  14.6259  PSI 

VENTURI  INLET  GAGE  PRESSURE  =  .87784  PSI 

VENTURI  INLET  ABSOLUTE  PRESSURE  =      13.74802  PSI 

VENTURI  DIFFERENTIAL  PRESSURE  =  .00093  PSI 

VENTURI  INLET  AIR  DENSITY  =  .06981  LBM/FT3 

VENTURI  PRESSURE  DROP  RATIO  (X)  =            .00007 

EXPANSION  FACTOR  (Y)  =  1.00 

ASSUMED  REYNOLDS  NUMBER  =  5300 

DISCHARGE  COEFFICIENT  <C>  =  .904 

VELOCITY  OF  APPROACH  FACTOR  <E)  =           1.03280 

FLOW  COEFFICIENT  <K)  =  .93365 

AIR  MASS  FLOW  RATE  =  .01795  LBM/SEC 

AIR  VOLUMETRIC  FLOW  RATE  =  15.431  CFM 

VENTURI  INLET  REYNOLDS 'NUMBER  =  5281 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

50  PERCENT  READING 


**EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =   27.30  IN  H20 

VENTURI  DIFFERENTIAL  PRESSURE  =  .037  IN  RED  OIL 

SPECIFIC  GRAVITY  OF  RED  OIL  =  .834 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  = 
VENTURI  INLET  GAGE  PRESSURE  - 
VENTURI  INLET  ABSOLUTE  PRESSURE  - 
VENTURI  DIFFERENTIAL  PRESSURE  = 
VENTURI  INLET  AIR  DENSITY  - 
VENTURI  PRESSURE  DROP  RATIO  (X)  - 
EXPANSION  FACTOR  (Y)  = 
ASSUMED  REYNOLDS  NUMBER  = 
DISCHARGE  COEFFICIENT  (C)  = 
VELOCITY  OF  APPROACH  FACTOR  <E)  = 
FLOW  COEFFICIENT  <K>  = 
AIR  MASS  FLOW  RATE  = 
AIR  VOLUMETRIC  FLOW  RATE  - 
VENTURI  INLET  REYNOLDS  NUMBER  = 


14.6259  PSI 

.98621  PSI 

13.63965  PSI 

.00111  PSI 

06926 

LBM/FT3 

. 00008 

1.00 

5750 

.906 

1 . 03280 

. 93592 

01958 

LBM/SEC 

16, 

.967  CFM 

5760 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

55  PERCENT  READING 


**EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =   31.85  IN  H20 

VENTURI  DIFFERENTIAL  PRESSURE  =  .044  IN  RED  OIL 

SPECIFIC  GRAVITY  OF  RED  OIL  =  .834 


♦♦CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  = 
VENTURI  INLET  GAGE  PRESSURE  = 
VENTURI  INLET  ABSOLUTE  PRESSURE  - 
VENTURI  DIFFERENTIAL  PRESSURE  = 
VENTURI  INLET  AIR  DENSITY  = 
VENTURI  PRESSURE  DROP  RATIO  (X)  = 
EXPANSION  FACTOR  <Y)  = 
ASSUMED  REYNOLDS  NUMBER  = 
DISCHARGE  COEFFICIENT  (C)  - 
VELOCITY  OF  APPROACH  FACTOR  (E)  = 
FLOW  COEFFICIENT  <K)  = 
AIR  MASS  FLOW  RATE  = 
AIR  VOLUMETRIC  FLOW  RATE  = 
VENTURI  INLET  REYNOLDS  NUMBER  = 


14.6259  PS I 

1. 15058  PS I 

13.47528  PS I 

.00133  PSI 

06842 

LBM/FT3 

.00010 

1  .  00 

6250 

.909 

1 . 03280 

. 93830 

02128 

LBM/SEC 

18. 

662  CFM 

6259 
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**ROTOMETER  CALIBRATION** 
ASME  HERSCHEL-TYPE  VENTURI 
4.26/2. 13  IN  DIA 

60  PERCENT  READING 


**EXPERIMENTAL  DATA** 


ATMOSPHERIC  PRESSURE  =  29.76910  IN  Hg 

VENTURI  INLET  TEMPERATURE  =  72.0  DEG  F 

VENTURI  INLET  STATIC  GAGE  PRESSURE  =   37.00  IN  H20 

VENTURI  DIFFERENTIAL  PRESSURE  =  .051  IN  RED  OIL 

SPECIFIC  GRAVITY  OF  RED  OIL  =  .334 


**CALCULATED  VALUES** 


ATMOSPHERIC  PRESSURE  =  14.6259  PS  I 

VENTURI  INLET  GAGE  PRESSURE  =  1.33663  PS I 

VENTURI  INLET  ABSOLUTE  PRESSURE  =      13.2S923  PSI 

VENTURI  DIFFERENTIAL  PRESSURE  =  .00154  PSI 

VENTURI  INLET  AIR  DENSITY  =  .06748  LBM/FT3 

VENTURI  PRESSURE  DROP  RATIO  (X)  =             .00012 

EXPANSION  FACTOR  (Y)  =  1.00 

ASSUMED  REYNOLDS  NUMBER  =  6700 

DISCHARGE  COEFFICIENT  (C)  =  .911 

VELOCITY  OF  APPROACH  FACTOR  <E>  =            1.03230 

FLOW  COEFFICIENT  <K)  =  .94036 

AIR  MASS  FLOW  RATE  =  .02230  LBM/SEC 

AIR  VOLUMETRIC  FLOW  RATE  =  20.277  CFM 

VENTURI  INLET  REYNOLDS  NUMBER  =  •     6707 


133 


APPENDIX  C.   EXPERIMENTAL  HEAT  TRANSFER  DATA 
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APPENDIX    D. 
SAMPLE    THERMOCOUPLE    CALIBRATION    DATA 

THERMOCOUPLE    #    14 
t  +  :fi  ************* 


TEMP 

PLAT 

TC 

BATH 

9£T 

PES 

>IST    0 

JTPUT 

TEMP 

**** 

****** 

**** 

1 

V  "? 

.1  ccc 

4  _i  j  _i 

"7  "7  "7 

13.  59 

£ 

■y  "7 

-_i  ._i  y ._, 

774 

1  9  .  6  1 

~7 

w  y 

1379 

1 

024 

-.  tr  r,   -. 

4 

hi y 

3  3  3  3 

1 

221 

3  0  '.  5  3 

5 

2  9 

4  5  6  9 

1 

549 

3  3  3  9 

K 

V  Q 

9557 

1 

759 

43  !  35 

f 

30 

7  ^  "7  ■  J 

ji 

107 

51  .  49 

C1 

31 

3107 

2 

340 

56  35 

Q 

J*  £, 

13  37 

"% 

i  vi  r 

65  .  33 

10 

--'  hi 

7447 

~. 

9  74 

7  1  .  20 

1 1 

"7  "7 

6364 

3 

401 

30  .  35 

12 

T  *7 

9  7  2  3 

7 

531 

83  53 

13 

77 

3729 

7 

■-■  tr  -7 

77  .50 

14 

"7  ~ 

4592 

~i 

344 

. 68 . 33 

15 

31 

9132 

"• 

605 

32  .  92 

16 

31 

04  3  1 

~i 

•-,  •-■  cr 

54  .  24 

17 

30 

r  ■-. .;.  -7 

~i 

0  0  2 

4  9  .  0  5 

18 

2  '~> 

707  0 

I 

.  6  5  2 

40  33 

19 

2  9 

2117 

1 

444 

35.9  6 

20 

w  y 

4275 

1 

121 

23  .  19 
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i* 
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,*-' 

r. 
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ai 

E-" 

jC 
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4J 

111 
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4  1     - 


34 


i!  b 


1  8 


1  7  - 
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3  4  - 

7  0  ■■ 
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.* 
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A"' 


** 


;*' 


,* 
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r- 


CVJ 


'."•J 

r- 


— f- 


—I 1 h- 

l*'J  ,-H  i—* 

IT'  <"H  .T, 
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—i     c-.j     ■■•j     '"'j     evi     ro     f-"' 
Thermocouple   Output    (mv) 
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*  *  *  P  0  L  Y  H  0  M  I  H  L    F  U  H  C  T  I  0  H  *  *  * 
(ORDER    =4) 


*X-1 


...    4 


*XA    2    + 
-  .  £1.3948987888 


*** 


774 

1 

62  4 

1 

221 

1 

549 

1 

759 

£ 

107 

2 

349 

2 

t   L_  1 

9  74 

"7 

401 

"7 

531 

T 

-■  C  -J 

£~ 

344 

,i 

605 

2 

"'  "■  c- 

•" 

0  0  2 

1 

652 

1 

444 

1 

121 

M 

i~1 

OF 

Y     (HOT;1 

******* 

13    59 

1  9  .  6  1 

-~i  C    .-i  i-i 

lL   J  .  O  lL 

30.53 
3  8  3  9 
43  .  35 
51  49 
56  .  85 
6  5  6  3 
71  .  £  0 
S  0    6  5 


.  ?  .  5ti 

6  8  3  3 
62  92 
54  24 
4  9  0  5 
4  0    8  8 


26  .  19 

SUM    OF    SQUARE    ERRORS* 

012263527709 


Y    ' 

C  H  L  C  > 

******** 

18 

5  9  6  6 

19 

6  198 

w  J 

7  9  0  3 

36 

5731 

"7  O 

4259 

43 

3  8  O  7 

5  1 

4394 

56 

854  9 

65 

6616 

71 

2135 

3  O 

67  17 

O  ~7 

5115 

i     i 

5034 

6  3 

2984 

o  4 

8  9  9  3 

54 

2127 

4  9 

M5  53 

4  0 

3  624 

"7  c* 

9  2  9  O 

w  O 

1554 
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APPENDIX  E. 
SAMPLE  HEAT  TRANSFER  DATA  REDUCTION  COMPUTER  PRINTOUTS 

1.   Non-f luidized  6.25  in.  configuration. 


♦^EXPERIMENTAL  DATA** 

RUN  #  06B20 

06/08/82  -  0425 


RUN  NUMBER  = 
BED  WIDTH  = 
STATIC  BED  HEIGHT  = 
AIR  FLOW  RATE  = 
AMBIENT  TEMPERATURE  = 
RIGHT  HEATER  VOLTAGE  = 
RIGHT  HEATER  CURRENT  = 
LEFT  HEATER  VOLTAGE  = 
LEFT  HEATER  CURRENT  - 
BED  EXPANSION  = 


06B20 

6.25  INCHES 

12.5   INCHES 

7.29  CFM 

74.50  DEG  F 

44.0  VOLTS 

.82  AMPS 

45.0  VOLTS 

1.00  AMPS 

0.0   INCHES 


COMMENTS:  NON-FLUIDI ZED/NON-EXPANDED 
BED  WITH  SMALL  LOCAL  SPOUT  AT  10-10 
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♦♦EXPERIMENTAL  DATA** 

RUN  #  06B20 

06/08/82  -  0425 

THERMOCOUPLE  READINGS 


TC# 

OUTPUT 

TEMP 

(mV) 

(DEG  P) 

1 

3.  332 

174.65 

2 

3.331 

174.55 

3 

_•> .  3o6 

174.76 

4 

.977 

79.23 

5 

3 .  330 

174.55 

6 

3 .  337 

174.31 

7 

3.343 

175.05 

8 

o  •  335 

174.71 

9 

3.312 

173.78 

10 

3.322 

174. 17 

11 

3.319 

174.07 

12 

3.316 

173.94 

13 

3.278 

172.36 

14 

3.282 

172.50 

15 

3.287 

172.72 

16 

3.286 

172.60 

17 

3.230 

1 70 . 36 

18 

170.55 

19 

3.238 

170.72 

20 

O  ■  £.  O  -i 

170.58 

21 

3.  129 

166.36 

3.  132 

166.51 

23 

3.127 

166.26 

24 

3.  129 

166.35 

25 

3.  130 

166.41 

26 

3.  131 

166.47 

27 

3.  128 

166.34 

28 

3.  128 

166.31 

29 

3.110 

165.59 

30 

3.113 

165.70 

31 

3.  Ill 

165.63 

32 

3.  108 

165.43 

•-J'O 

3.077 

164.22 

34 

3.075 

164. 14 

35 

3.075 

164. 11 

36 

3.  072 

163. 96 

37 

3 .  0 1  6 

161 .74 

33 

3 .  02 1 

161.93 

39 

3 .  0  i  9 

161.83 

40 

3 .  020 

162. 13 

TC# 

OUTPUT 

TEMP 

(mV) 

(DEG  P) 

41 

2.  110 

124.98 

42 

2.076 

123.54 

43 

1.958 

118.62 

44 

2.087 

123.87 

45 

2.  112 

124.92 

46 

2.002 

1 20 . 34 

47 

.911 

73.51 

48 

.911 

73.49 

49 

.910 

73.46 

50 

1  .  430 

96.22 

51 

1 .  408 

95.30 

52 

1.385 

94.33 

53 

.953 

75.34 

54 

.948 

75.  11 

55 

1  .  202 

86.  16 

56 

1.422 

95.35 

57 

1.417 

95.62 

68 

2.035 

121. 70 

69 

2.483 

140. 40 

70 

1.952 

113. 20 

71 

2.  Ill 

124.90 

72 

1 .  890 

115.58 

73 

2.072 

123.26 

74 

.936 

74.29 

75 

.  940 

74.47 

76 

2.  119 

125.24 

77 

2.331 

136. 13 

78 

1.872 

1  14.31 

79 

2.023 

121.41 
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♦♦EXPERIMENTAL  DATA** 
RUN#  06B20 
06/08/32  -  0425 


RIGHT-HAND  HEATER 
TEMPERATURE  PROPILE 


COL  1 


ROW 

1 

174. 

65 

ROW 

2 

174. 

55 

ROW 

O 

173. 

78 

ROW 

4 

172. 

36 

ROW 

170. 

36 

COL 


COL 


COL  4 


174.55  174.76  ****** 

174.81  175.05  174.71 

174.17  174.07  173.94 

172.50  172.72  172.60 

170.55  170.72  170.58 


AVE 
174.66 
174.78 
173.99 
172.54 
170. 55 


NOTE:     TEMPERATURES    SHOWN    ARE    IN    DEGREES    FAHRENHEIT 
AND    ARE    ARRANGED    AS    SEEN    FROM    WITHIN    THE    BED 


RIGHT-HAND     HEATER 
TEMP     vs     DEPTH 


t l0 

Z     9 

* 

8 

S  r 

-     6 

* 

£  5 
0)  4 

- 

* 

* 

> 

- 

* 

sa 

in 
go 

in 

CO 

in 

TEMP  (DEG  F) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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♦♦EXPERIMENTAL  DATA** 
RUN#  06B20 

06/08/32  -  0425 


LEFT-HAND  HEATER 
TEMPERATURE  PROFILE 

COL  1    COL  2    COL  3    COL  4        AVE 

ROW  1  !  166.36   166.51   166.26   166.35  !   166.37 

ROW  2  !  166.41   166.47   166.34   166.31  !   166.33 

ROW  3  !  165.59   165.70   165.63   165.43  !   165.59 

ROW  4  !  164.22   164.14   164.11   163.96  !   164.10 

ROW  5  !  161.74   161.93   161.83   162.13  !   161.91 


NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


z  1  0 


CO 

c 


OH 


LEFT-HnND  HEATER 
TEMP  v  s  D E  P T H 


in          s 

j-) 

m          u3 

■-3 

TEMP 

(  OE 

F  ) 


in 


a 

CD 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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♦♦EXPERIMENTAL  DATA^ 
RUN#  06B20 
06/08/82  -  0425 


PLUIDIZED  BED 
TEMPERATURE  PROFILE 


95.85    75.34    96.22 


95.  30 


95.62    75.11    94.33 


NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED 
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♦♦CALCULATED  RESULTS** 
RUN  #  06B20 
06/08/82  -  0425 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  =  173.30  DEG  F 

AVERAGE  PLATE  SURFACE  TEMPERATURE  =  173.28  DEG  F 

AVERAGE  BED  TEMPERATURE  =  89.68  DEG  F 

ELECTRICAL  ENERGY  INTO  HEATER  =  123.11  BTU/HR 

LOSS  FROM  HTR  BACKING  TO  ATM  =  11.22  BTU/HR 

h(HTR  BACKING-ATM)  =             .52  BTU/HR-FT^-DEG  F 

LOSS  FROM  HTR  PERIMETER  TO  BED  =  13.50  BTU/HR 

HEAT  FLUX  THRU  PLATE  =  98.39  BTU/HR 

h  (PLATE-BED)  =                  3.39  BTU/HR-FT"-2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  =  164.87  DEG  F 

AVERAGE  PLATE  SURFACE  TEMPERATURE  =      164.84  DEG  F 
AVERAGE  BED  TEMPERATURE  =  89.68  DEG  G 

ELECTRICAL  ENERGY  INTO  HEATER  =  153.54  BTU/HR 

LOSS  FROM  HTR  BACKING  TO  ATM  =  9.20  BTU/HR 

h(HTR  BACKING-ATM)  =  .42  BTU/HR-FT"2-DEG  F 

LOSS  FROM  HTR  PERIMETER  TO  BED  =  19.77  BTU/HR 

HEAT  FLUX  THRU  PLATE  =  124.58  BTU/HR 

h (PLATE-BED)  =  4.77  BTU/HR-FTN2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE  =  73.49  DEG  F 

AIR  OUTLET  TEMPERATURE  =  86.16  DEG  F 

AVERAGE  BED  TEMPERATURE  =  89.68  DEG  F 

q(AIR  OUT)  =  96.97  BTU/HR 

q(LOSS  F/R  WALL)  =  .50  BTU/HR 

q (TOTAL  OUT  OF  BED)  =  97.47  BTU/HR 

q (TOTAL  INTO  BED)  =  256.23  BTU/HR 

SUPERFICIAL  VELOCITY  =  .47  FT/SEC 
SUPERFICIAL  MASS  VELOCITY  =         121.92  LBM/HR-FT  2 

PARTICLE  REYNOLDS  NUMBER  =  2.44 
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Fluidized  6.24  in.  configuration 


♦♦EXPERIMENTAL  DATA** 

RUN  #  06B26 

06/08/82  -  0510 


RUN  NUMBER  = 
BED  WIDTH  = 
STATIC  BED  HEIGHT  = 
AIR  FLOW  RATE  ■ 
AMBIENT  TEMPERATURE  = 
RIGHT  HEATER  VOLTAGE  = 
RIGHT  HEATER  CURRENT  = 
LEFT  HEATER  VOLTAGE  = 
LEFT  HEATER  CURRENT  = 
BED  EXPANSION  = 


06B26 

6.25  INCHES 

12-5   INCHES 

9.23  CFM 

74.50  DEG  F 

60.0  VOLTS 

1. 12  AMPS 

60.0  VOLTS 

1.28  AMPS 

.1    INCHES 


COMMENTS:  FLUIDIZATION  WITH  2-3  INCH 
DIA  BUBBLES  ERUPTING  AT  LOCATION  5-5 
AND  VICINITY/ACTIVE  BOT  LEFT 
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**EXPERIMENTAL    DATA** 

RUN    #    0&B26 

06/08/82    -    0510 

THERMOCOUPLE    READINGS 


TC# 

OUTPUT 

TEMP 

(mV) 

(DEG    F) 

1 

2.729 

150.44 

2 

2.728 

150.36 

o 

2.724 

1 50 . 1 3 

4 

.961 

78.  19 

5 

155.53 

6 

2.859 

155.67 

7 

2.858 

155.63 

8 

2.848 

1  55  ■  jl  X 

9 

2.991 

1 60 . 95 

10 

3 .  000 

161.31 

11 

2.996 

161.  15 

12 

2.988 

1 60 . 35 

13 

3.  107 

165. 53 

14 

3 .  1 06 

165.48 

15 

3.  106 

165.49 

16 

3 .  1 06 

165.43 

17 

3.  141 

166.30 

18 

3.  144 

166.92 

19 

3.  140 

166.31 

20 

3 .  1 37 

166.79 

21 

130.99 

Tj 

2.290 

132.28 

2.291 

132.23 

24 

2.297 

132.51 

25 

2 .  340 

134.30 

26 

2.336 

134. 16 

27 

335 

134. 10 

28 

~y 

343 

134.63 

29 

"? 

346 

134.56 

30 

*7» 

335 

1 34 . 06 

31 

ry 

339 

134.26 

32 

T* 

354 

134.79 

.J ._ > 

2.293 

132.51 

34 

2.285 

131.94 

35 

2.291 

132. 13 

36 

2.310 

132.97 

37 

~?       '~>"T-r 

12^.73 

33 

2.213 

129. 16 

3^ 

129. 44 

40 

*■} 

259 

131.11 

TC# 

OUTPUT 

TEMP 

(mV) 

(DEG    F) 

41 

2.  146 

126.48 

42 

2.  116 

1 25 .21 

43 

1.996 

120.22 

44 

2.  160 

126.91 

45 

2.  154 

126.67 

46 

2 .  04 1 

121.97 

47 

.922 

74.  00 

48 

.922 

73.98 

49 

.922 

73.99 

50 

1.564 

101 .98 

51 

1   ccc 

101.63 

52 

1  .  55 1 

101.49 

53 

1.556 

101.66 

54 

1.588 

103. 00 

55 

1.526 

100. 12 

56 

1.551 

101. 40 

57 

1  .  550 

101. 35 

68 

1.927 

117. 14 

69 

1.958 

118.46 

70 

1 .  930 

117.27 

71 

2.363 

135.64 

72 

1.916 

116.63 

73 

2.  039 

121.37 

74 

1 .  045 

79.  10 

75 

1  .  260 

38.51 

76 

1.339 

113.41 

77 

1.356 

114. 14 

73 

1.673 

106.34 

79 

1.925 

117. 06 
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♦♦EXPERIMENTAL    DATA** 
RUN#    06B26 
06/08/32    -    0510 


RIGHT-HAND    HEATER 
TEMPERATURE    PROFILE 


COL    1 


COL 


COL    3 


COL    4 


AVE 


ROW 

1 

150. 

44 

ROW 

2 

155. 

53 

ROW 

o 

160. 

95 

ROW 

4 

165. 

53 

ROW 

5 

166. 

80 

150.36  150.18  ******  !  150.33 

155.67  155.63  155.21  !  155.51 

161.31  161.15  160.85  !  161.07 

165.48  165.49  165.43  !  165.48 

166.92  166.81  166.79  !  166.83 


NOTE:     TEMPERATURES    SHOWN    ARE    IN    DEGREES    FAHRENHEIT 
AND    ARE    ARRANGED    AS    SEEN    FROM    WITHIN    THE    BED 


RIGHT-HAND     HEATER 
TEMP     v  s     DEPTH 

0 
9 
8 
7h 


2 
c 
-     6 


CO 

o 


in 


s 


vJ3 


33 

N. 


LO 


TEMP     ( DEG     F ) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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♦♦EXPERIMENTAL    DATA^ 
RUN#    06B26 

06/08/82    -    0510 


ROW  1 

ROW  2 

ROW  3 

ROW  4 

ROW  5 


LEFT-HAND    HEATER 
TEMPERATURE    PROFILE 


COL    1 

130.99 

1 34 . 30 

134.56 

132.51 

129.78 


COL 


COL    3 


COL    4 


132.28  132.28  132.51 

134.16  134.10  134.63 

1 34 . 06  134.26  1 34 . 79 

131.94  132.18  132.97 

129.16  129.44  131.11 


AVE 
132.02 
1 34 . 30 
134.42 
132. 40 
129.37 


NOTE:     TEMPERATURES    SHOWN    ARE    IN    DEGREES    FAHRENHEIT 
AND    ARE    ARRANGED    AS    SEEN    FROM    WITHIN    THE    BED 


LEFT-HRND     HERTER 
TEMP     vs     DEPTH 


*-■ 

0 
9 

* 

8 

- 

o 

7 

% 

6 

1— 

5 

■ 

*   . 

CO 

4 

■ 

a. 

3 
2 

* 

i— 

i 

0 

* 

Ll! 

iS 

Li 

CSS 

ll 

<s 

D 

n> 

OJ 

CM 

ro 

n 

>* 

T 

T  F  r !  F   ( DEG  F  ) 


NOTE:  PLGTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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♦^EXPERIMENTAL  DATA** 
RUN#  06B26 
06/08/82  -  0510 


FLUIDIZED  BED 
TEMPERATURE  PROFILE 


101.40   101.66   101.98 


101.63 


101.35   103.00   101.49 


NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED 
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**CALCULATED  RESULTS** 
RUN  #  06B26 
06/08/32  -  0510 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  =  159.84  DEG  F 

AVERAGE  PLATE  SURFACE  TEMPERATURE  =      159.79  DEG  F 
AVERAGE  BED  TEMPERATURE  =  101.79  DEG  F 

ELECTRICAL  ENERGY  INTO  HEATER  =  229.29  BTU/HR 

LOSS  FROM  HTR  BACKING  TO  ATM  =  7.92  BTU/HR 

h(HTR  BACKING-ATM)  =  .36  BTU/HR-FT"2-DEG  F 

LOSS  FROM  HTR  PERIMETER  TO  BED  =  20.53  BTU/HR 

HEAT  FLUX  THRU  PLATE  =  200.79  BTU/HR 

h  (PLATE-BED)  =  9.97  BTU/HR-FT"V2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  =  132.60  DEG  F 

AVERAGE  PLATE  SURFACE  TEMPERATURE  =      132.54  DEG  F 
AVERAGE  BED  TEMPERATURE  =  101.79  DEG  G 

ELECTRICAL  ENERGY  INTO  HEATER  =  262.05  BTU/HR 

LOSS  FROM  HTR  BACKING  TO  ATM  =  1.54  BTU/HR 

h(HTR  BACKING-ATM)  =  .07  BTU/HR-FT"2-DEG  F 

LOSS  FROM  HTR  PERIMETER  TO  BED  =  33.49  BTU/HR 

HEAT  FLUX  THRU  PLATE  =  227.02  BTU/HR 

h (PLATE-BED)  =  21.26  BTU/HR-FT" 2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE  = 

AIR  OUTLET  TEMPERATURE  = 

AVERAGE  BED  TEMPERATURE  = 

q (AIR  OUT)  = 

q (LOSS  F/R  WALL)  = 

q (TOTAL  OUT  OF  BED)  = 

q (TOTAL  INTO  BED)  = 

SUPERFICIAL  VELOCITY  = 

SUPERFICIAL  MASS  VELOCITY  =         150, 

PARTICLE  REYNOLDS  NUMBER  = 


73 

,9<5 

»  DEG 

F 

1  00 

.  11 

:  DEG 

F 

101 

.7S 

>  DEG 

F 

247. 

04 

BTU/HR 

26. 

■re 

BTU/ 

HR 

*^  'f  .■ 

39 

BTU/ 

HR 

481  . 

33 

BTU/ 

HR 

. 

59 

FT/ 3 

EC 

65  LBM/HP-FT 

'■■  *■? 

152 


3.   N'on-  f  luidized  12.0  in.  configuration. 


♦^EXPERIMENTAL  DATA** 
RUN  #  12B30E 
06/17/82  -  0055 


RUN  NUMBER  = 
BED  WIDTH  - 
STATIC  BED  HEIGHT  = 
AIR  FLOW  RATE  = 
AMBIENT  TEMPERATURE  = 
RIGHT  HEATER  VOLTAGE  = 
RIGHT  HEATER  CURRENT  = 
LEFT  HEATER  VOLTAGE  = 
LEFT  HEATER  CURRENT  = 
BED  EXPANSION  = 


COMMENTS:  NON-FLUIDI ZED/NON-EXPANDED 
BED/NO  BUBBLES  BOTTOM  LEFT 


12B30E 

12.  OC 

)  INCHES 

12.5 

INCHES 

1C 

).53  CFM 

78.60  DEG  F 

46. 

0  VOLTS 

• 

85  AMPS 

56. 

0  VOLTS 

1. 

02  AMPS 

0.0 

INCHES 
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♦♦EXPERIMENTAL  DATA** 
RUN  #  12B30E 
06/17/82  -  0055 

THERMOCOUPLE  READINGS 


TC# 

OUTPUT 

TEMP 

(mV) 

(DEG  P) 

1 

3.  556 

183.47 

2 

3 .  550 

183.20 

3 

3.  553 

183.34 

4 

1 .  1 40 

89.84 

5 

3.554 

183.39 

6 

3.559 

183.55 

7 

3 .  565 

183.83 

8 

3.556 

183.41 

9 

3.541 

182.84 

10 

3.  551 

183.22 

11 

3.545 

183.03 

12 

3.541 

182.80 

13 

3.512 

181.62 

14 

3.515 

181.69 

15 

3.516 

181.78 

16 

3.  514 

181.61 

17 

3.468 

179.84 

18 

3.471 

179.90 

19 

3.468 

179.81 

20 

3.458 

179.52 

21 

2.886 

156.60 

22 

2.886 

1 56. 63 

23 

2.883 

156. 46 

24 

2.886 

156.58 

25 

2.  893 

156.89 

26 

2.894 

-  156.94 

27 

2.891 

156.82 

28 

2.890 

156.75 

29 

2.886 

156.59 

30 

2.888 

156.65 

31 

2.887 

156.62 

32 

2.886 

156.51 

2.869 

1 55 . 35 

34 

2.867 

1  wi^j .  76 

2.869 

155.81 

36 

2.866 

155.67 

37 

2.828 

154. 16 

38 

2.836 

154.47 

39 

1 54 . 30 

40 

2.831 

154.51 

TC# 

OUTPUT 

TEMP 

(mV> 

(DEG  P) 

41 

2.  123 

125.52 

42 

2.  145 

126.42 

43 

2 .  04 1 

122. 10 

44 

1.785 

1 1 1 . 20 

45 

1.862 

114.45 

46 

1.795 

111.65 

47 

.967 

75.99 

48 

.966 

75.92 

49 

.966 

75.94 

50 

1.438 

96.56 

51 

1.429 

96.21 

52 

1.403 

95.  11 

53 

1.273 

89.41 

54 

1.277 

89.56 

55 

1.257 

88.55 

56 

1.261 

88.86 

57 

1 .  260 

88 .  80 

58 

1.265 

89.01 

59 

1.274 

89.40 

60 

1.255 

88.87 

61 

1.257 

88.95 

62 

1.262 

89.  15 

63 

1.239 

88.  13 

68 

1.954 

118.29 

69 

2.477 

140. 15 

70 

2.  173 

127.50 

71 

2.527 

142.22 

72 

1.974 

119. 13 

73 

2 .  200 

128.63 

74 

1 .  095 

81.29 

75 

1.  165 

84.36 

76 

1.977 

119.26 

77 

2.  092 

124. 10 

78 

1.773 

110.61 

79 

1  .  740 

109. 20 
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**EXPER I MENTAL  DATA** 

RUN*  12B30E 

06/17/82  -  0055 


ROW  1 

ROW  2 

ROW  3 

ROW  4 

ROW  5 


RIGHT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1 


COL  2 


COL 


COL  4 


183.47  183.20  183.34  ****** 

183.39  183.55  183.83  183.41 

182.84  183.22  183.03  182.80 

181.62  181.69  181.78  181.61 

179.84  179.90  179.81  179.52 


AVE 

1 83 . 34 
1 83 . 55 
182.97 
181.68 
179.77 


NOTE:     TEMPERATURES    SHOWN    ARE    IN    DEGREES    FAHRENHEIT 
AND    ARE    ARRANGED    AS    SEEN    FROM    WITHIN    THE    BED 


RJGHT-HftND     HEfiTER 
TEMP     vs     DEPTH 


= 

\  0 

9 

* 

■z. 

D 

3 

7 

* 

►— < 

6 

t— 

5 

* 

CO 

4 

C' 

3 

* 

Li_ 

2 

- 

1— 

1 

* 

Q£ 

0 

Lu 

<S 

m 

s 

li'< 

S" 

in 

D> 

rv 

(S. 

CO 

id 

<Tl 

cr> 

TEMP  CDEG  F) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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**EXPER I MENTAL    DATA** 
RUN#    12B30E 

>055 


LEFT-HAND    HEATER 
TEMPERATURE    PROFILE 


CDL    1 


COL    2 


COL    3 


COL    4 


AVE 

ROW  1  !  156.60  156.63  156.46  156.58  !  156.57 
ROW  2  !  156. B9  156.94  156. 82  156.75  1  156.85 
ROW  3  !  156.59  156.65  156.62  156.51  !  156.59 
ROW  4  !  155.85  155.76  155.81  155.67  !  155.77 
ROW    5     !     154.16       154.47       154.30       154.51        !        154.36 


NOTE:     TEMPERATURES    SHOWN    ARE    IN    DEGREES    FAHRENHEIT 
AND    ARE    ARRANGED    AS    SEEN    FROM    WITHIN    THE    BED 


LEFT-HRNC     HERTER 
TEMP     ms     DEPTH 


X  i 

d 

►— * 

N-/ 

9 

8 

* 

•z. 

7 

* 

l»1 

6 

• 

■— 

5 

X 

*— * 

U*J 

4 

O 

3 

. 

% 

0. 

2 

■ 

y- 

1 

, 

* 

OL 

e 

, 

Ld 

\n 

IS 

in 

<St 

\T> 

Si 

Z> 

•*r 

m 

LO 

vu 

lO 

IV. 

TEMP  CDEG  F) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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^^EXPERIMENTAL  DATA** 

RUN#  12B30E 

06/17/82  -  0055 


FLUIDIZED  BED 
TEMPERATURE  PROFILE 


38.95    89.01    88.36    39.41    96.56 


89.15    89 . 40 


96.21 


88.13         88.87         88.80         89.56         95.11 


NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED 
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♦^CALCULATED  RESULTS** 
RUN  #  12B30E 
06/17/82  -  0055 


RIGHT-HAND  HEATER- 


AVERAGE  HEATER  TEMPERATURE  = 
AVERAGE  PLATE  SURFACE  TEMPERATURE  = 
AVERAGE  BED  TEMPERATURE  = 
ELECTRICAL  ENERGY  INTO  HEATER  = 
LOSS  FROM  HTR  BACKING  TO  ATM  = 


132.26  DEG  F 

182.23  DEG  F 

90.62  DEG  F 

133.41  BTU/HR 

12.85  BTU/HR 


h(HTR  BACKING-ATM)  = 


.61  BTU/HR-FT--2-DEG  F 


LOSS  FROM  HTR  PERIMETER  TO  BED  = 

HEAT  FLUX  THRU  PLATE  = 

h (PLATE-BED)  =  3.; 


14.77  BTU/HR 

105.80  BTU/HR 

BTU/HR-FT"2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  = 
AVERAGE  PLATE  SURFACE  TEMPERATURE  = 
AVERAGE  BED  TEMPERATURE  = 
ELECTRICAL  ENERGY  INTO  HEATER  = 
LOSS  FROM  HTR  BACKING  TO  ATM  = 


156.03  DEG  F 
155.98  DEG  F 

90.62  DEG  G 
194.90  BTU/HR 

9.90  BTU/HR 


h(HTR  BACKING-ATM)  = 


.63  BTU/HR-FT-2-DEG  F 


LOSS  FROM  HTR  PERIMETER  TO  BED  = 

HEAT  FLUX  THRU  PLATE  = 

h (PLATE-BED)  =  7.29 


19.57  BTU/HR 

165.43  BTU/HR 

BTU/HR-FT"2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE  = 

AIR  OUTLET  TEMPERATURE  = 

AVERAGE  BED  TEMPERATURE  = 

q (AIR  OUT)  = 

q (LOSS  F/R  WALL)  = 

q (TOTAL  OUT  OF  BED)  = 

q (TOTAL  INTO  BED)  = 

SUPERFICIAL  VELOCITY  = 

SUPERFICIAL  MASS  VELOCITY  = 

PARTICLE  REYNOLDS  NUMBER  = 


91 


75.95  DEG  F 

88.55  DEG  F 

90.62  DEG  F 

138.71  BTU/HR 

16.49  BTU/HR 

155.20  BTU/HR 

305.56  BTU/HR 

.35  FT/SEC 

38  LBM/HR-FT'2 

1  .83 
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Fluidized  12.0  in.  configuration. 


♦♦EXPERIMENTAL  DATA** 
RUN  #  12B48E 
06/17/82  -  0155 


RUN  NUMBER  = 
BED  WIDTH  = 
STATIC  BED  HEIGHT  = 
AIR  FLOW  RATE  - 
AMBIENT  TEMPERATURE  = 
RIGHT  HEATER  VOLTAGE  = 
RIGHT  HEATER  CURRENT  = 
LEFT  HEATER  VOLTAGE  = 
LEFT  HEATER  CURRENT  = 
BED  EXPANSION  = 


12B48E 

12.00  INCHES 

12-5   INCHES 

16.37  CFM 

81.50  DEG  F 

46.0  VOLTS 

.85  AMPS 

46.0  VOLTS 

1.02  AMPS 

0. 1   INCHES 


COMMENTS:  FLUIDIZATION  WITH  2-3  INCH 
BUBBLES  ERUPTING  AT  LOCATION  4-4 
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♦♦EXPERIMENTAL  DATA** 
RUN  #  12B48E 
06/17/82  -  0155 


THERMOCOUPLE 

READINGS 

TC# 

OUTPUT 

TEMP 

TC# 

<mV> 

(DEG  P) 

1 

3.471 

1 80 .  1 4 

41 

2 

3.474 

180.21 

42 

3 

3.478 

180.39 

43 

4 

1 .  083 

86.  14 

44 

5 

3.472 

180. 16 

45 

6 

3.479 

1  SO . 4 1 

46 

7 

3.484 

1 80 . 63 

47 

3 

3.476 

180.27 

48 

9 

3.445 

179.05 

49 

10 

3.455 

179.43 

50 

11 

3.450 

179.27 

51 

12 

3.445 

179.04 

52 

13 

3.394 

176.96 

53 

14 

3.399 

177. 13 

54 

15 

3.398 

177. 12 

55 

16 

3.  393 

176.84 

56 

17 

3 .  338 

174.67 

57 

18 

3.  343 

174.84 

58 

19 

3.339 

174.72 

59 

20 

3.327 

174.35 

60 

21 

2.831 

154.38 

61 

22 

2 .  833 

154.49 

62 

23 

2.833 

154.44 

63 

24 

2.839 

154.68 

68 

25 

2.866 

155.80 

69 

26 

2.866 

155.81 

70 

27 

2.864 

155.73 

71 

28 

2.864 

155.70 

72 

29 

2.869 

1 55 . 90 

73 

30 

2.867 

1 55 . 80 

74 

31 

2.863 

155.65 

75 

32 

2.863 

1 55 . 58 

76 

33 

2.833 

1 54 . 39 

77 

34 

2.827 

154. 14 

78 

35 

2.821 

153.87 

79 

36 

2.817 

153. 69 

37 

2.766 

151 . 65 

38 

2.764 

151.55 

39 

2.756 

151.22 

40 

2.755 

151.43 

OUTPUT     TEMP 

(mV)  (DEG  F) 

2.279  132.01 

2.244  130.53 

2.097  124.44 

2.121  125.29 
2.136  125.92 
2.024  121.26 

.959  75.64 

.959  75.62 

.960  75.67 

1.073  80.65 

1.122  82.84 
1.044  79.41 

.989  76.94 
1.089  81.33 
1.128  82.94 
1.097  81.67 
1 . 1 00  81.79 
1 . 1 00  81.79 
1.105  82.01 
1.096  81.90 
1.167  85.01 
1.096  31.37 
1.095      31.82 

2.133  125.82 

2.600  145.23 

2.088  123.93 

2.571  144.03 

1.946  117.95 

2.063  122.38 
1.031  78.48 
1.021      73.04 

2.024  121.24 

2.246  130.55 

1 . 306  112.01 

2.033  121.62 
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♦♦EXPERIMENTAL  DATA** 

RUN#  12B48E 

06/17/82  -  0155 


RIGHT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1 


COL 


COL  3 


COL  4 


ROW  1 

ROW  2 

ROW  3 

ROW  4 

ROW  5 


180.14  180.21 

180.16  180.41 

179.05  179.43 

176.96  177.13 

174.67  174.84 


180.39  ♦♦♦♦♦♦ 

180.63  180.27 

179.27  179.04 

177.12  176.84 

174.72  174.35 


AVE 
180.24 
1 80 . 37 
179.20 
177.01 
174.65 


NOTE:     TEMPERATURES    SHOWN    ARE    IN    DEGREES    FAHRENHEIT 
AND    ARE    ARRANGED    AS    SEEN    FROM    WITHIN    THE    BED 


R  I  GHT-HBHD     HERTER 
TEMP     v  s     CEPTH 


1  0 

9 

* 

o 

8 
7 

" 

* 

•— ■ 4 

6 

■ 

NM 

5 

• 

* 

LA 

4 

O 

3 
2 

♦ 

1 

e 

- 

* 

LU 

in 

s 

in 

ca 

in 

CO 

Z> 

*£> 

(^ 

IN. 

CO 

CO 

ETi 

TEMP  (DEG  F) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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**EXPER I MENTAL  DATA** 

RUN#  12B4BE 

'.'o/  i  /  /  a*.  —  Ol  jj 


LEFT-HAND  HEATER 
"EMPERATURE  PROFILE 


COL  1 


COL  2 


ROW  1  !  154.38  154.49 

ROW  2  I  155.30  155. SI 

ROW  3  !  155.90  155.80 

ROW  4  !  154.39  154. 14 

ROW  5  !  151.65  151.55 


COL  3 
154. 44 
1 55 . 73 
1 55 . 65 
153.87 
151.22 


COL  4 
154. 68 
155.70 
1 55 . 58 
153.69 
151.43 


c.i  ir 

r-ivc 

154. 50 
1 55 . 76 
155.73 
154.02 
151.46 


NOTE:     TEMPERATURES    SHOWN    ARE    IN    DEGREES    FAHRENHEIT 
AND    ARE    ARRANGED    AS    SEEN    FROM    WITHIN    THE    BED 


LEFT-HAND     HE  RTER 
TEMP     vs      DEPTH 


z  1 

0 

■ 

\^ 

9 

* 

o 

8 

7 

■ 

* 

~ 

6 

5 

* 

tn 

4 

o 

a. 

3 

2 

■ 

* 

i— 

1 

0 

* 

UJ 

in 

<s 

in 

s 

in 

S 

Z> 

T 

in 

in 

03 

<D 

(V. 

TEMP 

CDEG 

F  ) 

"" ' 

""" ' 

NOTE:     PLOTTED    TEMPERATURES    ARE    AVERAGES    OF    EACH 
HORIZONTAL    ROW 
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♦^EXPERIMENTAL  DATA** 

RUN#  12B48E 

06/17/82  -  0155 


FLUIDIZED  BED 
TEMPERATURE  PROFILE 


S5.01  81.79  81.67  76.94  80.6! 


81.87         82.01 


82.84 


81.82         81.90         81.79         81.33         79.41 


NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED 
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♦♦CALCULATED  RESULTS** 

RUN  #  12B48E 

06/17/S2  -  0155- 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  =  178.29  DEG  F 

AVERAGE  PLATE  SURFACE  TEMPERATURE  =  178.27  DEG  F 

AVERAGE  BED  TEMPERATURE  =  81.46  DEG  F 

ELECTRICAL  ENERGY  INTO  HEATER  =  133.41  BTU/HR 

LOSS  FROM  HTR  BACKING  TO  ATM  =  10.82  BTU/HR 

h(HTR  BACKING-ATM)  =             .51  BTU/HR-FT~ 2-DEG  F 

LOSS  FROM  HTR  PERIMETER  TO  BED  =  16.12  BTU/HR 

HEAT  FLUX  THRU  PLATE  =  106.48  BTU/HR 

h (PLATE-BED)  =                  3.17  BTU/HR-FT"2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  =  154.29  DEG  F 

AVERAGE  PLATE  SURFACE  TEMPERATURE  =  154.26  DEG  F 

AVERAGE  BED  TEMPERATURE  =  81.46  DEG  G 

ELECTRICAL  ENERGY  INTO  HEATER  =  160.10  BTU/HR 

LOSS  FROM  HTR  BACKING  TO  ATM  =  6.44  BTU/HR 

h(HTR  BACKING-ATM)  =             .35  BTU/HR-FT~2-DEG  F 

LOSS  FROM  HTR  PERIMETER  TO  BED  =  24.67  BTU/HR 

HEAT  FLUX  THRU  PLATE  =  128.98  BTU/HR 

h  (PLATE-BED)  =                  5.10  BTU/HR-FT -'-2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE  = 

AIR  OUTLET  TEMPERATURE  = 

AVERAGE  BED  TEMPERATURE  = 

q (AIR  OUT)  = 

q(LOSS  F/R  WALL)  = 

q (TOTAL  OUT  OF  BED)  = 

q (TOTAL  INTO  BED)  = 

SUPERFICIAL  VELOCITY  = 

SUPERFICIAL  MASS  VELOCITY  =         14: 

PARTICLE  REYNOLDS  NUMBER  = 


75 

\.bA 

■  DEG  F 

82 

\.9A 

•  DEG  F 

81 

.46 

>  DEG  F 

126. 

26 

BTU/HR 

0. 

00 

BTU/HR 

126. 

26 

BTU/HR 

276. 

25 

BTU/HR 

. 

hJlJ 

FT/SEC 

4  LBM/HR-FT 

2.37 
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5.   A  typical  data  run.   (Bed  thermocouple  probe  #57  reads 
significantly  higher  than  -other  bed  thermocouples. 
Reading  exceeds  highest  heater  thermocouple  reading . ) 


NOTE:   ALL  VALUES  WERE  RECALCULATED  AFTER  OMITTING 

THERMOCOUPLE  NR.  57  READING.  REVISED  RESULTS 
ARE  PRESENTED  IN  TABLES  AND  CURVES  PREVIOUSLY 
SHOWN . 


**EXPERIMENTAL  DATA** 
RUN  #  06B36M 
06/08/82  -  1630 


RUN  NUMBER  = 
BED  WIDTH  = 
STATIC  BED  HEIGHT  = 
AIR  FLOW  RATE  = 
AMBIENT  TEMPERATURE  = 
RIGHT  HEATER  VOLTAGE  = 
RIGHT  HEATER  CURRENT  = 
LEFT  HEATER  VOLTAGE  = 
LEFT  HEATER  CURRENT  = 
BED  EXPANSION  = 


06B36M 

6.25  INCHES 

12.5   INCHES 

12.48  CFM 

74.50  DEG  F 

60.0  VOLTS 

1. 12  AMPS 

60.0  VOLTS 

1.28  AMPS 

0.25  INCHES 


COMMENTS:  VIOLENT  FLUIDIZATION  AT  CARRYOVER 
LIMIT 
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♦^EXPERIMENTAL  DATA** 
RUN  #  06B36M 
06/08/82  -  1630 

THERMOCOUPLE  READINGS 


TC# 

OUTPUT 

TEMP 

(mV) 

(DEG  P) 

1 

1.930 

117.46 

2 

1.936 

117.68 

3 

1.951 

118.31 

4 

.980 

79.43 

5 

2.004 

120.53 

6 

2 .  004 

120.50 

7 

2.005 

120.57 

8 

2 .  000 

120.34 

9 

2.096 

124.34 

10 

2.089 

124.06 

11 

2.081 

123.72 

12 

2.088 

123.96 

13 

2.  180 

127.75 

14 

2.  185 

127.92 

15 

2.  164 

127.07 

16 

2.  166 

127. 11 

17 

2.209 

1 28 . 38 

18 

2.219 

129.27 

19 

2.214 

1 29 . 09 

20 

2.208 

128.96 

21 

1.936 

117.53 

22 

1.948 

118.03 

23 

1.947 

117.96 

24 

1.944 

117.82 

25 

2.005 

120.38 

26 

1.987 

119.65 

27! 

1.997 

1 20 . 04 

28 

2.017 

120.87 

29 

Z.  06 1 

122.74 

30 

2.053 

122.35 

31 

2.051 

122.30 

32 

2.056 

122.44 

2.071 

1 23 . 09 

34 

2.067 

122.38 

35 

2.075 

123.21 

36 

2.079 

123.38 

37 

2.  080 

123.42 

33 

2.  053 

122.50 

39 

2.055 

122.37 

40 

2.094 

124.25 

TC# 

OUTPUT 

TEMP 

(mV) 

(DEG  P) 

41 

1.773 

110.81 

42 

1.774 

110.85 

43 

1.726 

108.33 

44 

1.966 

118.32 

45 

2.  108 

124.75 

46 

2.007 

120.55 

47 

.936 

74.62 

48 

.935 

74.56 

49 

.935 

74.57 

50 

1.631 

104.85 

51 

1.592 

103.22 

52 

1.621 

1 04 . 49 

53 

1 .  608 

103. 89 

54 

1.609 

103.90 

55 

1.579 

102.38 

56 

1.629 

1 04 . 75 

57 

2.260 

131.38 

68 

1.582 

102.44 

69 

1.658 

105.70 

70 

1.942 

117.73 

71 

2.025 

121.28 

72 

1.565 

101.71 

73 

1.626 

104.32 

74 

1.313 

90 .  8 1 

75 

1.497 

98.78 

76 

1.768 

110.39 

77 

1.705 

107.71 

78 

1.682 

106.72 

79 

1 .  740 

109.20 
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* ♦EXPERIMENTAL  DATA** 

RUN#  06B36M 

06/08/S2  -  1630 


RIGHT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1    COL  2 


COL 


COL  4 


ROW    1  !  117.46  117.63  118.31  ****** 

ROW    2  !  120.53  120.50  120.57  120.34 

ROW    3  !  124.34  124.06  123.72  123.96 

ROW    4  !  127.75  127.92  127.07  127.11 

ROW    5  !  128.88  129.27  129.09  128.96 


AVE 
117.82 
120.48 
124.02 
127.46 
129.05 


NOTE:     TEMPERATURES    SHOWN    ARE    IN    DEGREES    FAHRENHEIT 
AND    ARE    ARRANGED    AS    SEEN    FROM    WITHIN    THE    BED 


RI  GHT-HAND     HERTER 
TEMP     vs     DEPTH 

f  1  0 

Z     9 

8 

I     7 

~     6 


en 

O 

CL 


on 

UJ 


a 


m 


CvJ 


TEMP     (DEG     F) 


NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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♦♦EXPERIMENTAL  DATA^ 
RUN#  06B36M 

06/08/32  -  1630 


LEFT-HAND  HEATER 
TEMPERATURE  PROFILE 


COL  1 


COL 


COL  3    COL  4 


ROW 

1 

ROW 

■■J 

ROW 

o 

ROW 

4 

ROW 

5 

17.53  113.03  117.96  117.32 

20.38  119.65  120.04  120.87 

22.74  122.35  122.30  122.44 

23.09  122.88  123.21  123.38 

23.42  122.50  122.37  124.25 


AVE 

117. 

34 

1 20 . 

24 

1  ^T1 

46 

1  *^T 

14 

1  ">~^ 

13 

NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED 


z  1  0 
Z     3 

8 

i  7 

~  6 

I:  5 

4 
3 
2 
1 

a 


o 

a. 


(X. 


LEFT-HRND  HERTER 
TEMP  vs  DEPTH 


♦ 
♦ 


<3 


IT) 


in 

CM 


TEMP  (DEG  F) 


in 


NOTE:  PLOTTED  TEMPERATURES  APE  AVERAGES  OF  EACH 
HORIZONTAL  ROW 
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♦^EXPERIMENTAL  DATA** 

RUN#  06B36M 

06/08/32  -  1630 


FLUIDIZED  BED 
TEMPERATURE  PROFILE 


104.75   103.89   104.85 


103.2: 


131.38   103.90   104.49 


NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT 
AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED 
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♦^CALCULATED  RESULTS** 
RUN  #  06B36M 
06/08/S2  -  1630 


RIGHT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  =  123.77  DEG  F 

AVERAGE  PLATE  SURFACE  TEMPERATURE  =      123.71  DEG  F 
AVERAGE  BED  TEMPERATURE  =  108.07  DEG  F 

ELECTRICAL  ENERGY  INTO  HEATER  =  229.29  BTU/HR 

LOSS  FROM  HTR  BACKING  TO  ATM  =  2.96  BTU/HR 

h(HTR  BACKING-ATM)  =  .19  BTU/HR-FT~ 2-DEG  F 

LOSS  FROM  HTR  PERIMETER  TO  BED  =  8.13  BTU/HR 

HEAT  FLUX  THRU  PLATE  =  218,20  BTU/HR 

h  (PLATE-BED)  =  40.18  BTU/HR-FT""2-DEG  F 


LEFT-HAND  HEATER 


AVERAGE  HEATER  TEMPERATURE  =  121.36  DEG  F 

AVERAGE  PLATE  SURFACE  TEMPERATURE  =      121.29  DEG  F 
AVERAGE  BED  TEMPERATURE  =  108.07  DEG  G 

ELECTRICAL  ENERGY  INTO  HEATER  -         262.05  BTU/HR 
LOSS  FROM  HTR  BACKING  TO  ATM  =  .00  BTU/HR 

h(HTR  BACKING-ATM)  =  .00  BTU/HR-FT^-DEG  F 

LOSS  FROM  HTR  PERIMETER  TO  BED  =  4.41  BTU/HR 

HEAT  FLUX  THRU  PLATE  =  257.64  BTU/HR 

h (PLATE-BED)  =  56.11   BTU/HR-FT"2-DEG  F 


ENERGY  BALANCE 


AVERAGE  AIR  INLET  TEMPERATURE  = 

AIR  OUTLET  TEMPERATURE  = 

AVERAGE  BED  TEMPERATURE  = 

q (AIR  OUT)  = 

q (LOSS  F/R  WALL)  = 

q (TOTAL  OUT  OF  BED)  = 

q (TOTAL  INTO  BED)  = 

SUPERFICIAL  VELOCITY  = 

SUPERFICIAL  MASS  VELOCITY  = 

PARTICLE  REYNOLDS  NUMBER  = 


74 

.51 

3  DEG  F 

102 

.31 

3  DEG  F 

108 

.  0" 

7  DEG  F 

353. 

81 

BTU/HR 

on 

30 

BTU/HR 

376. 

11 

BTU/HR 

488. 

33 

BTU/HR 

. 

30 

FT/SEC 

202.78  LBM/HR-FT~2 

4.05 
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APPENDIX  F . 
EQUIPMENT  LISTING 


DEVICE 

Axial  Flow 
Turbo  Compressor 

Rotometer 

30- Inch 

Water  Manometer 

Computer 


Data  Acquisition 
System 

Dual  Flexible  Disk 
Drive 

Plotter 


Printer 


Copper  Constantan 
Thermocouples 

Regulated  DC 
Power  Supply 
60V,  2  A  (Left) 

Regulated  DC 
Power  Supply 
60V,  2  A  (Right) 

Computating 
Bridge 

Calibration 
Bath 

Strip  Heaters 
250  Watts 


MANUFACTURER 
Spencer  Turbine 

Fisher  §  Porter 
Me  r  i  am 


Hewlett-Packard 
Company 

Hewlett -Packard 
Company 

Hewlett-Packard 
Company 

Hewlett-Packard 
Company 

Hewlett -Packard 
Company 


Lambda 

Hewlett-Packard 
Company 

Rosemount 

Rosemount 

Watlo 


MODEL  # 
5  50  SCFM 

34  SCFM 
33KA35 


SERIAL  # 


Omega  Engineering     10  Gauge 


A69513 


HP  85        2151A50182 

HP  3497A 

HP  8  2  901M 

HP  7225A 

HP  82905B    2143J12965 


LK345A 


6296A 


9  2  0A 


913A 


250  Watts    8044 


C1821 


1552A02618 


110 
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APPENDIX  G. 
EXPERIMENTAL  UNCERTAINTY  ANALYSIS 

In  order  to  arrive  at  an  estimate  of  the  difference 
between  measured  or  calculated  values  and  the  true  values 
of  various  parameters,  the  following  uncertainty  analysis 
was  accomplished: 
1.   TEMPERATURE  MEASUREMENTS 

A.  Thermocouple  Readings  -  Since  the  thermocouple  out- 
put voltage  was  measured  automatically  by  the  Data  Acquisition 
System,  human  error  in  reading  the  output  was  avoided;  however, 
the  values  arrived  at  by  the  System  were  still  subject  to 
random  error  as  well  as  inaccuracies  in  the  empirical  corre- 
lation developed  during  the  calibration  of  the  thermocouples. 
Because  of  these  effects,  it  is  believed  that  the  thermo- 
couple readings  are  accurate  to  within  ±  F. 

B.  Ambient  Temperature  Reading  -  This  measurement  was 
accomplished  with  a  mercury- in-glass  thermometer  having  mini- 
mum graduations  of  1°F  and  a  range  of  -20°F  to  +140°F.   It 

is  believed  that  the  thermometer  could  be  read  correctly  to 
within  ±0.5°F.   Allowing  for  the  error  within  the  device, 
the  uncertainty  of  the  thermometer  reading  is  considered 
to  be  ±1°F. 
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2.  PRESSURE  MEASUREMENTS 

A.  Bed  Pressure  Drop  -  The  minimum  manometer  scale  read- 
ing was  0.1  inch  with  a  total  range  of  30  inches.  Readings 
of  0.05  inches  were  possible.  For  static  measurements,  the 
pressure  readings  are  considered  to  be  accurate  within  0.21; 
however,  during  conditions  under  which  the  water  column  was 
fluctuating,  the  uncertainty  was  increased  to  as  much  as  2% 
of  the  full  scale  reading. 

B.  Venturi  Inlet  Static  Pressure  -  The  minimum  scale 
reading  was  0.1  inch  with  a  total  range  of  60  inches.   Readings 
of  0.05  inches  were  possible.   The  pressure  readings  are  con- 
sidered to  be  accurate  to  within  0.1%  of  full  scale. 

C.  Venturi  Differential  Pressure  -  The  minimum  scale 
reading  was  0.001  inches  over  a  range  of  2  inches.   The 
reading  was  considered  to  be  accurate  to  within  0.1%  of 
full  scale. 

D.  Atmospheric  Pressure  -  Minimum  scale  readings  of 
0.001  in.  Hg  over  a  range  of  approximately  32  inches  were 
observed.   As  a  consequence,  the  uncertainty  of  the  reading 
was  considered  to  be  0.1%  of  full  scale. 

3.  PI  STANCE/ SIZE  MEASUREMENTS 

A.  Bed  Width/Depth  -  Measurements  were  considered  . 
accurate  to  within  ±0.0625  in.  or  0.5%  of  the  apparatus 
depth. 
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B.  Bed  Probe  Depth  -  Due  to  the  bending  motion  imparted 
to  the  probes  by  the  particle  activity,  the  depth  measure- 
ments were  considered  accurate  to  within  the  0.25  in.  or 

1  .  S%    of  the  bed  depth. 

C.  Bed  Expansion  -  Due  to  the  rapidly  fluctuating  ex- 
panded bed  surface  level,  the  bed  depth  is  considered  to 
have  an  uncertainty  of  ±0.5  in.  or  5  %    of  bed  depth. 

D.  Venturi  Inlet/Throat  Diameter  -  These  measurements 
are  considered  accurate  to  within  ±0.001  in.  or  0.2%  of 
diameter . 

E.  Particle  Diameter  -  Determined  by  optical  microscope 
examination  conducted  by  Morgan,  [Ref.  2],  D   =  0.0122  (310 
microns)  ±  1% . 

F.  Plate/Wall  Thickness  -  Measurements  were  considered 
accurate  to  within  0.02  in.  or  5%  of  the  wall  thickness. 

4.   VOLTAGE/ CURRENT  MEASUREMENTS 

Minimum  ammeter  scale  graduations  were  0.1  Amp  with  a 
total  range  of  2  Amps.   Because  the  power  supplies  were  not 
calibrated  or  compared  with  a  known  standard  an  uncertainty 
of  ±10%  of  the  full  scale  deflection  is  assigned  to  the 
current  reading.   The  minimum  voltmeter  scale  reading  was 
1  volt  with  a  range  of  60  volts.   Once  again,  due  to  the 
lack  of  calibration,  an  uncertainty  of  5%  is  assigned  of 
full  scale  deflection. 
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5.  THERMAL  CONDUCTIVITIES 

Due  to  the  variation  in  conductivities  found  in  the 
literature,  the  following  uncertainties  are  assigned: 

A.  k      232  Btu/Hr-Ft-°F  ±0.51 

B.  k  ,  =  0.112  Btu/Hr-Ft-°F  ±10% 

mb 

C.  k.   =  0.025  Btu/Hr-Ft-°F  ±10% 

l 

D.  k£   =  0.112  But/Hr-Ft-°F  ±10% 

f  r 

6.  DENSITY  MEASUREMENTS 

Particle  Density  -  Established  by  Morgan  to  be  154.6 
Lbm/Ft2  ±3°6. 
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APPENDIX  H.  ROTOMETER  CALIBRATION  COMPUTER  PROGRAM  LISTING 


10  PRINTER  IS  701 
20  S-.S34 

30  CLEAR  <?  DISP  " 
40  DISP  "  " 
50  DISP 
60  PAUSE 

70  clear 

80  INPUT 
*0  CLEAR 
100  INPUT 


1  00 


ROTOMETER  CALIBRATION* 


THIS  PROGRAM  CONVERTS  VENTURI  PRESSURE  DATA  INTO  ROTOMETER  PLOW  RATES" 

*  DISP  "ENTER  Patm(IN  Hg>" 
P 

&    DISP    "ENTER    Tm(DEG    F>  " 
T 

Y" 


110  CLEAR  &  DISP  "ENTER 

120  INPUT  Y 

130     !     CALCULATE    Patin(PSI) 

140  Pl-Ptl4.7Z29.92 

150  !  CALCULATE  Tin(DEG  R) 

160  T1«T*459.69 

170  CLEAR  e  DISP  "PRESS  CONTINUE  TO  START  PROCESSING* 

180  PAUSE 

190  CLEAR 

200  INPUT 

210  CLEAR 

220  INPUT 

230  CLEAR 

240  INPUT 


&  DISP 

F 

e  disp 

12 

&  DISP 

13 

250  !   CALCULATE  FL0WRATE5 
260  !  CALCULATE  Pig (PS I) 
270  P2«I2*14. 7/33. 91/12 
280  !  CALCULATE  Pi(PSI> 
290  P3-P1-P2 

300  !  CALCULATE  DELTA  P(PSI) 
310  PA«I3*62. 4*. 834/12/144 


ENTER  FLOWMETER  READING" 


"ENTER  Pig (IN  K20) 


"ENTER  DELTA  P(IN  RO) 


!  CALCULATE  X 

X-P4/P3 

!  CALCULATE  RHOi 

D-P3* 144/53. 34 /Tl 

!  CALCULATE  MASS  FLOW (LBM/SEC) 

DISP  "ENTER  ASSUMED  REYNOLDS  NUMBER" 


320 

330 

340 

350 

360 

370  CLEAR 

380  INPUT 

390  CLEAR 

UMBER" 

400  INPUT 

E-1.0 

K«C*E 

M-K*3. 56327/ 144*Y*SQR (2*32. 174*D*P4*144) 

!  CALCULATE  VOLUME  FLOW  RATE (FT3/SEC) 


DISP  "ENTER  DISCHARGE  COEFFICIENT  ASSOCIATED  WITH  ASSUMED  REYNOLDS  N 


C 
52796 


410 

420 

430 

440 

450  Vl-M/D 

460  !  CALCULATE  VOLUME  FLOW  RATE (FT3/MIN) 

470  V2-V1*60 

480  !  CALCULATE  REYNOLDS  NUMBER 

490  R2-Vl/14.25309*144*D»4.26/12/.000000379/32. 174 

500  CLEAR  &  DISP  "ASSUMED  REYNOLDS  NUMBER  -",R 

510  ! 

520  DISP  "CALCULATED  REYNOLDS  NUMBER  «",R2 

530  ! 

540  DISP  "DO  YOU  WISH  TO  REPEAT  CALCULATION  USING  A  DIFFERENT  RE ASSUMED  REYNOLDS 

NUMBER?  (Y/N)" 
550  INPUT  Wl« 
560  IF  W1*«"Y"  THEN  360 
570  PRINT  CHR*(12) 
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530  PRINT  USING  590 

590  IMAGE  2/ . 2SX, " **ROTOMETER  CALIBRATION**" 

oOO  PRINT  USING  610 

610  IMAGE  23X."ASME  HERSCHEL-TYPE  VENTURI" 

620  PRINT  USING  630 

=30  IMAGE  33X, "4.26/2. 13  IN  DIA"620  PRINT  USING  630  ;  F 

640  PRINT  USING  650  :  F 

630  IMAGE  / .32X.DD. X, "PERCENT  READING" 

660  PRINT  USING  670 

670  IMAGE  32X," " 

680  PRINT  USING  690 

690  IMAGE  2/. 31X, " * *EXPERIMENTAL  DATA**" 

700  PRINT  USING  710  ;  P 

710  IMAGE  2/, 15X. "ATMOSPHERIC  PRESSURE  ■" . 1 4X , DD . DDDDD. X . " IN  Hg " 

720  PRINT  USING  730  :  T 

730  IMAGE  15X "VENTURI  INLET  TEMPERATURE  ■" , 13X . DD. D . X . "DEG  F" 

740  PRINT  USING  750  ;  12 

750  IMAGE  15X,"VENTURI  INLET  STATIC  GAGE  PRESSURE  »" , 2X . DD . DD. X . " IN  H20" 

760  PRINT  USING  770  ;  13 

770  IMAGE  15X."VENTURI  DIFFERENTIAL  PRESSURE  -" , 3X . D. DDD. X . " IN  RED  OIL" 

780  PRINT  USING  790  :  S 

790  IMAGE  15X, "SPECIFIC  GRAVITY  OF  RED  OIL  «" , 1 6X , D. DDD 

800  PRINT  USING  810 

810  IMAGE  2/.31X. "**CALCULATED  VALUES**" 

820  PRINT  USING  830  ;  PI 

830  IMAGE  2/. 15X. "ATMOSPHERIC  PRESSURE  «" . 17X , DD. DDDD, X , "PSI " 

840  PRINT  USING  850  ;  P2 

850  IMAGE  13X."VENTURI  INLET  GAGE  PRESSURE  «". 10X, D. DDDDD. X. "PSI " 

860  PRINT  USING  870  ;  P3 

870  IMAGE  15X, "VENTURI  INLET  ABSOLUTE  PRESSURE  «", 5X . DD. DDDDD, X , "PSI " 

880  PRINT  USING  890  :  P4 

8<?0  IMAGE  15X.  "VENTURI  DIFFERENTIAL  PRESSURE  «",  SX  ,  D.  DDDDD,  X  ,  "PSI  " 

900  PRINT  USING  910  ;  D 

910  IMAGE  15X. "VENTURI  INLET  AIR  DENSITY  -", SX . D. DDDDD, X , "LBM/FT3" 

920  PRINT  USING  930  ;  X 

930  IMAGE  15X. "VENTURI  PRESSURE  DROP  RATIO  <X)  -", 10X, D. DDDDD 

940  PRINT  USING  950  ;  Y 

950  IMAGE  15X. "EXPANSION  FACTOR  (Y)  =".24X,D.DD 

960  PRINT  USING  970  ;  R 

970  IMAGE  15X. "ASSUMED  REYNOLDS  NUMBER  -".21X.DDDD 

980  PRINT  USING  990  ;  C 

990  IMAGE  15X, "DISCHARGE  COEFFICIENT  (C)  «".1SX,D.DDD 

1000  PRINT  USING  1010  ;  E 

1010  IMAGE  15X. "VELOCITY  OF  APPROACH  FACTOR  (E)  -",1  OX , D. DDDDD 

1020  PRINT  USING  1030  :  K 

1030  IMAGE  15X,"FL0W  COEFFICIENT  (K)  -", 21X, D. DDDDD 

1040  PRINT  USING  1030  ;  M 

1050  IMAGE  13X."AIR  MASS  FLOW  RATE  «",  15X,  D.  DDDDD,  X,  "LBM/SEC 

1060  PRINT  USING  1070  ;  V2 

1070  IMAGE  15X,"AIR  VOLUMETRIC  FLOW  RATE  -" , 14X, DD. DDD. X, "CFM" 

1080  PRINT  USING  1090  ;  R2 

1090  IMAGE  1 5 X, "VENTURI  INLET  REYNOLDS  NUMBER  -"13X.DDDD 

1100  CLEAR  «■  DISP  "DO  YOU  WISH  TO  PROCESS  ANOTHER  SET  OF  READINGS?  (Y/N) 

1110  INPUT  W2* 

1120  IF  W2*«HY"  THEN  170 

1130  END 
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APPENDIX  I. 
THERMOCOUPLE  CALIBRATION  COMPUTER  PROGRAM  LISTING 


20 
30 
40 


50 


38 


90 

1  0  0 


1 1 0 
12Q 

130 
140 

150 

1  c"0 
170 
180 

190 

2  8  8 

210 


220 
2  3  0 
2  ^  0 
2  5  0 

260 
27  0 

2  S  0 
2 '?  0 

3  0  0 
310 
3  2  0 


146 
:'50 


DI3P  "  DA'.vCAL" 

n  J  3  p  "  " 

DISP  "THIS  CALIBRATION  PROGR 
AM  IS  FOR  USE  WITH  ThE  HP-34 
97  DATA  ACQUIS  I TI OH  SYSTEM" 
DISP  "  " 

DISP  "UP  TO  50  THERMOCOUPLES 
MAY  BE  CALIBRATED  SIMULTAHI 
OUSLV" 
DISP  "  " 

DISP  "UP  TO  20  DIFFERENT  BAT 
H  TEMPERATURE  SETTINGS  MAY  B 
E  USED" 
DISP  "  " 

DISP  "THERMOCOUPLE  NUMBERS  S 
HOULD  AGREE  WITH  CHANNEL  HUM 
BERS" 
PAUSE 
DISP  u  " 
OPTION  BASE  1 
MASS  STORAGE  IS  "  =  D701  ,; 
SHORT  R<20,  50>  ,  M-:20,  2> 
DIM  T<2Q,2> 


DIS 

DIS 

AR 

I  HP 

D I  S 
D  I  S 

TUR 


p  ..  .. 

P  "HOW  MANY  THERMOCOUPLES 

E  BEING  CALIBRATED?" 

UT  K 
pt  ■■  it 

P  "HOW  MANY  BATH  THEMPERA 
E  SETTINGS  ARE  BEING  USED 


INPUT  H 

H*="AF" 

DISP  "  " 

DISP  "WHAT  IS  FIRST  CHANNEL 

NUMBER?" 

DISP  "  (TWO  DIGITS)" 

INPUT  Y$ 

Hf C33=Y$ 

C^VALCy*  > 

HfC53="AL" 

DISP  "  " 

DISP  "WHAT  IS  LAST  CHANNEL  N 

UMBER?" 

DISP  "(TWO  DIGITS;-" 

INPUT  HSC73 

H$C93="yTlVR5" 

C*  =  H$ 
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370 
336 

399 

4  0  0 
410 
420 
430 

4  40 
4  50 

460 
470 


4  90 

5  q  m 

5  i  0 

520 
530 
54  0 
5  5  0 
560 
5  7  0 
5  8  0 

5  9  0. 

6  0  0 

610 

62  0 
6  3  0 

64  0 
6  5  0 

6  6  0 
6  7  0 

6  8  0 
69  0 


DISP  "  " 

DISP  "WHAT  DO  YOU  WISH  TO  CA 

LL  YOUR  DATA  FILE'?" 

DISP  "  <SIX  CHARACTERS  MflX>" 

INPUT  D* 

FOR  1=1  TO  N 

DISP  "  " 

DISP  "    TEMPERATURE  SETTING 

#"  }  I 
DISP  "  " 

DISP  "WHAT  IS  THE  PLATINUM  B 
ATH  SENSOR  RESISTANCE?'1 

INPUT  T<I, 1> 

DISP  "  " 

DISP  "WHAT  IS  THE  CORRESPOND 

ING  BATH  TEMPERATURE  vC>?" 

INPUT  T<I,2> 

DISP  "  " 

DISP  "TO  START  SAMPLING  PRES 

S  CONTINUE" 

PAUSE 

clear  709 

OUTPUT  ?fr3    ;  C$ 

FOR  J=l  TO  K 

OUTPUT  709  ; "AS" 

ENTER  709  ;  B 

R  <  I ,  J  >  =  B  *  1  O  O  0 

NEXT  J 

BEEP 

MEXT  I 

DiSP  "  " 

CLEAR  »*  CISP  "       SAMPLING 

COMPLETE" 
DISP  "  " 

DISP  "RESULTS  WILL  BE  PRINTE 
D  MOMENTARILY" 
DISP  "  " 
DISP  "THERMOCOUPLE  OUTPUT  IS 

IN  ft  I  LI VOLTS" 
DISP  "  " 
DISP  "BATH  TEMP. IS  IN  CELSIU 


700  DISP  "  " 

710  PRINT  "DATA  IS  STORED  IN  FIL 

E  NAMED  " 
720  PRINT  USING  730  ;  CHR*(34>,D 

* ,  CHRi-<'34> 
""30  IMAGE  11>wA.6A,A 

46  CREATE  Q$, 1 , 165Q9 

50  ASSIGNS  l  TO  Dt 

60.  E  =  C-i 

70  FOR  1=1  TO  K 
780  FOR  .  J=l  TO  N 
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796  M(  J,  1  )=R< J,  I  :■ 

See  M<J,2>=TCJ,2) 

810  PRINT*  1  .;  M<  J,  1  )  >  MC  J,  2> 

320  NEXT  J 

330  F  =  E+I 

346  PRINT  USING  350 

35-0  I M  R  G  E  2  • '  >  "  X  *  *  *  X  X  X  *  *  *  *  *  4:  *  *  *  *  X 
J**:*::*:***:*:.***" 

360  PKlNT  USING  370 

370  IMAGE  "  *:******:«:******#**#** 

■T-  X  -T-  X  -T-  ^t-  +  *•  X        ^  /  ^  A 

830  PRINT  USING  850  ;  F 

890  IMAGE  /,6X, "THERMOCOUPLE  #", 

X  .>  2  D 

900  PRINT  USING  910 

910  1  M  A  G  E  6X  *  "  *  *  4: 4: 4  *  X  X  X  X  X  X  X  X  X  X  X  " 

,2/,X 

920  PRINT  USING  930 

?3^     I  MAGE  "  TEMP  "  ,  4X  ,  "  PLAT  "  ,  t~X  ,  "  T 

C",5X, "BATH" 

940  PRINT  USING  950 

950  IMAGE  X, "SET", 3X, "RESIST" , 3X 

,  "  OUTPUT  U,3X,  "TEMP" 
368  PRINT  USING  970 
970  I  MAGE  " XXXX " >  3X ,  " XXXXXX " , 3X , . 

"XXX XXX" ; 3X,  "XXXX" , / 
930  FOR  L=l  TO  N 
SSib    PRINT  USING  1000  ;  LiKLiDi 

m  ci_ ,  i  > ..  n  <:  l  >  2 :- 

1  0  O  O  I  M  AGE  X  ,  D  0  ,  4  X  ..  D  D '.  0  D  D  D  >  2  X  ,  D  D 

. DC  J, 3X, 00 . 00 
1010  NEXT  L 
1620  CLEAR  15  OISP  "      THERMOCO 

UPLE  S"  ,F 
1636  OISP  "READY  FOR  LEAST  SQUAR- 
ES CURVE  FITTING" 
1 0  4  6  OISP  "  ■ 
1656  OISP  "ENTER  O  FOR  A  POWER  F 

UNCTION  -:V  =  AX-B>" 
1  6  6  6  OISP  "  " 
1676  OISP  "ENTER  1  FOR.  AN  EXPONE 

NTIAL  FUNCTION  <Y  =  Ae^BX:>" 
1  6  8  6  0  ISP  u  " 

1630  OISP  "ENTER  2  THRU  16  FOR  A 
POLYNOMIAL  OF  N  TERMS  ''.DEb 

PEE  =  H-D" 
1106  DISP  "  >'.  Y  =  C1  i-C2XJ'C3X'-2  +  C4X'v3 

+  .  .  .  .  C «;  N  >  X  y-  <  N  -  1  >  >  " 
1110  0  I  M  ' r '  >:.  2  6  '■'  .-  X  <  2  6  >  •  A  <  1  6  .•  1  6 ')  >  E  < 

16.  1  O  >  ..  C  <  1  O  )  >  0  <  1  0  ) 
1120  MAT  Y  =  2ER>:N> 
1136  MAT  X=ZERCN> 
1146  FOR  J=l  TO  H 
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1150 

1  160 

1170 

1186 

1  130 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1230 

12  30 

1  3  0  U 

1310 

132  0 

1330 

1340 

1350 

1  360 

13  70 

1  380 

1 3  y  0 

14  00 

1410 

14  20 

14  30 

1  4  40 

1450 

!  460 

14  70 

1480 

14  30 

1  5  0  0 

1510 

1520 

1530 

1  540 

1550 

156  0 

1  ^  ~» a 

1530 

1530 

1  6  0  0 

}<<J}=M<rj,  1) 
Y< J)=M< J,2> 

next  J 

INPUT  Q 

REM  "CALCULATE  LOGARITHMS  G 

F  X-  HMD  Y-VALUES  IF  NECESS 

fl  R  Y  " 

IF  Q>=2  THEN  1280 

FOR  J=l  TO  N 

LET  Y<  J>=LuG<M(  J.  2  :•  > 

NEXT  J 

IF  Q=l  THEN  1280 

FOR  J=l  TO  N 

LET  X(J>=LOG(f'HJ,  1  >) 

NEXT  J 

REM  "CALCULATE  ELEMENTS  OF 

A-MATRIX  AND  D-VECTOR" 

LET  N1=Q 

IF  Nl;=2  THEN  1320 

LET  Nl=2 

MAT  A=ZER<N1,N1> 

MAT  D=ZER<N1) 

FOR  J»l  TO  HI 

FOR  K=l  TO  HI 

IF  J+K>2  THEN  1390 

LET  A<J,K)=H 

GHTO  1420 

FOR  L=l  TO  H 

LET    fl<  J,  K>=A<  J,  K>+X<L>M  J+fc 

-2  > 

NEXT  L 

NEXT  K 

FOR  L=l  TO  N 

IF  J> 1  THEN  1470 

LET  D<J>=D<J>+Y<L) 

GOTO  1480 

LET  DCJ>=D<J>+Y<L>tXCL>A<J- 

n 

NEXT  L 

NEXT  J 

REM  "SOLVE  SIMULTANEOUS   LI 

NEAR  EQUATIONS" 

MAT  B=INV<fl> 

MAT  C  =  B*G 

REM  "PRINT  EQUATION  FOR  CUP 

','E  FIT" 

PRINT  USING  1550 

IMAGE  /," 


IF  Q> 1  THEN  1670 

LET  Cl=EXP-::C':i  )  > 

IF  Q=l  THEN  1630 

PRINT  USING  1600 

IMAGE  5  X ,  " * % # p 0 W E  R  F  U H C T I  0 N 

*.%*",sf  X 
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1  6  1  9 
1628 
1638 
1640 

1659 

1660 
1670. 

1639 
1630 

1780 
1710 

1720 
1730 
1740 
1750 

1760 
17  70 

17S9 


17 

30 

IS 

:00 

18 

18 

IS 

!  2  0 

IS 

3  8 

IS 

:40 

1  s 

:=;0 

IS 

60 

IS 

70 

IS 

S  0 

IS 

30 

1  3 

0  0 

13 

10 

13 

20 

13 

30 

13 

4  0 

«  r 

50 

13 

60 

13 

70 

13 

S  8 

13 

30 

28 

80 

28 

10 

20, 

20 

28 

30 

28' 

4  0 

PRINT  "Y="  iClJ  "  TX-",  C<2> 

GOTO  1360 

PRINT  USING  1640 

IMAGE  3X, "***EXPONENTIAL  FU 

NOTION***" ,/* X 

PRINT  "Y=" ;C1; "*EXP< ";C<2>; 

"*X>" 

GOTO  1360 

IF  C-2>>=0  THEN  174  8 

PRINT  USING  1630 

IMAGE  3X,  ,,***POLYNOMIAL  FUN 

CTION***" 

PRINT  USING  1710  ;  Q-l 

IMAGE  18X, " (ORDER  =",0,")", 

/  >  X 

PRINT  "Y=" ;C(l)iC<2);MX"; 

GOTO  1738 

PRINT  USING  1758 

IMAGE  3X, "***POLYNOMIfiL  FUN 

CTION***" 

PRINT  USING  1778  ;  Q-l 

IMAGE  18X,"  <  ORDER  =  U,D,")% 

/iX    . 

PR  I  NT  " Y= " i CC 1 >  i  " +" ; C<  2> ; " * 

Xu; 

IF  Q=2  THEN  1368 

FOR  J=3  TO  Q 

IF  C(J)>=0  THEN  1340 

PRINT  C<J>.;  n*X~"  .;  J-l; 

GOTO  1350 

PR  I  NT  "  +  "  .;  C  C  J  >  ;  " *XA  M  ;  J-l  J 

NEXT  J 

REM  "PRINT  INPUT  VALUES  OF 

X  AND  Y  AND  CALCULATED  VRLU 

ES  OFY" 

IF  Q>=2  THEN  1358 

FOR  J=l  TO  N 

LET  Y<J>=EXP<Y<J)> 

NEXT  J 

IF  Q=l  THEN  1350 

FOR  j=l  TO  N 

LET  X<J>=EXP<XCJ>) 

NEXT  J 

PRINT 

PRINT  .USING  1370 

IMAGE  /, 4X, "X" 

- 5X , " Y  < C A L C )  " 

PRINT  USING  1338 

IMAGE  3X,  "  ***"  .  5! 

, 5X> "********" 

LET  3=0 

FOR  J=l  TO  N 

IF  Q>=2  THEN  2038 

IF  G!=l  THEN  2-868 

LET  Yl=Cl*XCJ>*CC2> 


ACT  ':■'" 
":*:*..*:*:*.*:*" 
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2  6  5  6 
2  6  €<  @ 
2  G70 
2039 

2090 

2100 
2 1 1  0 
2120 
2130 

214  0 

2150 
2160 
2170 

2150 


2 

159 

2 

200 

^.' 

210 

2 

220 

o 

230 

2 

240 

^ 

250 

2 

2t"0 

- 

2  SO 

-. 

230 

- 

3  0  0 

2310 


d  s  2  O 


Su 


2340 

235U 


£  --■  €•  0 


i .-.  7  0 
2380 


GOTO  2120 

LET  Y 1 =C 1 *EXP <C<£> *X C J > > 

GOTO  2120 

LET  Y1=C(1) 

FOR  K=2  TO  Q 

LET  Y 1 =Y 1 +C< K ) *X < J ) '• < K-l > 

NEXT  K 

LET  S=S+CY< J>-Y1 >-2 

PRINT  USING  2140  ;  X C J > , Y  < J 

>  ,  Y  1 

I  MAGE  X  ■  DD  .  ODD ,  5X  ,  DD  .  DO  ,  6X  > 

DO  CiDOn 

NEXT  J 

PRINT 

PRINT  "SUM  OF  SQUARE  ERRORS 

=  U;S 

CLEAR  (g  OISP  "DO  YOU  DESIRE 

A  PLOT?" 
DISP  "CENTER  Y  OR  N)° 
INPUT  PS- 

IF  P*="N"  THEN  3690 
GOSUB  2250 
PAUSE 
GOTO  3690 

D  I M  Y  K4)ifll(5);«U5) 
VK1)/Y1(3)»IHF  6  Y1C2;,Y1< 
4>=-lNF  @  fll(l)..fll(2),fll(3) 
,  RK4).;  Al  <5>=6 
FOR  J=l  TO  N 

IF  yi<2XX<J)  THEN  Y1(2)=X( 
J> 
IF  yi(l»X(J)  THEN  yi(l)=X( 

•J> 

IF  V1(4XY(J)  THEN  Vl(4)  =  r( 

J> 

IF  yi(3»V(J)  THEN  yi(3)=Y( 

J> 
Al < 1 >=A1 CI )+X< J>  @  Rl<2>=fil 

C  2  >  +  X  C  J  )  t  X  en  i*  R 1  ■: '3  '■<  =  A 1  C  3  ) 

+  YCJ> 

Hl(4)=fll(4)  +  v(J)l:Y(J)  g  fllC 

5)=A1  <5>+X<J>*Y<J>. 

NEXT  J 

m<t)»fll(i)/N  e  Ml  C2>  =  CA1<2 

>  -A  1  <  1  >  *  A  1  C  1  >  /-  N  )  •  C  N- 1  >  (5  M 1 

<3>=A1  C3>.-'N 

Ml  <4>  =  CA1  C4>-A1  C3>*A1  <3>-'N> 

/CN-1>  \'i    Ml  (5>  =  CA1  C5)-A1C1  > 

*  A  1  C  3  )  .'N  >  /  c  A  1  C  2  >  -  A 1  C  1  ,'  *  A  1  <  1 

R5  =  M1  C3>-H1  <  1  >*M1  C5>  6  D5  =  0 
CLEAR  e  D5=l  g  DISP  "AUTO  X 
-SCALING? = Y/N" 
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2390 
2460 
2410 

2420 

2430 
2440 

2450 
2460 

2470 


:49Q 


2510 
2520 

2530 

2540 

2550 


o 

60 

2 

c 

70 

2 

30 

2 

90 

2 

£ 

00 

£ 

6 

10 

2 

6 

20 

.-*• 

6 

30 

2 

6 

40 

o 

6 

50 

2 

K 

60 

2 

6 

70 

^ 

6 

80 

INPUT  Rf 

ON  FNR  GOTO  23S0 , 24 1 0 , 2420 

xi=yki>  »2  ;:2=yi(2)  e  goto 

2470 

CLEAR  6  DISP  "ENTER  MINIMUM 

VALUE  FOR  X-AXiS" 
INPUT  XI 
CLEAR  8  DISP  "ENTER  MAXIMUM 

VALUE  FOR'  X-AXIS"  . 
INPUT  X2 

IF  X1>=X2  THEN  BEEP  fi  GOTO 
24  20 

CLEAR  15  DISP  "DO  YOU  WISH  V 
ERTICAL  OR  HORIZONTAL  LABEL 
S  ON  THE  X-AXIS?  =  V.-'H" 
INPUT  Rf 

IF  NOT  LEN<R*>  THEN'  BEEP  @ 
G  0  T  0  247  0 
S 3 -  1  5  IF  u -  C f  < '  R t L  1  ,  1  3  :■  =  "  V  " 

THEN  S3=0  G  GOTO  2520 
IF  UPC* '.R-flli,  1II>#"H"  THEN  B 
EEP  5  GOTO  2470 
CLEAR  5  DISP  "DESIRED  NUMBE 
R  OF  INTERVALS  ON  X-AXIS?  < 
<=16>" 
INPUT  LI 

IF'LKl  OR  Ll>16  OR  L1#IHT< 
LI)  THEN  BEEP  6  GOTO  2520 
CLEAR  S  DISP  "DESIRED  NUMBE 
R  OF  SUBDIVISIONS  BETWEEN  L 
ABELED  INTERVALS?" 
INPUT  L3 

SS  =  0  S  IF  L3  =  0  THEN  S8=l  @ 
L3=l 

IF  L3<1  OR  L3>L1  THEN  BEEP 
S  GOTO  2550 
CLEAR  i?  DISP  "DO  YOU  DESIRE 

AUTOMATIC  SCALING  OF  THE  Y 
-AXIS?: Y/H" 
INPUT  R* 

ON  FNR  GOTO  2530 , 2670 ,  2620 
CLEAR  S  DISP  "ENTER  MINIMUM 
DESIRED  VALUE  FOR  Y-AXIS" 
TNPUT  Yl 
CLEAR  @  DISP  "ENTER  MAXIMUM 

DESIRED  VALUE  FOR  Y-AXIS" 
INPUT  Y2 

IF  Y1>Y2  THEN  BEEP  S  GOTO  2 
620  ELSE  26S0 
Y  -1  =  Y 1  <  3  >  \'i    Y  2  =  Y  1  <  4  > 
CLEAR  S  DISP, "ENTER  DESIRED 

N Utt BER *  OF  INTER V A L S  0 N  THE 

Y-AXIS  <:<=12>" 
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2590 

2700 

2710 


272A 
2730 

2740 

2759 

2760 

2770 

2730 
2750 
2300 
2S1Q 
2320 
2330 

234  0 
2350 
2  3  6  0 

2  3  7  0 
2330 
2  3  5  0 
2900 
2910 
2920 
2530 
2540 

2550 


2560 
2970 
2930 
2550 

3000 
3  0  1  0 

3  0  2  0 
3030 

INPUT  L2 

IF  L2<1  OR  L2M2  OR  L2#IHT< 
L2>  THEN  E:EEP  £  GOTO  2680 
CLERR  G  DISP  "ENTER  DESIRED 
NUMBER  OF  SUBDIVISIONS  BET 
WEEN  LABELED  INTERVALS"   . 
INPUT  L4 

S7  =  0  (2  IF  L4  =  0  THEN  S7=l  @ 
L4  =  l 

IF  L4<1  OR  L4>12  THEN  BEEP 
ik    GOTO  2710 

Z1=X2-X1  i?  Z3=IMT<20@/L1)*L 
1  £  Z2  =  Y2-Y1  ii  Z4  =  INT(144/L 
2>*L2 

D1=Z1-'Z3  £  D2=Z2'Z4  G  Z5  =  X2 
+  <207-Z3):i:Ol  5  X0=X1-48*D1 
5  Y0=Y 1-40*02 

GCL-ERR  6  SCRLE  XO,Z5>Y0,Y2+ 
U51-Z4>*D2 

XRX IS  Y 1 , Z 1 /L 1 , X 1 , X2  +  D 1 
YflX I S  X 1 , Z2/L2 , Y 1 , Y2  +  D2 
IF  S3  THEN  2569 
W 1 =L GT  <  RBS ( Z 1 /L 1 *L3  > ) 
W=LGT<RBS<X1+<X1=0)>)+1 
IF  W>5-(.SGNCX1)=-1>  OR  Wl<- 
3  THEN  2520 
IF  39  THEN  3520 
LDIR  50 

FOR  J  =  X1  TO  X2+D1  STEP  Zl/L 
l*L3 

MOVE  J-4i:Dl,Y0 
GOSUB  3320  • 
LRBEL  V*C1,V3 
NEXT  J 
GOTO  2560 
MOVE  X0,Y0 
LDIR  O 
I1=X1  15  GOSUB  3450  @  I  1  =  Z 1  -^ 

li  e  z$=v$  e  gosub  3450 

LRBEL  " XM I N= " &Z3X1 ■ J &3& " : T I 
CS=M«<V$C1,  103 

IF  S7  THEN  3140 

Wl=LGT<flBSv Z2/L2*L4>> 

W=LGT  (.  RBS  (  Y  1  +  <  Y 1  =0  >  )  )  + 1 

IF  W>5-<SGN<Y1 >=-l>  OR  Wl<- 

3  THEN  3070 

LDIR  O 

FOR  J=Y1  TO  Y2+D2  STEP  Z2/L 

2*L4 

MOVE  X8jJ-4*02 

GOSUB  3320 
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364'? 
3  8  5  C 
3u6  C 
387  5 
3  8  S  8 
3  8  9  9 
3180 
31  lO 
7128 
3138 
3148 
3158 
3160 


3170 

3  ISO 

3  1  9  0 
3  2  8  0 
3210 
3  2  2  0 
3230 
3240 
3258 
7  ?  t-,  0 
327  0 
3  2  S  8 
3  2  9  0 
3  3  8  0 
3318 
3320 

3348 
3350 
3368 


3370 

3330 
3398 
3400 
341  0 

3  4  2  8 
3430 


4  4  0 

450 
:4£fl 


LABEL  VSC1,V3 

NEXT  J 

GOTO  3140 

LDIR  30 

I  1 =  Y 1  5  GOSUB  3450 

MOVE  X0+12*B1,Y1 

LABEL  "  YM I  N=  "  StUfC  1  .183 

MOVE  X0+24*D1,Y1 

I1=Z2/L2  G  GOSUB  3458 

LABEL  "TICS»B&V*C1. 183 

PEHUP  @  LDIR  0 

MOVE  X<1>,YU) 

FOR  -J=l  TO  H  5  MOVE  X<J>-2* 

D1,Y<J>-4*D2  a  LABEL  "*M  8 

NEXT  J  @  BEEP 

FOR  L=0  TO  188 

X4=Y1 < 1 >  +L*<  <  Y 1 < 2 >  -Y  1  <  1  )  >  / 1 

80  > 

IF  Q>=2  THEN  3250 

IF  Q=l  THEN  3238 

LET  Y4=C1*X4^C<2> 

GOTO  3290 

LET  Y4=C1*EXP<C<2>*X4> 

G  0  T  0  329  0 

LET  Y4=C<1> 

FOR  K=2  TO  Q 

LET  Y4  =  Y4  +  COO*X4'-<K-1> 

NEXT  K 

PLOT  X4,Y4 

NEXT  L 

BEEP  i2  GOTO  3390 

||  .j-  _  ii  n       £       y  _  j 

V#=VflL*<X> 

r>    POS<V*i"E")    THEN    3408 

G9  =  LGT<ABS(X+<X  =  0>>  > 

IF    LEN-:V$>>5    AND    RBS<G9?'>4- 

<SGN<X>»-n    THEN    V  =  5    G    V$  =  " 

»    e    RETURN 

IF  LEN<V*><5  THEN  V=LEN<VS> 
@  RETURN 

V  *  C  1  ,  5  3  =  V  R  L  $  <  X  )  e  V  =  5 
GRAPH  6  RETURN 
E0=POS<V* , "E"X 

IF  VtCl, 13="-u  THEN  V1C33=V 
fCEGG  G  GOTO  3430 
W*C23=V*CE03 

V  =  LEH<V$>  (5  IF  V>5  THEN  PR  I 

NT  USING  3560  ;  J  6  V=5  t*  V 
■t  —  " — " 

RETURN 
U  -$  =  "  " 

0*=vnC*<m 
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3470 
3480 
3490 

3509 
3519 
3520 
3530 

3540 
3559 


3560 

3570 
3  5  8  0 

3590 

3S&& 
3610 
3629 
363% 
364  0 

3650 


IF    P0S<: ','$,  "E"  >     THEN    3500 
KtSLl,  101=VflL$<  ID 
-RETURN 

V*C6, 103=V*CPOS<V*, "Eu>3 
RETURN 

LDIR  0  i5  L9=-INF 
FOR  J  =  X1  TO  X2  +  D1  STEP  Z 1  •  L 
1*L3 

G  G  S  U  8  332  O 

IF    L?>J-<V*4  +  6V*01    OP    L9>Z5 

+<1-V*3>*D1  THEN  PRINT  US  IN 

G  3560  ik    GOTO  3610 

IMAGE  "LABEL  DELETED  AT  ",? 

D.4D 

M  0  V  E  J  +  C  2  -  V  *  4  >  #  D  1  ,  Y  1  - 1 2  *  D  2 

L3=J+<V*4+2)*01 

IF  L9>Z5  THEN  MOVE  Z5+<2-V* 

8>*D1 ,Y1-12*D2 

LABEL  V*C1;V3 

NEXT  J 

GOTO  2968 

DEF  FNR 

IF  NOT  LEN<R$>  THEN  11=1  @ 

GOTO  3660 

I  1 =POS  < " VN" , UPC* < R*C1 .13))  + 


;■ 

fci 

6  0 

j 

6 

70 

"7 

if 

80 

-■ 

h. 

90 

3 

7 

00 

*7 

7 

10 

"7 

7 

20 

"7 

i' 

30 

~2 

-> 
r 

40 

~7 

7 

50 

"7 

7 

60 

*7 

7 

70 

IF  11=1  THEN  BEEP 

FNR=I1 

FN  END 

CLEAR  5  DISP  "DO  YOU  WISH  T 

0  FIT  fl  DIFFERENT  TYPE  CURV 

E?<Y/N) " 

INPUT  US 


3780 


IF  U*="Y"  THEN 
GOTO  3740 
GOTO  1020 
NEXT  I 

ASSIGN*  1  TO  * 
BEEP 

CLEAR  15  DISP  " 
ION  COMPLETE" 
END 


'30 


CALIBRAT 
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APPENDIX  J. 
HEAT  TRANSFER  DATA  ACQUISITION  COMPUTER  PROGRAM  LISTING 

:•:>  MASS  STORAGE  IS  ":D701" 

20  SHORT  V(BO) 

30  DIM  C(5).7(80).R<5).L(5).01*C1 00 3 

40  Dl*«"  " 

50  CLEAR  g  DISP  "  QUICK" 

60  DISP  "  " 

70  DISP  "THIS  PROGRAM  IS  FOR  USE  IN  SAMPLING  FOR  FUTURE  PROCESSING  " 

SO  DISP  "OPERATING  PARAMETERS  AND  TZ    OUTPUTS  ARE  STORED  ON  DISK  FOR  FUTURE  REDUC 

TION" 

90  PAUSE 

100  CLEAR  *  DISP  "WHAT  IS  THE  DATE?  (MO/DA/YR)" 

110  INPUT  Dl* 

120  CLEAR  6-  DISP  "WHAT  IS  THE  TIME?  (MILITARY)" 

130  INPUT  T1S 

140  CLEAR  C  DISP  "WHAT  IS  THE  BED  WIDTH?  (INCHES)" 

ISO  INPUT  X 

160  CLEAR  &  DISP  "WHAT  IS  THE  STATIC  BED  HIEGHT?  (INCHES)" 

170  INPUT  SI* 

ISO  CLEAR  &    DISP  "WHAT  WAS  THE  HIGHEST  NUMBERED  PROBE  UTILIZED?  (TWO  DIGITS)" 

190  INPUT  L* 

200  L-YAL(L«) 

210  CLEAR  &  DISP  "WHAT  IS  THE  VOLTAGE  TO  THE  RIGHT-HAND  HEATER?  (VOLTS)"  ' 

220  INPUT  VI 

230  CLEAR  &  DISP  "WHAT  IS  THE  CURRENT  TO  THE  RIGHT-HAND  HEATER?  (AMPS)" 

240  INPUT  II 

250  CLEAR  6-  DISP  "WHAT  IS  THE  VOLTAGE  TO  THE  LEFT-HAND  HEATER?  (VOLTS)" 

260  INPUT  V2 

270  CLEAR  e  DISP  "WHAT  IS  THE  CURRENT  TO  THE  LEFT-HAND  HEATER?  (AMPS)* 

280  INPUT  12 

290  CLEAR  e  DISP  "WHAT  IS  THE  AMBIENT  TEMPERATURE?  (DEG  F) " 

300  INPUT  Al 

310  CLEAR  &  DISP  "WHAT  IS  THE  AIR  FLOW-RATE?  (*/.)" 

320  INPUT  Fl 

330  F2-F1J.373 

340  CLEAR  &  DISP  "WHAT  IS  THE  BED  EXPANSION"?  (INCHES)" 

350  INPUT  E1S 

360  CLEAR  €•  DISP  "ENTER  ANY  COMMENTS  YOU  WISH  RECORDED  (MAX  100  CHARACTERS)" 

370  INPUT  01*C13 

3B0  CLEAR  &  DISP  "WHAT  DO  YOU  WISH  TO  CALL  YOUR  DATA  FILE?  (SIX  CHARACTERS  MAX)" 

390  INPUT  F* 

400  H*-"AF01AL" 

410  H*,C73-"79" 

420  H*C9D«"VT1VR5TE0- 

430  CLEAR  709 

440  OUTPUT  709  ;  H* 

450  CREATE  F«, 1.1500 

460  ASSIGN*  1  TO  F« 

470  ASSIGN*  2  TO  " C0EFF2 : D700 " 

480  PRINT*  1  ;  D1«,T1*,X,S1*,L*.L,V1, II, V2, 12, Al , Fl , F2, El*, 01* 

490  PRINT  USING  5O0 

500  IMAGE  'tttttttXt**ttttttt*XXXtttXtttt" 

510  PRINT  USING  520 

520  IMAGE  3/,3X, "**EXPER I MENTAL  DATA«»" 

530  PRINT  USING  540  ?  F* 

540  IMAGE  BX, "RUN  #",X,6A 

550  PRINT  USING  560  t  D1*,T1* 

560  IMAGE  6X.BA, X, "-", X , 4A 

570  PRINT  USING  580 

580  IMAGE  2/, 2X, "TO*", 3X, "OUTPUT", 5X, "TEMP" 
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ra0    PPIM    USING    600 

rv.-     IMAGE    QX.  "  idiV)  "  .  5X,  "  !DEG    F) "Li 57*00 

elC  PRINT  USING  c20 

b20  IMAGE  2X«" ",3X." ".5X." " 

c~:  CLEAR  >•'  DISP  "PRESS   CONTINUE  WHEN  READV  TO  BEGIN  SAMPLING" 
640  PAUSE 

650  OUTPUT  709  :"TE2" 
60O  POP  1=1  TO  79 
670  OUTPUT  709  :"AS" 
6S0  ENTER  709  :  R 
690  V(I)»R*1000 
700  NEXT  I 

710  OUTPUT  709  ; "TE1TE" 
720  ENTER  709  ;  S 
730  BEEP 

740  CLEAR  (?  DISP  "SAMPLING  COMPLETE" 
750  DISP  "RESULTS  WILL  BE  PRINTED  MOMENTARILY" 
760  READ*  2, 1 
770  FOR  1*1  TO  79 

780  READ*  2  :  C  (  1  )  ,  C  (2)  ,  C  (3)  .  C  <4)  .  C  (5) 
790  IF  IX.  THEN  BIO 
800  GOTO  320 
810  IF  K-67  THEN  940 

820  TCD-Cll  )+C<2)  *V(I)-t-C(3)  *V  <  I )  *V<I)-t-C(4)  *V  ( I )  tV  (  I  )  *V(I)-t-C(5>  *V  ( I  )  *V(I  )  iV  <  I  )  tv 
(I) 

830  T(I)»T (I) f 9/5+32 
840  PRINT*  1  ;  V(I)  ,T(I) 
850  IF  1*41  THEN  920 
860  PRINT  USING  870 

870  IMAGE  3/,2X. "TC*".3X. "OUTPUT ",5X, "TEMP" 
880  PRINT  USING  890 
890  IMAGE  9X.  " (mV)  n.5X,  "  (DEG  F)  " 
900  PRINT  USING  910 

910  IMAGE  2X,"- — ".3X," ",5X." " 

920  PRINT  USING  930  ;  I .  V  ( I )  ,  T  ( I ) 

930  IMAGE  2X,DD,3X.DD.DDD.4X,DDD.DDDD 

940  NEXT  I 

950  ASSIGN*  1  TO  * 

960  ASSIGN*  2  TO  * 

970  R<1)«(T<1)+T(2)+T<3) >/3 

980  R(2)-(T(5)-t-T(6)+T(7)+T(8)  )  /4 

990  R(3)»(T(9)+T(10)+T(11)+T<12) ) /4 

1000  R<4)*(T<13>-t-T<14)+T<15)+T<16>  )  /4 

1010  R(5)-<T<17)-t-T  (18>+TU9)-t-T(20>  )  /4 

1020  L(1)=(T(21)+T(22)+T(23)+T(24) )/4 

1030  L(2)  =  (T<25)-t-T<26)-t-T(27)-t-T(28)  )  /4 

1040  LC3)-(T<29>«-TC30>+T<31>+T<32>  )  /4 

1050  L<4>-<TC33)+T(34>+T<35)+T(36>  > /4 

1060  L(5)"(T(37)-t-T(3B)-t-T(39)+T<40))/4 

1070  Gl-(T<41)t-T(42)+T<43)  )/3 

1080  G2-(T(44)+T(45)+T(46) )/3 

1090  Pl-<T<47>t-T(48>+T<49)  ) /3 

1100  Bl-0 

1110  Xl-0 

1120  FOR  1-50  TO  L 

1130  IF  1-55  THEN  1150 

1140  GOTO  1170 

1150  Xl-1 

1160  GOTO  1260 

1170  IF  1-65  THEN  1190 


189 


:  i  bo 

SOTO 

1210 

::  =o 

A  1  =2 

1 200 

GOTO 

12=0 

•  —•  i  ' 

IF  I =oo  THEN  1230 

■  ~2C 

GOTO  . 

1250 

1 230 

XI  =3 

124C 

3CT0 

1260 

1  250 

B1=B1- 

►T(I) 

1260 

NEXT 

r 

1270 

ni=l-: 

30+1 -XI 

1280 

H1  =  (R  (  1)  +R  (2)  -t-R<3)+R(4)  +R  (5)  )  /5 

1290 

H2=  (L  <  1  )  -H_  (2)  -H_  (3)  +L  (4) +L (5)  )  /5 

1  300 

B1=B1/N1 

1310 

E1=V1*I1*3.4121 

1320 

E2=V2i 

H2*3.4121 

1330 

if  xi : 

>1  THEN  1360 

1 340 

P2=T(55) 

1350 

GOTO  : 

1400 

I3o0 

if  xi: 

:2  THEN  1390 

1370 

P2=(T(55)+T (65) ) /2 

1380 

GOTO  : 

1400 

1390 

P2=(T  (55)  +T  (65) -^T  (66)  )  /3 

14O0 

PRINT 

USING  1410  :  X 

1410 

IMAGE 

3/, "BED  WIDTH  =  " . DD. DD, X , " IN" 

1420 

PRINT 

USING  1430  :  SIS 

1430 

IMAGE 

"STATIC  BED  HEIGHT  -  ",5A."  IN" 

1440 

PRINT 

USING  1450  ;  F2 

1450 

IMAGE 

"AIR  FLOW  RATE  *  " , DDD.DD."  CFM 

1460 

PRINT 

USING  1470  ;  Al 

1470 

IMAGE 

" AMB I ENT  TEMPERATURE-  " . DD . DD , " 

1480 

PRINT 

USING  1490  :  VI 

1490 

IMAGE 

"RIGHT  HEATER  VOLTAGE-  ",DD.D," 

1500 

PRINT 

USING  1510  ;  11 

1510 

IMAGE 

"RIGHT  HEATER  CURRENT  -  ",D.DD. 

1520 

PRINT 

USING  1530  :  V2 

1530 

IMAGE 

"LEFT  HEATER  VOLTAGE  -  " .DD.D." 

1540 

PRINT 

USING  1550  :  12 

1550 

IMAGE 

"LEFT  HEATER  CURRENT  =  ",D.DD," 

1560 

PRINT 

USING  1570  :  El* 

1570 

IMAGE 

"BED  EXPANSION  -  " . 4A. "  INCHES" 

1580 

PRINT 

USING  1590  :  01* 

1590 

IMAGE 

"COMMENTS:  ",100A 

6160 

END 

DEG  F' 


VOLTS' 


AMPS' 


VOLTS' 


AMPS" 
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APPENDIX  K. 
HEAT  TRANSFER  DATA  REDUCTION  COMPUTER  PROGRAM  LISTING 


10  SHORT  VOO) 

20  PRINTER  IS  701.  100 

21  PRINT  CHR«<27)&"StllL" 

30  *>r:nt  chr*  c2T)  vsa  oD" 

40  DIM  C (5) .T (80) ,R(5) ,L(5) ,01*C1003 

50  Ol*»"  " 

60  E<*«0 

70  B9=0 

80  CLEAR  i?  DISP  "  PROCESS" 

90  DISP  "  " 

100  DISP  "THIS  PROGRAM  IS  FOR  USE  IN  PROCESSING  DATA  WHICH  HAS  BEEN  PREVIOUSLY  R 

ECORDED. " 

110  DISP  "  " 

120  DISP  "HEAT  TRANSFER  COEFFICIENTS  ARE  CALCULATED,  DISPLAYED,  AND  PRINTED." 

130  PAUSE 

140  CLEAR  <?  DISP  "WHAT  IS  THE  NAME  OF  THE  DATA  FILE?  (SIX  CHARACTERS  MAX)" 

130  INPUT  F* 

160  CLEAR  e  DISP  "PRESS  CONTINUE  TO  BEGIN  PROCESSING" 

170  PAUSE 

180  ASSIGNt  1  TO  F» 

190  READ*  1  ;  D1S, Tl», X, SIS. L*,L, VI , 1 1 , V2. 12, Al , Fl , F2, El*, 01* 

200  F2-F1*. 32454*. 793135 

210  CLEAR  &  DISP  "DO  YOU  WISH  A  PRINT  OUT  OF  THE  EXPERIMENTAL  DATA?  CY/N) " 

220  INPUT  RB« 

230  IF  R8*-"Y"  THEN  260 

240  B9-1 

230  GOTO  1110 

260  CLEAR  e  DISP  "         COMMENT  STATEMENT:" 

270  DISP  "  " 

280  DISP  01* 

290  DISP  "  " 

300  DISP  "DO  YOU  WISH  TO  MAKE  ANY  CORRECTIONS?  <Y/N) " 

310  INPUT  R9* 

320  IF  R9*-"N"  THEN  360 

330  DISP  "  " 

340  DISP  "ENTER  CORRECTED  COMMENT  STATEMENT" 

350  INPUT  01 « 

360  ! 

370  ! 

380  !       ***PRINT  EXPERIMENTAL  CONDITIONS*** 

390  ! 

400  PRINT  CHR* (12) 

410  PRINT  USING  420 

420  IMAGE  3/,30X, "**EXPER I MENTAL  DATA**" 

430  PRINT  USING  440  ;  F* 

440  IMAGE  35X,"RUN  #",X.6A 

450  PRINT  USING  460  |  Dl*,Tl* 

460  IMAGE  33X.8A.X, "-",X,4A 

470  IF  E9-1  THEN  1030 

480  PRINT  USING  490  j  F* 

490  IMAGE  5/,22X,"RUN  NUMBER  »",20X,6A 

500  PRINT  USING  510  ;  X 

510  IMAGE  /,22X,"BED  WIDTH  -", 14X,DD. DD, X, "INCHES" 

520  PRINT  USING  530  ;  SI* 

530  IMAGE  /,22X, "STATIC  BED  HEIGHT  «". 6X, 5A, X, "INCHES" 

540  PRINT  USING  550  5  F2 

330  IMAGE  /,22X,"AIR  FLOW  RATE  -*' ,  12X,  DDD.DD.  X,  "CFM" 

560  PRINT  USING  370  j  Al 
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f  =  ; 

a  X 

=  .  : 

=  2'- 

b3C 

=.40 

650 

soO 

670 

690 

690 

700 

710 

720 

730 

740 

750 

760 

-70 

7S0 

790 

900 

910 

920 

930 

940 

950 

960 

970 

990 

990 

900 

910 

920 

930 

9A0 

950 

960 

970 

990 

O90 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1090 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 


"EMPERATURE  =' 


AMBIEN"! 
590  :  v; 

"RIGHT  ~i£— — E?  ^C_ "AGE  =' 

si 0  :  I  '. 

"RIGHT  HEATEF  CURRENT  =  •*  ,  6X.  D.  DD.  X.  "AMPS" 

630  :  -Z 

"LEFT  HEATER  VOLTAGE  ="  .  oX  .  DD.  2.  X  .  "VCL.Tg  " 
650  :  12 

'VET-  HEA~En  2URRENT  =  " . 7X , D . DD . X . - AMPS" 

670  :  El* 

"BED  EXPAN9I0N  =". 1 1 X . 4A. X ." INCHES" 


I  '•liiGc   .  22  X 
PRINT  USING 

:--3£   .22; 

PRINT  USING 

:*iAGE   .  22x 

PRINT  _s:ng 

IMAGE  • . 22X 

PRINT  USING 

I-AGE   .  22x 

PRINT  USING 

I  MAGE  /  .  22  X 

D9«LEN  [01  s 

□1SCD9+1.  100:  =  "  " 

IF  D9>30  THEN  740 

PRINT  USING  720  ;  01SCI.D93 

IMAGE  /.22X. "COMMENTS: ", X.30A 

GOTO  970 

K9-25 

IF  01«CK9,K93«a  "  THEN  790 

k'9«K9-l 

GOTO  750 

PRINT  USING  790  :  01SC1.K93 

IMAGE  /.22X. "COMMENTS: ". X.40A 

D2-K9+35 

IF  D9>D2  THEN  970 

K9-K9*l 

IF  K9>D9  THEN  970 

PRINT  USING  950  :  01SCK9-1.D93 

IMAGE  22X.50A 

GOTO  970 

;_9=D2 

IF  C1SCL9,L9D«" 

L9-t_9»l 

GOTO  990 

PRINT  USING  920 

IMAGE  22X.50A 

L8-L9+1 

IF  LS>D9  THEN  D9=L9 

PRINT  USING  960  :  01 «CL9+1 . D93 

IMAGE  22X.50A 

PRINT  CHRS(12) 

E9-1 

GOTO  360 


;x.  DD.  DD.  X.  "DEG  F' 
. 5X.DD.D.  K,  "  -'3_~: 


THEN  910 


01SCK9*1.L9D 


***PRINT  THERMOCOUPLE  READINGS*** 

I 

»RINT  USING  1040 

IMAGE  /. 30 X. "THERMOCOUPLE  READINGS" 

=>RINT  USING  1060 

IMAGE  /, 13X, "TC#",3X. " OUTPUT ". 4 X , "TEMP", 14X, "TC#".3X, "OUTPUT", 4X, "TEMP" 

PRINT  USING  1090 

IMAGE  20X,  "  (mV) ",4X,  " (DEG  F) " .  19X , " (mV) " , 4X . "  (DEG  F)  " 

PRINT  USING  1100 

IMAGE  13X.  " ".3X,  " ",4X,  " ",  14X,  " U,3X.  " "  .  4X,  " " 

FOR  1-1  TO  79 

IF  I>L  THEN  1140 

GOTO  1150 

IF  K-67  THEN  1160 

READ*  1  :  V ( I ) , T ( I > 

NEXT  I 
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1170  IF  B9=l  THEN  1390 

: : 30  wi=o 

1 1 90  Z 1 =0 

1200  FOR  (  =1  TO  40 

1210  IF  Wl  =  l  THEN  J-J+l  6'  GOTO  1230 

1220  J»K+40 

1230  IF  J--L  THEN  1250 

1240  GOTO  1260 

1250  IF  J -.=67  THEN  1340 

1260  IF  J :79  THEN  Z 1 = 1 

1270  IF  Zl=l  THEN  1310 

1280  PRINT  USING  1290  :  K . V <K>  . T (K)  , J . V ( J )  . T  ( J ) 

1290  IMAGE  13X.DD.3X.DD.DDD.4X.DDD.DD. 13X . DD, 3X . DD. DDD. 4X . DDD. DD 

1300  GOTO  1330 

1310  PRINT  USING  1320  :  K,V(K),T(K) 

1320  IMAGE  13X.DD.3X, DD. DDD, 4X . DDD. DD 

1330  GOTO  1360 

1340  J=Jt-l  6'  Wl*l 

1350  GOTO  1250 

1360  NEXT  K 

1370  ASSIGN*  1  TO  » 

1380  PRINT  CHRSU2) 

1390  R(1>=(T<1)*T (2)+T(3) ) /3 

1 400  R(2)*(T(5>+T(6)+T(7)+T(8))/4 

1410  R(3)«<T<9>+T(10)^T(H)+T(12)  )/4 

1420  R<4)  =  (T(13>+T(14)-t-T(l5)-»-T<16)  )  /4 

1430  R(5)=(T(17)+T(18>+T<19>+T(20) ) /4 

1440  L(l)  =  (T(21)«-T(22)-t-T(23)+T(24)  )  /4 

1450  L<2)»(T(25)+T(26)+T(27)-t-T(28)  )  /4 

1460  L(3)=(T(29)+T(30)+T(31)*T(32) ) /4 

1470  L<4)*(T(33)+T(34)+T(35)-t-T<36) ) /4 

1480  L(5)=(T(37)+T (38) *T (39) +T (40) ) /4 

1490  Hl-(R(l)-t-R(2)t-R(3)fR<4)fR(5)  )  /5 

1500  H2«(L(1)-H_(2)+L(3)-M_(4)t-L(5)  ) /5 

1510  IF  B9-1  THEN  3500 

1520  ! 

1530  !        «**DISPLAY  HEATER  TEMP  SNAPSHOT (R-H) *** 

1540  ! 

1550  PRINT  USING  1560 

1560  IMAGE  3/,30X. "**EXPERIMENTAL  DATA*!" 

1570  PRINT  USING  1580  ;  F« 

1580  IMAGE  35X. "RUN#".X,6A 

1590  PRINT  USING  1600  ;  D1*,T1* 

1600  IMAGE  33X.8A, X, "-", X , 4A 

1610  PRINT  USING  1620 

1620  IMAGE  3/, 32X, "RIGHT-HAND  HEATER" 

1630  PRINT  USING  1640 

1640  IMAGE  31 X, "TEMPERATURE  PROFILE" 

1650  PRINT  USING  1660 

1660  IMAGE  /.25X,"C0L  1".3X,"C0L  2",3X,"C0L  3".3X,"C0L  4" , 7X, "AVE" 

1670  PRINT  USING  1680 

1680  IMAGE  23X," " 

1690  PRINT  USING  1700  :  T ( 1 ) , T(2) . T (3) ,R ( 1 ) 

1700  IMAGE  17X,"R0W  1 " , X, " ! " . X. DDD. DD, 2X, DDD.DD.2X, DDD. DD, 2X. "**X*S*" , 2X, M ! " , 2X . 

DDD.DD 

1710  PRINT  USING  1720 

1720  IMAGE  23X, " ! ",33X. "! " 

1730  PRINT  USING  1740  ;  T  (5)  ,  T  (6)  .  T  (7)  .  T  (8)  ,  R  (2) 

1740  IMAGE  17X,"R0W  2"  ,  X,  "  !  "  ,  X,  DDD.  DD.  2X,  DDD.  DD,  2X,  DDD.  DD,  2X  .  DDD.  DD,  2X,  "'  i  "  ,  2X  .  DD 

D.DD 
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1"*50  PRINT  USING  1760 

'.  760  I  MAGE  23X  ,  "  !  "  ,  ZZV. .  "  .'  " 

: — 0  ='RINT  USING  1730  :  T(9)  .7(10!  ,T(11)  .7(12)  ,RC) 

1~30  IMAGE  17X.  "ROW  3".X. " i ",  X.  DDD. DD. 2X , DDD.  DD. 2X . DDD. DD.  2X.DDD.DD.  2X.  "  i  "  .  2>  .  DD 

D.DD 

1790  PRINT  USING  1300 

1  BOO  I MAGE  23X ,  "  !  "  , ZZX ,  "  !  " 

1810  PRINT  USING  1820  ;  T \ 13) , T < 14) , T ( 15) , T < 1 6) , R (4) 

1820  IMAGE  17X."R0W  4 " . X . " I " . X . DDD . DD. 2X . DDD. DD. 2X . DDD . DD. 2X . DDD. DD . 2X . " ! " . 2> . DC 

D.DD 

1830  PRINT  USING  1840 

1840  IMAGE  23X. " I". 33X. "I " 

1850  PRINT  USING  1860  ;  T(17),T(1B),T(19).T(20).R(5) 

1860  IMAGE  17X,"R0W  5" . X , " : " . X . DDD. DD, 2X , DDD. DD. 2X . DDD. DD, 2X . DDD. DD, 2X . " i " . 2X . DD 

D.DD 

1870  PRINT  USING  1880 


1880  I  MAGE  23  X  ,  " " 

1890  PRINT  USING  1900 

1900  IMAGE  /, 16X, "NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT" 

1910  PRINT  USING  1920 

1920  IMAGE  22 X, "AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED" 

1930  PRINT  USING  1940 

1940  IMAGE  21/. 17X. "NOTE:  PLOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH" 

1950  PRINT  USING  1960 

1960  IMAGE  23X, "HORIZONTAL  ROW" 

1970  PRINT  CHR*(12) 

1 930  ! 

1990  !   ***PLOT  TEMP  GRAPH (R-H) ** 

2000  ! 

2010  PRINTER  IS  2 

2020  PRINT  USING  2030  :  F* 

2030  IMAGE  4/.10X."RUN#  " . X , 6A 

2040  PRINT  USING  2050 

2050  IMAGE  3/ 

2060  PRINTER  IS  701,100 

2070  Z9=IP(R(3) ) 

2030  ZB=RMD(Z9.5) 

2090  IF  Z8=0  THEN  2120 

2100  Z9=Z9-1 

2110  GOTO  2080 

2120  Z6-Z9-15 

2130  Z7»Z9+15 

2140  GCLEAR  6"  SETGU 

2150  LOCATE  25,118.20.30 

2160  !  FRAME 

2170  SCALE  Z6.Z7.0.  10 

2130  FXD  0 

2190  LAXES  1,1. Z6. 0,5,1 

2200  Z5-Z7-16 

2210  MOVE  Z5,-3 

2220  DEB 

2230  LDIR  0  3  LORG  4 

2240  LABEL  "TEMP  (DEB  F) " 

2250  Z4-Z6-3 

2260  MOVE  Z4,5 

2270  LDIR  90  3  LORG  4 

2230  LABEL  "VERT  POSITION  (IN)" 

2290  Z3=Z6+14 

2300  MOVE  Z3, 12 

2310  LDIR  0  3  LORG  4 
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•40 


2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
25*0 
2600 
2610 
2620 
2630 
D.DD 
2640 
2650 
2660 
2670 
D.DD 
2680 
2690 
2700 
2710 
D.DD 
2720 
2730 
2740 
2750 
D.DD 
2760 
2770 
2780 
2790 
D.DD 
2800 
2810 
2820 
2830 
2840 
2850 
2860 


LhBEl.  "RIGHT -HAND  HEATER" 

12=:=-: 4 

MOVE  12, 1 : 

LDIR  0  <•■    L3RG  4 

i-ABE-  "TEMC  vs  DEPTH" 

SETUU 

Y=9 

FOR  1=1  TO  5 

MOVE  R  (  I  )  . Y 

Y=Y-2 

LORG  5 

LABEL  "«" 

NEXT  I 

COPY 

!      ***DISPLAY  HEATER  TEMP  SNAPSHOT  <L-H)<»* 

PRINT  USING  2490 

IMAGE  3/.30X, "**EXPER I MENTAL  DATA*«" 

PRINT  USING  2510  :  F$ 

IMAGE  35X, "RUN#" . X.6A 

PRINT  USING  2530  :  D1*,T1* 

IMAGE  33X, 8A. >."-". X.4A 

PRINT  USING  2550 

IMAGE  3/. 33X. "LEFT-HAND  HEATER" 

PRINT  USING  2570 

IMAGE  31 X. "TEMPERATURE  PROFILE" 

PRINT  USING  2590 

IMAGE  /.25X."C0L  1",3X,"C0L  2".3X."C0L  3",3X,"C0L  4".7X,"AVE" 

PRINT  USING  2610 

IMAGE  23 X." " 

PRINT  USING  2630  ;  T < 21 ) . T (22) , T (23) , T (24) . L < 1 ) 

IMAGE  17X. "ROW  1",X. " I " , X , DDD. DD. 2X , DDD. DD, 2X . DDD. DD, 2X . DDD. DD, 2X . 

PRINT  USING  2650 

IMAGE  23X. "  ! ",33X,  "  !  " 

PRINT  USING  2670  ;  T (25) , T (25)  , T (27) , T (28)  . L  (2) 

IMAGE  17X. "ROW  2",X, " ! " , X . DDD. DD, 2X , DDD. DD, 2X , DDD. DD, 2X , DDD. DD, 2X . 

PRINT  USING  2690 

IMAGE  23X. " ! ",33X,  " I " 

PRINT  USING  2710  :  T (29) , T (30) , T (31 ) , T (32) . L (3) 

IMAGE  17X, "ROW  3",X, " ! " , X . DDD. DD, 2X , DDD. DD. 2X , DDD. DD. 2X . DDD. DD, 2X . 

PRINT  USING  2730 

IMAGE  23X,  " ! ",33X.  "  !  " 

PRINT  USING  2750  ;  T  (33)  ,  T  (34)  .  T  (35)  ,  T  (36)  .  L  (4) 

IMAGE  17X, "ROW  4".X, " ! " , X , DDD. DD, 2X , DDD. DD. 2X , DDD. DD, 2X . DDD. DD, 2X . 

PRINT  USING  2770 

IMAGE  23X , " ! " , 33X ,  n !  " 

PRINT  USING  2790  :  T  (37)  ,  T  (38)  .  T  (39)  .  T  (40)  ,  L  (5) 

IMAGE  17X, "ROW  5", X, " I " . X . DDD. DD,2X , DDD. DD, 2X, DDD. DD. 2X . DDD. DD,2X, 

PRINT  USING  2810 

IMAGE  23 X." " 


2X.DD 


2X.DD 


.2X.DD 


.2X.DD 


',2X,DD 


PRINT  USING  2830 

IMAGE  /. 16X, "NOTEs  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT- 
PRINT  USING  2850 

IMAGE  22 X, "AND  ARE  ARRANGED  AS  SEEN  FROM  WITHIN  THE  BED" 
PRINT  USING  2870 
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2370  IMAGE  21/ ,  17X.  "NOTE:  ^'LOTTED  TEMPERATURES  ARE  AVERAGES  OF  EACH' 

2330  PRINT  USING  2390 

2390  IMAGE  23X. "HORIZONTAL  ROW" 

29O0  PRINT  CHRS(12> 

2-10  ' 

2-20  !      t  *='LOT  TEMP  GRAPH  <L-H)  ».*  " 

2930  ! 

2C40  PRINTER  IS  2 

295i."'  PRINT  USING  2960  :  Fs 

2^60  IMAGE  4/. 10 X. "RUN*  " , X . 6A 

2970  PRINT  USING  2980 

2<?30  IMAGE  3/ 

2990  PRINTER  IS  701,100 

3000  Z9=IP(L(3> ) 

3010  Z3=RMD(Z9.5) 

3050 


3020 

IF  Z8=0  THEN 

3030 

Z9=Z9-1 

3040 

GOTO  3010 

3050 

Z6=Z9-15 

3060 

Z7=Z9+15 

3070  GCLEAR  ?  SETGU 

3030  LOCATE  25.113,20,30 

3090  '  FRAME 

3100  SCALE  Z6.Z7.0,10 

3110  FXD  0 

3120  LAXES  1.1, Z6, 0,5.1 

3130  Z5=Z7-16 

3140  MOVE  Z5.-3 

3150  DEG 

3160  LDIR  0  6'  LORG  4 

3170  LA8EL  "TEMP  (DEG  F) " 

3130  Z4=Z6-3 

3190  MOVE  Z4.5 

3200  LDIR  90  &  LORG  4 

3210  LABEL  "VERT  POSITION  (IN)" 

3220  Z3=Z6+14 

3230  MOVE  Z3, 12 

3240  LDIR  0  &  LORG  4 

3250  LABEL  "LEFT-HAND  HEATER" 

3260  Z2=Z6+14 

3270  MOVE  Z2, 11 

3230  LDIR  0  &    LORG  4 

3290  LABEL  "TEMP  vs  DEPTH" 

3300  SETUU 

3310  Y=9 

3320  FOR  1=1  TO  5 

3330  MOVE  L ( I ) . Y 

3340  Y=Y-2 

3350  LORG  5 

3360  LABEL  "*" 

3370  NEXT  I 

3330  COPY 

3390  !  PRINT  PAGE  HEADING 

34O0  PRINT  USING  3410 

3410  IMAGE  3/.30X. "**EXPER I MENTAL  DATA»*n 

3420  PRINT  USING  3430  :  F* 

3430  IMAGE  35X , "RUN#" , X, 6A 

3440  PRINT  USING  3450  ;  D1*,T1$ 

3450  IMAGE  33X , 3A, X , "-" , X , 4A 

3460  PRINT  USING  3470 
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3470 
3430 
-.  4.90 

3500 
35 1 0 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
36 1  0 
3620 
3630 
3640 
3650 
3660 
3670 
36S0 
3690 
3700 
3710 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4030 
4O40 
4050 
4060 


IMAGE  2-'.  34  X.  "FLU  ID  I  ZED  BED" 

PRINT  USING  3490 

IMAGE  3 1 X. M TEMPERATURE  PROFILE" 

READ  BED  SIZE 
Q1*=P*l1 .23 
IF  Gl*="12" 
IF  Gl*="10" 
IF  G14="08" 


THEN  4590 
THEN  4220 
THEN  38('0 


IF  B9=l  THEN  3850 

!      **DISPLAY  6  IN  BED  TEMP  SNAPSHOT** 
i 

PRINT  USING  3600 

PRINT  USING  3620 

IMAGE  27X, " 

PRINT  USING 

IMAGE  27X, " 

PRINT  USING 

IMAGE  27X. " I ". X.DDD.DD.2X,DDD.DD.2X.DDD.DD.2X. " ! 

PRINT  USING 

IMAGE  27X. " 

PRINT  USING 

IMAGE  27X. " 

PRINT  USING 

IMAGE  27 X. " 

PRINT  USING 

IMAGE  27 X. " 

PRINT  USING 

IMAGE  27 X.  " 

PRINT  USING 

IMAGE  27X. " I ".X.DDD.DD.2X.DDD.DD.2X,DDD.DD,2X, " ! 

PRINT  USING 

IMAGE  27X, " 

PRINT  USING  3830 

PRINT  USING  3840 

IMAGE  27X, " ! ",25X, " ! 

IMAGE  27X,  "• 

!  CALCULATE  AVERAGE  BED  TEMP 

Bl»  (T  (50) +T  (51 )  *T  (52)  t-T  (53) +T  (54) -t-T  (56) +T  (57)  )  /7 

GOTO  4900 


".25X, 

"  ;  " 

3640 

".25X, 

" :  •• 

3660  : 

T(56) . 

T(53) 

.T( 

50) 

". X.DDD.DD.2X.DDD. 

DD. 

2X, 

DDD. 

DD. 

3680 

" . 25 X , 

" :  •• 

3700 

".25X, 

H .  H 

3720  ; 

T(51) 

" , 1 7X , 

DDD.DD, 

2x. " : 

" 

3740 

" ,  25X, 

" :  ■• 

3760 

",25X. 

3780  : 

T(57) . 

T(54) 

,T( 

52) 

",X,DDD.DD.2> 

..DDD. 

DD, 

2X, 

DDD. 

DD, 

3800 

M,25X, 

•• ;  •• 

IF  B9»l 


**DISPLAY  8  IN  BED  TEMP  SNAPSHOT** 
THEN  4180 


PRINT  USING  3930 

IMAGE  2/,25X, " 

PRINT  USING  3950 


IMAGE  25 X, " 
PRINT  USING 
IMAGE  23X," 
PRINT  USING 
IMAGE  25X, " 
PRINT  USING 
IMAGE  25X, " 
PRINT  USING 
IMAGE  25X, " 
PRINT  USING 
IMAGE  25 X, " 


".29X,  •• : " 

3970 

".29X,".'" 

3990  :  T(56) ,T(53) ,T(50) 

", 5X. DDD. DD.2X, DDD. DD.2X. DDD.DD 

4010 

",29X,  •• :  ■• 

4030 

" .29X, " ! " 

4050  ;  T(51) 

",21X,DDD.DD,2X, " ! " 


,2X, 


PRINT  USING  4070 
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-•.'  >.' 
-030 
iO=»0 

- :  c»o 

4110 

4 1 20 

4  1  30 

-J  1  40 

4150 

4160 

4170 

4130 

4  190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4230 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

44O0 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

+T(6: 

4620 

4630 

4640 

4650 


57  -  .  T  <  54 )  .  T  < 52 ) 

DD.2X.DDD. DD.2X.DDD.DD.2X, 


iriAGE  25X.  "  :  " .  2°>. . 

PRINT  USING  4090 

IMAGE  25X. " ! " , 29X. 

PRINT  USING  4110  :  ~C 

IMAGE  25X.  '•  :  ".  5X.  DDD. 

PRINT  USING  4130 

IMAGE  25X. " ! " .29X. " 

PRINT  USING  4150 

IMAGE  25X. " ! ".29X. " 

PRINT  USING  4170 

1  CALCULATE  AVERAGE  BED  TEMP 

B 1 =  <  T  <  50  >  +T  <  5 1 ) +T ( 52 ) +T  <  53 ) +T ( 54  >  +T  <  56 ) +T ( 57 )  ) /7 
GOTO  4900 

!        ***DISPLAY  10  IN  BED  TEMP  SNAPSHOT*** 
IF  B9=l  THEN  4550 


PRINT  USING  4260 

IMAGE  2/.23X. " 

PRINT  USING  4280 


IMAGE  23X, " 
PRINT  USING 
IMAGE  23X. " 
PRINT  USING 
IMAGE  23X. " 
PRINT  USING 
IMAGE  23 X. " 
PRINT  USING 
IMAGE  23X, " 
PRINT  USING 
IMAGE  23X, " 
PRINT  USING 
IMAGE  23X, " 
PRINT  USING 
IMAGE  23X. " 
PRINT  USING 
IMAGE  23X, " 
PRINT  USING 
IMAGE  23 X, H 
PRINT  USING 
IMAGE  23X, " 
PRINT  USING  4500 
IMAGE  23 X. " 


«  -^-^X  «  "  i 
4300 

" . 33X. "  !  " 

4320  ;  T (58) ,T (56) .T(53) ,T(50) 

",X,DDD.DD.2X.DDD.DD.2X,DDD.DD.2X.DDD.DD.2X. " ! " 
4340 

" . 33X , " : " 
4360 

" , 33X, " ! " 

4380  ;  T(59) ,T(51) 
",X,DDD.DD. 18X,DDD.DD,2X. " I " 
4400 

" i 33X. "  !  " 
4420 

" .  zzx  .'•:•' 

4440  :  T(60) ,T(57) ,T(54) ,T(52) 

".X.DDD.DD.2X,DDD.DD.2X,DDD.DD,2X,DDD.DD.2X. " 1 " 

4460 

" . 33X , " ! " 

4480 

■-33X.-I" 


PRINT  USING  4520 

IMAGE  /. 15X. "NOTE:  TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT" 

PRINT  USING  4540 

IMAGE  2 IX, "AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED" 

!  CALCULATE  AVERAGE  BED  TEMP 

Bl-<T(50)+T(51>*T(52)+T(53)+T(54)-t-T<56)+T(57)+T(58>t-T  (59)  t-T  (60)  )  /  10 

GOTO  4950 

I 

!      ***DISPLAY  12  IN  BED  TEMP  SNAPSHOTttt 

!  CALCULATE  AVERAGE  BED  TEMP 

Bl-<T(50)+T<51)+T(52)*T(53)*T<54)fT(56)t-T(57)+T<58)-t-T(59)+T(60)*T(61)+T<62) 

)  )/13 

IF  B9-1  THEN  4950 

i 

PRINT  USING  4650 

IMAGE  2/.19X," " 
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4660 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 
4  750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5110 
5120 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 


PRINT 

IMAGE 

PP.  I  NT 

IMAGE 

PR  I  NT 

IMAGE 

PRINT 

IMAGE 

PRINT 

IMAGE 

PRINT 

IMAGE 

PRINT 

IMAGE 

PRINT 

IMAGE 

PRINT 

IMAGE 

PRINT 

IMAGE 

PRINT 

IMAGE 

PRINT 

IMAGE 

IF  B9» 

PRINT 

IMAGE 

PRINT 

IMAGE 


"(61)  .T  (58)  .  T(5e>  .TC 
DD.2X.DDD.DD.2X.DDD. 


DD.2X. 


iO) 

DDD. 


T (62) .TC59) .T(51 ) 
DD.2X.DDD.DD. 18X, 


DDD.DD.2X. 


USING 
19X.  " 
USING 
19)  .  " 
USING 
19X.  " 
USING 
19X,  " 
USING 
19X  .  " 
USING 
19X.  " 
USING 
19X.  " 
USING 
19X,  " 
USING 
19X,  " 
USING 
19X.  " 
USING 
19X.  " 
USING 

19X. " 

'1  THEN  4960 

USING  4920 

/, 15X. "NOTE 

USING  4940 

2 IX. "AND  ARE  ARRANGED  AS  SEEN  FROM  ABOVE  THE  BED" 


4o7  0 

". 41X.  " 

4o90 

"  .  4  1  X  .  " 

4710  : 

".X, DDD 

4730 

"  .  4  1  X  .  " 

4750 

",41X, " 

4770 

" . X . DDD 

4790 

".41X.  " 

4810 

"  .  4  1  X  ,  " 

4830 

" . X , DDD 

4850 

"  .  4  1  X  .  " 

4870 

" . 41X. " 

4890 


DD.2X.DDD.DD.2X. 


T ( 63 )  . T ( 60 ) , T ( 57 )  . T ( 54 ) , T  <  52 ) 
DD ,  2X  .  DDD .  DD ,  2X  .  DDD . DD . 2X . DDD . DD , 2X , DDD . DD . 2X . 


TEMPERATURES  SHOWN  ARE  IN  DEGREES  FAHRENHEIT' 


***CALCULATE  RESULTS** 


1  CALCULATE  AVERAGE  AIR  INLET  TEMP 
P1»(T(47)+T(48)+T(49) ) /3 
1  CALCULATE  AIR  OUTLET  TEMP 
P2=T(55) 

!  CALCULATE  RH  HEATER  GUARD  TEMP 
Gl=(T(41)+T(42)-t-T(43)  )  /3 
!  PRINT  PAGE  HEADING 
PRINT  CHR«(12) 
PRINT  USING  5070 

IMAGE  3/.30X. "**CALCULATED  RESULTS**" 
PRINT  USING  5090  ;  F« 
IMAGE  35X ,  " RUN  # " , X . 6A 
PRINT  USING  5110  !  D1S.T1* 
IMAGE  33X.BA.X, "-" ,X,4A 
PRINT  USING  5130 

IMAGE  2/, 32X, "RIGHT-HAND  HEATER" 
!  PRINT  AVERAGE  HEATER  TEMPERATURE 
PRINT  USING  5160  :  HI 

IMAGE  2/, 15 X, "AVERAGE  HEATER  TEMPERATURE  -" , 12X . DDD. DD, X, "DEG  F" 
!  CALCULATE  AND  DISPLAY  ELECTRICAL  ENERGY  INTO  HEATER 
E1-V1*I1*3.4121 

!  CALCULATE  AND  DISPLAY  LOSS  OUT  OF  BACK  OF  HEATER 
G3«(H1-G1) *. 1703163 

!  CALCULATE  CONVECTION  HEAT  TRANSFER  COEFFICIENT  FROM  GUARD  TO  ATMOSPHERE 
C1=G3/ ( . 347222* (Gl-Al ) ) 

!  CALCULATE  LOSS  FROM  SIDES  OF  HEATER  TO  ATMOSHPERE 
S3-C1*. 069444* (T(72)-A1) 
!  CALCULATE  LOSS  FROM  TOP  STRIP  TO  ATMOSPHERE 
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5260  '_iC=Cl  *  .  0208333*  ■  T  ■  oo  <  -Al  ) 

ZZ.~r.  ZA,_CU_ATE  LOSS  FROM  BCT""Qr«  STRIP  TC  ATMOSPHERE 

5230  S7=C1*.020S333»  ■"■'"  I'-Ai) 

52*0  !  CALCULATE  ~~TA,_  _0S3ES  TC  ATMOSPHERE 

5300  r-4=53-53-^L2*S7 

5310  !  CALCULATE  HEAT  FLUX  THRU  RH  PLATE 

5320  P3"E1-H4 

5330  V9«0 

5340  '  CALCULATE  SURFACE  TEMPERATURE  Or  PLATE 

5350  C'5=H1-P3*. 0208333/ ( (232-. 032* (HI -70) >*. 347222) 

5360  !  CALCULATE  HEAT  TRANSFER  COEFFICIENT  FOR  PLATE  TO  BED 

5370  P7=P3/ (.347222* (P5-B1) > 

5380  '  CALCULATE  LOSS  FROM  SIDE  STRIPS  TO  BED 

5390  C3=P7*. 069444* (T (73) -Bl) 

5400  '  CALCULATE  LOSS  FROM  TOP  STRIP  TO  BED 

5410  S9=P7*. 0208333* (T(69)-B1) 

5420  !  CALCULATE  LOSS  FROM  BOTTOM  STRIP  TO  BED 

5430  S5=P7*. 0208333* (T(70)-B1) 

5440  '  CALCULATE  TOTAL  LOSSES  INTO  BED 

5450  H3-C3+S9+SS 

5460  !  CALCULATE  NEW  HEAT  FLUX  THRU  PLATE 

5470  X9=E1-H4-H3 

5480  V9»V9+1 

5490  !  COMPARE  NEW  PLATE  HEAT  FLUX  WITH  OLD  AND  REPEAT  CALCULATION  IF  DIFFERENCE 

EXCEEDS  0.01 
5500  X8-ABS(P3-X9) 

5510  IF  X8>=.01  THEN  P3-X9  6  GOTO  5340 
5520  !  PRINT  RESULTS 
5530  PRINT  USING  5540  :  P5 

5540  IMAGE  15X. "AVERAGE  PLATE  SURFACE  TEMPERATURE  «" . 5X . DDD. DD. X . "DEG  F" 
5550  PRINT  USING  5560  :  Bl 

5560  IMAGE  15X . "AVERAGE  BED  TEMPERATURE  «" , 15X . DDD. DD. X . "DEG  F" 
5570  PRINT  USING  5580  :  El 

5580  IMAGE  1 5X . "ELECTRICAL  ENERGY  INTO  HEATER  *" . 8X . DDD. DD. X , "BTU/HR" 
5590  PRINT  USING  5600  s  H4 

5600  IMAGE  15X."L0SS  FROM  HTR  BACKING  TO  ATM  «". 9X . DDD. DD. X , "BTU/HR" 
5610  PRINT  USING  5620  :  CI 

5620  IMAGE  15X."h(HTR  BACKING-ATM)  ■" , SX . DDD. DD. X , "BTU/HR-FT~2-DEG  F" 
5630  PRINT  USING  5640  ;  H3 

5640  IMAGE  15X."L0SS  FROM  HTR  PERIMETER  TO  BED  =". 7X . DDD. DD, X . "BTU/HR" 
5650  PRINT  USING  5660  :  P3 

5660  IMAGE  15X."HEAT  FLUX  THRU  PLATE  =", 17X , DDD. DD. X, "BTU/HR" 
5670  PRINT  USING  5680  ;  P7 

5680  IMAGE  15X, "h (PLATE-BED)  ■"  .  14X, DDD. DD, X, "BTU/HR-FT^2-DEG  F" 
5690  ! 

5700  !        CALCULATE  AND  DISPLAY  RESULTS  FOR  LEFT-HAND  HEATER 
5710  ! 

5720  PRINT  USING  5730 

5730  IMAGE  2/, 33X, "LEFT-HAND  HEATER" 
5740  !  CALCULATE  LH  HEATER  GUARD  TEMPERATURE  AVE 
5750  G2-(T(44)+T(45)+T(46) ) /3 

5760  !  CALCULATE  ELECTRICAL  ENERGY  INTO  HEATER 
5770  E2«V2*I2*3.4121 

5780  !  CALCULATE  LOSS  FROM  BACK  OF  HEATER  TO  ATMOSPHERE 
5790  G4-(H2-G2) *. 1703163 

5800  '  CALCULATE  CONVECTION  HEAT  TRANSFER  COEFFICIENT  FROM  GUARD  TO  ATMOSHPERE 
5810  C2=G4/ (.347222* (G2-A1) ) 

5820  !  CALCULATE  LOSS  FROM  SIDES  OF  HEATER  TO  ATMOSPHERE 
5830  S4«C2*.069444*(T(78)-A1> 
5840  !  CALCULATE  LOSS  FROM  TOP  STRIP  TO  ATMOSPHERE 
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5350  U3-C2*. 0208333* (T f"o -Al ) 

5S60  !  CALCULATE  LOSS  FROM  BOTTOM  STRIP  TO  ATMOSPHERE 

5870  S3-U3 

5330  !  CALCULATE  TOTAL  LOSSES  TO  ATMOSPHERE 

59«>0  H&-G4+54-HJ3+SB 

5900  !  CALCULATE  HEAT  FLUX  THRU  LH  PLATE 

5910  P4=E2-H6 

5920  '  CALCULATE  SURFACE  TEMPERATURE  OF  PLATE 

5930  P6-H2-P4* . 0208333/ ( (232-. 032* (H2-70) > *. 347222) 

5*40  !  CALCULATE  HEAT  TRANSFER  COEFFICIENT  FOR  PLATE  TO  BED 

5950  PB-P4/ t . 347222* (P6-B1 ) ) 

5960  !  CALCULATE  LOSS  FROM  SIDE  STRIPS  TO  BED 

5970  C4=*PB*.  0694441 (T(79) -Bl) 

5990  '  CALCULATE  LOSS  FROM  TOP  STRIP  TO  BED 

5990  U1-P8*. 0208333* (T(77)-B1) 

6000  !  CALCULATE  LOSS  FROM  BOTTOM  STRIF  TO  BED 

6010  S6=U1 

6020  '  CALCULATE  TOTAL  LOSSES  TO  BED 

6030  H5-C4+UH-S6 

6040  !  CALCULATE  NEW  HEAT  FLUX  THRU  PLATE 

6050  X7-E2-H6-H5 

6060  !  COMPARE  NEW  HEAT  FLUX  TO  OLD  AND  REPEAT  CALCULATION  IF  DIFFERENCE  EXCEEDS 

0.01 
6070  X6-ABS<P4-X7) 

6080  IF  X6>=.01  THEN  P4-X7  &    GOTO  5920 
6O90  !  PRINT  RESULTS 
6100  PRINT  USING  6110  :  H2 

6110  IMAGE  2/. 15X, "AVERAGE  HEATER  TEMPERATURE  -" . 12X, DDD. DD, X , "DEG  F" 
6120  PRINT  USING  6130  ;  P6 

6130  IMAGE  15X, "AVERAGE  PLATE  SURFACE  TEMPERATURE  ■" . 5X, DDD. DD, X. "DEG  F" 
6140  PRINT  USING  6150  :  Bl 

6150  IMAGE  15X. "AVERAGE  BED  TEMPERATURE  -" . 15X, DDD. DD. X . "DEG  G" 
6160  PRINT  USING  6170  :  E2 

6170  IMAGE  1 5 X, "ELECTRICAL  ENERGY  INTO  HEATER  ■" , 8X, DDD. DD. X . "BTU/HR" 
6180  PRINT  USING  6190  j  H6 

6190  IMAGE  15X,"L0SS  FROM  HTR  BACKING  TO  ATM  ■". 9X , DDD. DD. X , "BTU/HR" 
6200  PRINT  USING  6210  :  C2 

6210  IMAGE  15X,"h(HTR  BACKING-ATM)  ■" , 8X. DDD. DD, X, "BTU/HR-FT~2-DEG  F" 
6220  PRINT  USING  6230  5  H5 

6230  IMAGE  15X."L0SS  FROM  HTR  PERIMETER  TO  BED  -", 7X, DDD. DD, X. "BTU/HR" 
6240  PRINT  USING  6250  ;  P4 

6250  IMAGE  15X. "HEAT  FLUX  THRU  PLATE  ■", 17X. DDD. DD, X. "BTU/HR" 
6260  PRINT  USING  6270  :  PB 

6270  IMAGE  15X , "h (PLATE-BED)  ■" , 14X . DDD. DD, X , "BTU/HR-FT~2-DEG  F" 
6280  ! 

6290  !  CALCULATE  ENERGY  BALANCE  VALUES 
6300  ! 

6310  !  CALCULATE  AIR  DENSITY 
6320  X5=1.325*29.92/<P2+459.69) 

6330  I  CALCULATE  ENERGY  CARRIED  OUT  OF  BED  BY  AIR  STREAM 
6340  A5-F2*X5*60*.241* (P2-P1) 

6350  !  CALCULATE  LOSS  THRU  FRONT  AND  BACK  FACES 
6360  L6«.448*X*(T(75)-T<74> ) - 
6370  !  CALCULATE  TOTAL  ENERGY  LEAVING  BED 
6380  A6-A5+L6 

6390  !  CALCULATE  TOTAL  ENERGY  INTO  BED  FROM  RH  HEATER 
6400  B3-E1-H4 

6410  !  CALCULATE  TOTAL  ENERGY  INTO  BED  FROM  LH  HEATER 
6420  B4-E2-H6 
6430  !  CALCULATE  TOTAL  ENERGY  INTO  BED  FROM  ALL  SOURCES 
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6440  A7=B3*B4 

&450  !  CALCULATE  CONVECTION  HEAT  TRANSFER  COEFFICIENT  F/B  FACES-ATMOSPHERE 

6460  :>o=Lo,  ;24xX* (T(74) -Al) ) 

64"?0  !  CALCULATE  CONVECTION  HEAT  TRANSFER  COEFFICIENT  F/E  FACES-BED 

=43v  D7=L6>  C24*X* (Bl-T(75) ) ) 

64^0  !  CALCULATE  SUPERFICIAL  VELOCITY 

6500  U=F2/  '.  o*  X / 1 44 )  /  60 

6510  !  CALCULATE  SUPERFICIAL  MASS  VELOCITY 

6520  G=U*X5*3600 

6530  !  CALCULATE  PARTICLE  REYNOLDS  NUMBER 

6540  R=G/3600*. 0009166667/. 000000396/32. 174 

6550  !  PRINT  RESULTS 

6560  PRINT  USING  6570 

6570  IMAGE  2/ , 34X , "ENERGY  BALANCE" 

6580  PRINT  USING  6590  ;  PI 

6590  IMAGE  2/ , 15X , "AVERAGE  AIR  INLET  TEMPERATURE  -" , 9X , DDD. DD , X , "DEG  F" 

6600  PRINT  USING  6610  ;  P2 

6610  IMAGE  15X."AIR  OUTLET  TEMPERATURE  =" , 1 6X . DDD. DD . X , "DEG  F" 

6620  PRINT  USING  6630  ;  Bl 

6630  IMAGE  15X, "AVERAGE  BED  TEMPERATURE  ■" , 1 5X . DDD. DD, X , "DEG  F" 

6640  PRINT  USING  6650  :  A5 

6650  IMAGE  15X,"q(AIR  OUT)  =", 27X , DDD. DD, X , "BTU/HR" 

6660  PRINT  USING  6670  :  L6 

6670  IMAGE  15X,"q(L0SS  F/R  WALL)  =", 21 X , DDD. DD, X, "BTU/HR" 

6680  PRINT  USING  6690  ;  A6 

6690  IMAGE  15X,"q (TOTAL  OUT  OF  BED)  =", 18X , DDD. DD, X , "BTU/HR" 

6700  PRINT  USING  6710  ;  A7 

6710  IMAGE  15X."q (TOTAL  INTO  BED)  =", 20X , DDD. DD. X . "BTU/HR" 

6720  PRINT  USING  6730  ;  U 

6730  IMAGE  15X , "SUPERFICIAL  VELOCITY  -", 17X, DDD. DD, X , "FT/SEC" 

6740  PRINT  USING  6750  ;  G 

6750  IMAGE  15X , "SUPERFICIAL  MASS  VELOCITY  ■" , 7X, DDD. DD, X , "LBM/HR-FT~2" 

6760  PRINT  USING  6770  ?  R 

6770  IMAGE  15X ,  "PARTICLE  REYNOLDS  NUMBER  «"  ,  20X  .  DDD.  DD 

6780  END 
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